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The  purpose  of  this  study  was  to  design  longitudinal 
control  laws  for  the  analog  backup  flight  control  system  on 
the  X-29A  aircraft.  The  3 5 X  static  instability  encountered 
in  this  aircraft  suggested  a  unique  problem. 

A  stable  control  law  design  which  meets  Level  I  flying 
qualities  was  made  with  a  rigid  body  and  linear  equations 
of  motions  assumed.  The  Appendices  include  computer  output 
from  TTYLON  and  TOTAL,  the  computer  programs  used  to  find 
transfer  functions  and  do  root  locus  design,  respectively. 
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Abstract 


V 

X-29A  wind  tunnel  data  were  input  to  the  computer  pro¬ 
gram  TTYLON  to  produce  a/5c,  NZ/5C,  and  q/5c  bare  airframe 
transfer  functions  throughout  the  flight  envelope.  The 
short  period  mode  below  Mach  1  was  unstable.  Where  it  was 
stable,  the  response  was  slow. 

The  normal  acceleration  sensor  was  placed  so  that  the 
short  period  frequency  would  be  within  the  desired  range  of 
3  to  10  rad/sec.  Nz  feedback  was  used  to  stabilize  the 
short  period  mode.  Pitch  rate  feedback  was  also  used  to  add 
damping  for  Level  I  flying  qualities  in  the  short  period 
mode.  To  stabilize  the  phugoid  mode  and  achieve  Level  I 
flying  qualities,  a  lag  filter  design  was  developed.  Gain 
schedules  for  KNz  and  Kq,  which  were  found  using  classical 
multiloop  control  law  design  are  presented. 


Introduction 


I . 


Problem  Statement 

The  purpose  of  this  thesis  is  to  use  classical  control 
techniques  to  design  backup  analog  control  laws  for  the  X-29A 
Forward  Swept  Wing  aircraft.  The  3 5 Si  negative  static  margin 
inherent  in  the  airframe  makes  this  a  unique  problem.  The 
design  and  analysis  in  this  study  are  limited  to  the  longi¬ 
tudinal  axis.  Therefore,  all  references  to  control  laws  ar.d 
system  response  will  include  only  the  longitudinal  modes. 

Background 

Wing  sweep  was  introduced  to  aircraft  when  it  was  dis¬ 
covered  that  the  critical  Mach  number  could  be  increased, 
thus  delaying  the  transonic  drag  rise.  Sweepba ck,  however 
has  some  disadvantages,  such  as:  tip  stall,  which  reduces 
controllability;  poor  lateral  control;  and  increased  induced 
drag  at  high  incidence  angles.  Forward  sweep  while  still 
reducing  critical  Mach  number,  promises  to  reduce  trim  drag 
on  the  aircraft  and  increase  maneuverability  since  stall  be¬ 
gins  at  the  root  of  the  wing. 

The  canard-forward  swept  wing  configuration  being  used 
by  Grumman  has  the  above  mentioned  35%  negative  static  margin 
for  subsonic  flight.  It  becomes  stable  in  supersonic  flight. 
A  triplex  digital  flight  control  system  with  an  analog  back¬ 
up  will  be  used  to  handle  the  instability.  This  allows  the 
aircraft  to  be  operational  after  a  failure  in  one  or  all  of 


the  digital  flight  control  computers.  Modern  control  theory, 
specifically  implicit  model  following  is  being  used  as  the 
main  digital  control  law  design  tool  by  Grumman.  In  this 
study,  as  with  the  original  design,  classical  frequency  do¬ 
main  methods  will  be  used  for  the  analog  backup  control  law 
design. 

General  Approach 

The  aircraft  aerodynamic  characteristics  found  in  wind 
tunnel  testing  were  used  in  the  computer  program  TTYLON  to 
determine  the  aircraft  longitudinal  axis  transfer  functions. 
The  computer  program  TOTAL  was  then  used  to  do  frequency  do¬ 
main  analysis  and  design.  Feedback  states  were  decided  upon 
and  gains  and/or  gain  schedules  calculated  for  the  final 
system  design. 

Sequence  of  Presentation 

The  theoretical  basis  for  the  transfer  function  calcula¬ 
tions,  root  locus  design  and  analysis,  multiloop  closure, 
and  gain  scheduling  are  presented  in  Chapter  II.  The  actual 
steps  taken  in  design  can  be  found  in  Chapter  III.  Chapter 
IV  contains  the  conclusions  drawn  and  recommendations  made 
for  the  study. 

Assumptions 

The  longitudinal  equations  of  motion  for  the  aircraft 
are  assumed  to  be  the  linear  perturbation  equations  about 
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steady  level  flight.  The  aircraft  is  also  assumed  to  be  a 
rigid  body.  These  assumptions  must  be  used  to  allow  use  of 
existing  tools  and  to  limit  the  scope  of  the  problem. 
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Theory  and  Development 


II . 


Deriving  the  Transfer  Functions 

Wind  tunnel  data  were  obtained  from  the  X-29A  Advanced 
Development  Program  Office  (ADPO) .  These  data  are  in  the 
form  of  stability  derivatives,  lift  and  drag  coefficients, 
weights  and  moments  of  inertia.  The  data  may  be  found  in 
Appendix  A.  The  computer  program  TTYLON  -  A  Longitudinal 
Aircraft  Transfer  Function  Program,  which  is  maintained  by 
ASD/ENFTC  at  Wngnt-Patterson  AFB ,  was  used  to  calculate 
transfer  functions  from  the  wind  tunnel  data  for  selected 
aircraft  configurations  and  flight  conditions.  The  flow 
diagram  for  the  program  is  shown  in  Appendix  3  and  the  dotted 
lines  indicate  the  options  which  were  used.  The  stability 
derivatives  were  found  for  a  trimmed  aircraft  at  the  desired 
flight  condition.  The  longitudinal  perturbation  state  equa¬ 
tions  then  were  solved  for  the  transfer  functions.  The 
equations  used  are: 
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N  ,  the  normal  acceleration,  is  given  by  the  following 
equation: 


Nz  =  Uo  q/g  +  q/g  ~  Ur,  5/5  "  esinef 


As  can  be  seen  in  Fig  1.  the  X-29A  uses  three  types  of 
surfaces  for  longitudinal  trim  and  control. 


input.  This  assumption  is  valid  because  it  is  the  only 
active  control  surface;  the  others  are  passive.  Whenever 
an  elevator  input  is  referred  to  in  TTYLON,  the  canard  should 
be  substituted. 

The  flight  conditions  are  defined  by  Mach,  altitude, 
weight,  center  of  gravity,  and  flight  path  angle  m  TTYLON. 
For  this  study,  the  following  were  held  constant: 

1)  Weight  =  15,000  lb 

2)  Center  of  Gravity  =  451  in  (37.58  ft) 

3)  Flight  Path  Angle  =  0  degrees 

The  flight  conditions  chosen  for  control  law  design  are  shown 
in  Table  I. 


Table  I 

Flight  Conditions  Used  for  Control  System  Design 


Altitude  (ft)  Mach 

0  0.3 

0.6 
0.8 

20K  0.4 

0.6 
1.0 
1.2 

40K  0.5 


0.6 

1.0 

1.4 
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This  flight  envelope  is  limited  at  low  Mach  by  aerodynamics 
(stall) ,  and  at  high  Mach  by  structural  limits  (dynamic 
pressure) .  The  maximum  dynamic  pressure  is  1200  lb/ft  . 

The  roots  of  the  bare  airframe  canard  transfer  functions  at 
the  various  flight  conditions  are  shown  in  Table  II. 

Table  II 

Characteristic  Roots  of  Bare  Airframe 
Canard  Transfer  Functions 


Flight 

Condition 

Short 

Period 

Phugoid 

M 

h(ft) 

real 

imag 

real 

imag 

0.3 

0 

1.9478 

0 

-0 . 14311e-l 

0.14132 

-3.0539 

0 

0.6 

0 

4.0251 

0 

-0. 10817e-l 

0 . 69910e-l 

-6.6276 

0 

0.8 

0 

5.4774 

0 

-0 . 16487e-l 

0. 41562e-l 

-9.3031 

0 

0.4 

20K 

1.8793 

0 

-0.12824e-l 

0.11286 

-1.9912 

0 

0.6 

20K 

2.8675 

0 

-0.70548e-2 

0.76l64e-l 

-4.0591 

0 

1.0 

2QK 

3.9644 

0 

-0. 39750e-l 

0 

-6.0004 

0 

0.16974e-l 

0 

1.2 

2QK 

3.4147 

0 

-0 . 91117e-l 

0 

-5.7362 

0 

0 . 53264e-l 

0 

0.5 

40K 

1.5780 

0 

-0 . 15728e-l 

0 . 94651e-l 

-1.9912 

0 

0.6 

40K 

1.9091 

0 

-0 . 10247e-l 

0 . 80888e-l 

-2.4136 

0 

1.0 

40K 

2.6941 

0 

-0.67527e-2 

0 . 47972e-l 

-3.5394 

0 

1.4 

40K 

1.5136 

0 

-0. 82773e-2 

0. 39491e-l 

-2.5960 

0 
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Once  the  bare  airframe  transfer  functions  were  found, 
a  test  was  made  to  find  what  effect  the  canard  actuator 
dynamics  had  upon  the  bare  airframe  dynamics.  The  transfer 
function  for  the  canard  actuator  is  given  by: 


,  . _ (20.202)  (144.928)  (71.4) _ 

'a'  1  "  (s+20.202)  (s+144.928)  (s+52 . 55+j  48 . 34) 


Bode  plots  of  the  open  loop  N2/6c  transfer  function 
with  and  without  the  actuator  included  are  shown  in  Figs  2 
and  3.  One  can  see  that  the  response  is  affected  very  little 
and  only  above  40  rad/sec  which  is  not  in  the  frequency  range 
of  the  longitudinal  mode  dynamics.  The  effect  the  actuator 
had  on  dynamics  was  therefore  considered  negligible  and  the 
transfer  function  was  assumed  to  be  1. 

The  sensor  and  pilot  command  transfer  functions  were 
also  assumed  to  have  little  effect  on  the  aircraft  response. 
This  is  because  the  lag  in  both  cases  is,  in  general,  an 
order  of  magnitude  smaller  than  the  period  of  the  short  per¬ 
iod  mode. 

Root  Locus  Design 

The  computer  program  TOTAL,  which  is  maintained  by  Mr. 
John  Smith,  ASD/ENFTC ,  was  the  main  tool  in  this  design  and 
analysis  of  the  control  laws.  Once  the  bare  airframe  trans¬ 
fer  functions  were  found  using  TTYLON,  they  were  put  in 
TOTAL.  Option  41  produces  a  listing  of  the  root  locus. 
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Fig  3.  Nz/<Sc  with  Actuator  Dynamics 
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Looking  at  an  example  m  Appendix  E,  one  can  see  that  the 
gain  and  damping  ratio  at  each  point  on  the  locus  is  listed. 
By  looking  at  these  listings  for  a  series  of  flight  condi¬ 
tions,  one  may  find  the  gain  at  which  damping  ratio,  damped 
frequency,  or  stability  criteria  may  be  met.  Option  42  in 
TOTAL  allows  the  designer  to  obtain  the  closed  loop  roots 
for  the  system  with  a  specified  gain  on  the  state  being  fed 
back.  Option  43  allows  the  same  kind  of  root  finding  for  a 
specified  damping  ratio.  Examples  of  Option  42  and  43  list¬ 
ings  may  be  found  in  Appendix  E. 

Choosing  Outputs  to  Feed  Back 

Most  textbooks  which  address  feedback  control  law  de¬ 
sign  recommend  that  a  should  be  fed  back  to  stabilize  the 
short  period  mode  as  This  kind  of  stability  deri¬ 

vative  augmentation  causes  an  increase  in  the  short  period 
damping  ratio  and  undamped  natural  frequency.  The  phugoid 
mode  is  essentially  not  affected. 

For  this  case,  angle  of  attack  feedback  was  attempted 
first.  The  X-29A  ADPO  advised  against  this,  saying  that  the 
a  sensor  was  too  slow  for  a  fighter  type  aircraft.  In  cases 
where  the  a  sensor  lags  too  much  to  be  used  for  accurate 

feedback,  N  feedback  is  often  substituted.  Usually,  a  and 
z 

Nz  have  very  similar  root  loci  and  frequency  responses  within 
the  typical  range  of  aircraft  control  system  frequencies. 

In  this  case,  they  appear  to  be  "opposites."  Opposite  signs 
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on  the  gain  are  necessary  for  stability. 

In  Fig  4  one  can  see  that  in  the  180  degree  root  locus 

for  a/6,  the  locus  of  short  period  roots  begins  at  the  phu- 

goid  poles  and  goes  to  infinity,  while  in  Fig  5  the  N  /5  0 

z  c 

angle  locus  approaches  its  short  period  zeros  from  infinity. 

Because  the  N  accelerometer  is  more  accurate  for  feedback 
z 

than  the  a  sensor,  and  since  N  feedback  will  stabilize  the 

z 

short  period  mode  as  well  as  a  feedback,  Nz  was  chosen  for 
feedback  to  stabilize  the  short  period  mode. 

The  standard  method  of  increasing  damping  in  the  short 
period  mode  is  q  feedback,  and  that  is  the  method  applied 
here.  Once  the  Nz  loop  has  been  closed  and  the  short  period 
roots  are  stable,  one  can  see  from  Fig  6  that  feeding  back 
q  will  increase  the  damping  ratio  while  decreasing  the  damped 
frequency  of  the  short  period  mode. 

Placement  of  the  Normal  Acceleration  Sensor 

Once  it  was  decided  that  normal  acceleration  would  be 
fed  back  to  stabilize  the  short  period  mode,  the  acceler¬ 
ometer  location  had  to  be  decided  upon.  Looking  at  the  Nz/Sc 
transfer  functions  that  result  from  locations  between  the 
nose  and  the  tail,  it  was  found  that  the  zeros  of  the  short 
period  were  most  affected  by  sensor  placement. 

After  looking  through  Flying  Qualities  Handbooks  and 
data  from  aircraft  with  similar  purposes,  and  considering 
past  experience,  a  desired  damped  frequency  of  about  4  rad/sec 
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ANGLE  OF  ATTACK  AT  M=0.6.  H=20K  FT 


Id  MZ-H  *9 ’0= 


for  the  short  period  was  chosen.  Knowing  that  the  transfer 
function  for  the  short  period  approximation  is: 


Z-  -X  (M-  +Z-  M-)][S2  +  2;u)  s+gj  2) 
o  a  a  o  w  n  n 

c _ c  c _ 

2  2 
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sp  n  n 

sp  sp 


(4) 


One  can  see  that: 
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(10:450) 


The  damped  frequency  of  the  numerator  zeros  increase  as  X  , 

a 

the  distance  from  the  center  of  gravity  to  the  sensor  loca¬ 
tion,  increases.  One  can  see  from  the  actual  root  loci  of 
the  transfer  function  in  Appendix  C,  that  the  further  forward 
that  the  normal  acceleration  sensor  is  placed  with  respect 
to  the  center  of  gravity,  the  higher  the  damped  frequency  of 
the  short  period  open  loop  zeros.  By  looking  at  Figs  5  and 
6  one  can  see  that  closure  of  the  N2  loop  will  stabilize  the 
short  period  roots  and  make  them  oscillatory,  while  closure 
of  the  q  loop  will  decrease  frequency  and  increase  damping. 
Therefore,  the  location  of  the  accelerometer  was  chosen  so 

that  the  N  /5  zeros  for  the  short  period  mode  were  somewhat 
z  c 

higher  than  the  desired  short  period  root  frequency  of  4 
cycles/'sec.  With  the  Nz  sensor  located  at  the  20  ft  station, 
the  short  period  zeros  were  between  2.22  and  8.90  rad/sec. 

The  TTYLON  runs  for  Mach  =  0.6  at  20,000  ft  altitude  are 
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shown  m  Appendix  D.  With  the  center  of  gravity  at  37.53  ft 
aft  of  the  nose,  the  accelerometer  was  positioned  20  ft  aft 
of  the  nose  as  it  gives  the  short  period  mode  a  damped  fre¬ 
quency  of  4.2318  rad/sec  before  the  loop  closure. 

Gain  Scheduling 

For  the  desired  damping  ratio  of  the  closed  loop  roots, 
0.7  was  chosen  because  it  provides  a  fast  response  with  a 
reasonable  amount  of  overshoot.  More  damping  will  slow  the 
system  down,  and  less  will  increase  the  overshoot.  Gain 
scheduling  was  necessary  in  order  to  achieve  the  desired 
damping  ratio.  Looking  through  the  root  locus  for  each 
flight  condition,  the  gain  which  resulted  in  a  0.7  damping 
ratio  for  the  closed  loop  short  period  roots  was  determined. 
These  gains  were  plotted  versus  Q  at  each  altitude  and  Mach 
number.  A  relationship  between  q  and  can  then  be  sought 
for  the  gain  schedule. 

The  method  used  to  find  this  relationship  is  the  least 
squares  fitting  of  data.  Beginning  with  the  basic  system 
equation: 

A  x  =  b,  (6) 

T 

one  may  premultiply  both  sides  by  A  and  the  equation  is: 

AT  A  x  =  AT  b.  (7) 

In  this  particular  case,  x  is  the  vector  [c 

lD 
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where  C  and  D  are  the  coefficients  m  the  equation 


Kq  =  C  +  D  Q 

The  vector  x  may  be  found  by  using 

T  - 1  T 
x  =  (AXA)  1  Ad 


where 


and 


A  = 


(  3 ) 


(?) 


(10) 


(il) 


This  method  will  give  the  equation  for  the  best  straight  line 
with  the  least  amount  of  error  squared  between  the  data  points 
and  the  chosen  line. 

For  the  second  order  curve  fit,  one  is  solving  for  the 


vector  x  = 


equation : 


where  C,  D,  and  E  are  the  coefficients  m  the 


K  =  C  +  DQ  +  EQ^ 

q 


(12) 
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Eq  9  is  again  used  to  solve  for  x  with  the  sar.e  b  matrix  as 
in  Eq  11,  but  with  the  following  A  matrix: 


1  -i  V 
1  Q2  V 


1  '“M  V 


113} 


This  equation  will  result  in  the  Dest  paraDOla,  hyperbola, 
or  circle  to  fit  the  data  points. 


Multilooo  Closure 

A  method  of  multiloop  control  law  design  for  the  fre¬ 
quency  domain  is  described  by  McGruer,  Ashkenas,  and  Graham 
in  Ref  10.  This  method  allows  one  to  achieve  with  equations 
the  same  results  that  were  found  using  the  successive  loop 
closure  method  discussed  in  the  section,  "Choosing  Outputs 
to  Feed  Back." 

Coupling  numerators  are  necessary  for  the  multiloop 
technique  when  there  is  more  than  one  input.  For  this  study, 
there  is  only  the  canard  input  to  the  longitudinal  mode 
response,  therefore  coupling  numerators  will  not  be  discussed. 
A  thorough  explanation  may  be  found  m  Ref  10.  Fig  7  illus¬ 
trates  a  block  diagram  of  a  flight  control  system's  longitu¬ 
dinal  axis  with  simultaneous  Nz  and  q  feedback  through  the 

canard.  To  derive  the  transfer  functions  for  N  /5  and  q/5  , 

z  c  c 

one  must  begin  with  the  following  equation: 
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Nz  (s) 


L9 


£  (s) 


:  (s  )  G  (s  ) 

P  c 


Is  >  N 


z 


5  ) 


K  G  ( s  )  q  ( s  )  f 14) 

q  q 


where : 


N'z(s)  =  [Nz  (s) /5c(s)  ]0LE(s)Gs  (s)  (15) 


and 


q ( s )  =  [q (s) /5c (s)  ]QLE (s)Gs  (s) 


(16) 


Substituting  Eqs  15  and  16  into  Eq  14  yields: 

E(s) {1+Kn  Gn  (s)Gs(s) [Nz (s)/oc (s) ]ql 
+  KqGq  (s ) Gg (s )  [q  (s )  /  *c  (s )  ] =  t  ( s )  G^  ( s  )  (1™) 


The  command  and  servo  transfer  functions,  Gc(s)  and  Gg(s), 

have  been  assumed  to  be  1  as  was  explained  in  the  section 

"Deriving  the  Transfer  Functions."  This  means  that  E(s)  = 

■5  (s)  =  <5  (s)  .  The  two  multiloop  transfer  functions  then 
P  <- 

simplify  to: 


[N  (s) /5  (s)  ] 
z  c  CL 


[Nz(5)/5c(s)]ol 


1+kN  [NZ(S>/6C(S)]0L 
z 


+Kq[q(s)/5c{s)  ]0L 


(18) 


and 


[q(s)/5c(s)  1CL 


[q  (s)/5c(s) ]QL 

1+KN  Tnz(s)/5c(s) ]0L+Kq[q(s)/5c(s) ]QL 

2  (19) 
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III.  Design 


Multiloop  Design 


The  gain  chosen  for  Nz  was  fed  back  and  the  Nz  closed 
loop  roots  were  found  using  Option  42  of  TOTAL.  From  Fig  7, 


one  can  see  that  when  the  Nz  loop  is  closed,  this  transfer 
function  is: 


nJ  _  NlP2 

6  c  =  +  DxD2 

Jrr. 


_ ^2num _ 

l  z  ^num 


(20) 


_ (As4-rBs^+Cs -^-t-Ds  +  E) _ 

Kn  (As4+Bs^+Cs^+Ds+E) + (s4+as^+bs^+cs+d) 


The  highest  order  coefficient  in  the  denominator  becomes 

(K”  A+l) .  The  effective  q/o  open  loop  transfer  function 
z  *" 

is  then: 


<3  _  _ 3num _ 

<5C  (Km  A+l)  ( N 2  closed  loop  roots) 

OL 


(22) 


A  gain  Kq  is  then  found  which  meets  the  desired  damping  ratio 
requirements.  After  finding  the  q  closed  loop  roots  with 
Option  42,  the  closed  loop  transfer  function  is  found: 
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g 

_ ^nura _ 

(K^j  A+l )  (q  closed  loop  roots) 


(25) 


This  is  the  overall  closed  loop  q/6c  transfer  function.  The 

overall  N  /o  transfer  function  of  course  has  the  same  char- 
^  c 

actenstic  equation  and  the  original  Nz  numerator: 
f  n  N2 

2  =  _ ^num _  (26) 

&c  (K^  A+l) (q  closed  loop  roots) 

L  JCL  ‘2 


The  Nz  loop  was  closed  first  because  this  feedback  stabili2ed 
the  short  period  mode.  Any  kind  of  special  flight  control 
modes  such  as  altitude  hold,  or  speed  control  which  may  be 
desired  can  be  wrapped  around  the  Nz  closed  loop  now  that  the 
system  is  stable.  In  this  case  the  pitch  rate  feedback  was 
wrapped  around  the  Nz  feedback  to  make  the  system  meet  Level 
I  Flying  Qualities. 

Selection  of  Gains 

An  attempt  was  made  to  find  a  K^z  gain  schedule  which 
would  result  in  a  constant  damped  frequency  of  the  short 
period  roots.  The  highest  frequency  for  which  this  is 
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possible  with  a  stable  short  period  .-node  is  4.7  rad/sec. 

Mach  =0.5,  h  =  40,000  ft  is  the  limiting  flight  condition 

at  which  the  short  period  is  unstable  at  damped  frequencies 

higher  than  4.7  rad/sec.  Fig  8  illustrates  this.  At  all 

the  higher  Mach  flight  conditions,  the  open  loop  NZ/<5C  toot 

locus  is  limited  by  its  short  period  zero  locations  and  will 

not  go  to  a  damped  frequency  as  low  as  4.7  rad/sec.  At 

Mach  =  0.8,  sealevel,  the  lowest  frequency  attainable  is  6.9 

rad/sec.  Fig  9  illustrates  this  problem.  As  was  discussed 

in  the  section,  "Placement  of  Normal  Acceleration  Sensor, " 

the  N  accelerometer  location  is  the  factor  which  drives 
z 

this  frequency  flexibility. 

Since  constant  damped  frequency  was  not  possible, 
another  simple  approach  to  Nz  feedback  was  attempted.  It 
was  decided  that  a  constant  gain,  Kjg  ,  would  be  found  to 
stabilize  the  short  period  mode  for  all  flight  conditions 
and  that  the  short  period  damped  frequency  would  be  allowed 
to  vary.  K^z  =  -1.2  was  chosen  as  it  was  the  lowest  gain 
which  would  stabilize  the  limiting  condition  Mach  =  0.5,  h  = 
40,000  ft.  This  can  be  seen  in  the  root  loci  in  Appendix  E. 
With  the  Nz  loop  closed  and  the  accelerometer  at  the  20  ft 
station,  the  damped  frequencies  varied  from  4.01  rad/sec  to 
9.21  rad/sec.  These  results  of  the  Nz  loop  closure  are 
shown  in  Table  III. 
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OPEN  LOOP  NZ  FOR  M=0.5.  H=40K 


Fig  9.  N2/oc  at  M=0.8,  Sealevel 


Table  III 


Roots  Resulting  From  K^z  =  -1.2 
Flight  Condition  Short  Period  Phugoid 


M 

h  (ft) 

real 

imag 

real 

imag 

0.3 

0 

-0.29209 

4.0147 

0. 43183e-2 

0 . 60757e-l 

0.6 

0 

-0.82973 

6.7733 

-0 . 98176e-2 

0 . 91793e-2 

0.8 

0 

-1.2305 

9 . 2053 

-0 . 33280e-l 

0 

-0. 64486e-3 

0 

0.4 

20K 

-0.17979 

4.0600 

0 . 46977e- 2 

-0.5  2  39  5e-l 

0.6 

20K 

-0.37512 

4.9676 

-0 . 19634e-2 

-0 . 24173e-l 

1.0 

20K 

-0.56668 

6.5666 

-0 . 24510e-l 

0 

0 . 18225e-2 

0 

1.2 

20K 

-0.56932 

6.6404 

-0 . 34739e-l 

0 

0. 32160e-2 

0 

0.5 

40K 

-0 . 11778e-l 

4.7494 

-0 . 54698e-2 

-0. 51752e-l 

0.6 

40K 

-0.11849 

4.2240 

0 . 40106e-2 

0 . 38340e-l 

1.0 

40K 

-0.23649 

4.7065 

-0. 37698e-2 

0 . 12623e-l 

1.4 

40K 

-.26632 

4.9925 

-0 . 1 559 5e-l 

0 

-0. 36254e-3 

0 

One  can  see  that  the  short  period  mode  has  now  been  stabil¬ 
ized.  The  damping  ratios  range  from  C.002  to  0.13.  This  is 
much  too  low,  and  an  increase  in  damping  must  be  achieved. 

As  was  discussed  in  the  section,  "Gain  Scheduling, "  K 
was  chosen  to  try  to  achieve  6  =  0.7  for  all  short  period 
roots.  Table  IV  shows  the  Q  vs  schedule  that  provided 
the  0.7  damping  ratio  desired. 
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able  IV 


K 

q 

to  Achieve  ;  =  0.7  at 

Selected  Flight 

Condition. 

Flight  Condition 

Kq 

M 

h  (ft) 

Q(lb/ft2) 

0.5 

40K 

69 

2.08 

0.6 

40K 

99 

3.00 

0.4 

20K 

109 

3.18 

0.3 

0 

134 

3.65 

0.6 

20K 

246 

5.25 

1.0 

40K 

275 

5.98 

0.6 

0 

535 

8.55 

1.4 

40K 

539 

10.25 

1.0 

20K 

682 

10.08 

0.8 

0 

951 

12.21 

1.2 

20K 

982 

13.45 

Looking  at  the  data  points  plotted  m  Fig  10,  one  can 
see  that  there  are  a  few  curves  which  look  like  they  may  fit 
the  points.  In  this  study,  three  methods  were  tried.  First, 
a  linear  curve  fit  using  all  the  data  points  was  used.  Sec¬ 
ond,  a  quadratic  curve  fit  was  tried,  and  finally,  the  points 
were  broken  into  two  groups.  One  line  was  fit  through  the 
high  Q  group  and  another  through  the  low  Q  group.  These 
methods  shall  be  referred  to  as  Method  I,  II,  and  III,  respec¬ 
tively.  Using  the  method  of  least  squares  described  m  the 
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urve 


section,  "Gain  Scheduling,"  the  vector  of  coefficients 


\c 

or  x  = 

'C' 

Id 

D 

was  solved  for.  The  work  can  be  found  m 


Appendix  F.  In  Method  I,  the  equation  for  the  line  which 
results  is: 


K  =  2.081  +  0.01142Q  (27) 

Using  a  second  order  curve  fit  in  Method  II,  the  equation 
for  Kq  is: 

K  =  1.14  +  0.01866Q  -  0 . 000006849Q2  (28) 

4 

Grouping  the  data  points  in  Fig  10  into  the  six  lower  and 
the  five  higher  points,  the  two  lines  which  result  from 
Method  III  are: 

2 

For  Q  less  than  or  equal  to  405  lb/ft 


Kq  =  1.174  +  0.01717Q  (29) 


For  Q  greater  than  405  lb/ft^ 


Kq  =  4.839  +  0 . 008226Q  (30) 

The  damping  ratios  for  the  closed  loop  roots  using  the  gains 
found  by  each  of  these  methods  are  shown  in  Table  V. 
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Table  V 


Damping  Ratios  for  Closed  Loop  Roots  Resulting 
From  the  Three  Curve  Fitting  Methods 


Flight 

Condition 

M 

h 

Q (lb/ft2) 

Method  I 

Method  II 

Method  III 

0.3 

0 

134 

0.69 

0.72 

0.67 

0.6 

0 

535 

0.68 

0.75 

0.74 

CO 

• 

o 

0 

951 

0.73 

0.73 

0.72 

0.4 

20K 

109 

0.73 

0.73 

0.68 

0.6 

20K 

246 

0.66 

0.74 

0.72 

1.0 

20K 

682 

0.69 

0.73 

0.72 

1.2 

20K 

982 

0.69 

0.69 

0.68 

in 

• 

o 

40K 

69 

0.95 

0.88 

0.79 

0.6 

40K 

99 

0.74 

0.74 

0.67 

1.0 

40K 

275 

0.62 

0.70 

0.69 

1.4 

40K 

539 

0.58 

0.65 

0.64 

The  desired  damping  ratio  is  0-7.  Below  are  listed  the  mean 
and  standard  deviation  for  each  set  of  points. 

Table  VI 

Determination  of  Goodness  of  Fit  For 
Three  Curve  Fitting  Methods 

Method  I  Method  II  Method  III 

Mean  0.7054  0.7327  0.7018 

Standard  Deviation  0.0897842  0.0539509  0.0395006 
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One  can  see  that  Method  III  has  the  mean  which  is  closest 
to  0.7  and  has  the  least  standard  deviation  from  it.  There¬ 
fore,  the  best  linear  fit  to  the  K  curve  is: 

q 

2 

If  Q  is  less  than  or  equal  to  405  lb/ft 

Kq  =  1.174  +  0.01727Q  (29) 

And  if  Q  is  greater  than  405  lb/ft^ 

K  =  4.839  +  0 . 008226Q  (30) 

q 

These  equations  can  be  put  into  the  analog  computer  to 
control  the  gain  changes.  The  roots  that  result  from  this 
schedule  can  be  seen  m  Table  VII. 
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Flight  Phugoid  Short  Period 

Condition 


y 

h 

real 

imag 

real 

imag 

r 

0.3 

0 

0 . 269  56e- 2 

0 . 5 4711e- 1 

3.0201 

3.3035 

0.67 

0.6 

0 

-0. 10159e-l 

0 . 55097e- 2 

-5.8805 

5.3295 

0.74 

0.  s 

0 

-0. 3070 le- 1 

0 

-7.8115 

7.5555 

0.72 

-0. 51049e-3 

0 

0.4 

20K 

-0. 37562e-2 

0 . 48607e- 1 

-2.9694 

3.2268 

0.68 

0.6 

20K 

-0 . 24807e-2 

0 . 2l361e- 1 

-3.9429 

3.8222 

0.72 

1.0 

20K 

-0 . 24323e-l 

0 

-5.3280 

5.1358 

0 . 72 

0.14568e-2 

0 

1.2 

20K 

-0. 33995e-l 

0 

-5.1033 

5.5252 

0.68 

0 . 25804e-2 

0 

0 . 5 

40K 

0 . 50595e-2 

0. 49254e-l 

-3.9592 

3.0402 

0.79 

0.6 

40K 

0. 35963e-2 

0.36407e-l 

-3.0001 

3.2901 

0.67 

l.C 

40K 

~0. 39368e-2 

0. 11629e-l 

-3.5106 

3.6391 

0.69 

1.4 

40K 

-0 . 30  269e- 3 

0 

-3.4927 

4.1618 

0.64 

v 
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Table  VIII 

Shcrt-Term  Pitch  Response  to  Pitch  Controller 


(LjspTe2 

vS.  ;sp) 

Flight 

Condition 

wsp 

To 

^2 

U)spT02 

C 

Mach 

h  (ft) 

0 . 3 

- 

3.3035 

9 . 1889 

30 . 3555 

0.67 

0 . 6 

a 

VJ 

5.3295 

3 . 8709 

20.6298 

0 .74 

0.8 

A 

V 

7.5555 

2.5883 

19.5562 

0.72 

0.4 

20K 

3.2263 

13.378 

43 . 1686 

0.68 

0.6 

2  OK 

3.8222 

8.4619 

32.3429 

0.72 

1.0 

2  OK 

5.1358 

5.2871 

27.1535 

0.72 

1.2 

2  OK 

5.5252 

4.6868 

25.8958 

0.68 

0 . 5 

4  OK 

3.0402 

22.230 

67.5823 

0.79 

0.6 

4CK 

3.2901 

19.032 

62.6187 

0.67 

1.0 

*  A  \r 

3.6391 

12.306 

44.7836 

0.69 

1.4 

40K 

4.1618 

8.9739 

37.3477 

0.64 

Alien  these  values  are 

plotted 

as  in  Fig 

12,  one  can 

see  that 

the  short  period  perf 

ormar.ce  meets  Level 

I  standards 

• 

Phugoid 

Desicr.  Modification 

One  can  see  from  Table  VII  that  half  of  the  flight  con¬ 
ditions  with  the  Nz  loop  closed  have  an  unstable  phugoid 
mode.  This  occurred  when  was  fed  back  to  stabilize  the 

short  period  mode  because  the  Nz/:c  transfer  function  has 
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zeros  in  the  right:  half  plane.  At  that  time  it  was  decided 
that  the  stability  of  the  short  period  mode  was  of  greater 
importance  than  the  phugoid  mode  because  of  the  faster  re¬ 
sponse  of  the  former.  The  N2  feedback  which  resulted  in  the 
phugoid  instability  was  used.  The  Flying  Qualities  Specifi¬ 
cation  requires  the  time  to  double  amplitude  for  the  phugoid 
mode  be  greater  than  55  sec  for  Level  III.  Time  response 
plots  such  as  the  one  in  Fig  13  show  that  the  phugoid  mode 
time  to  double  amplitude  for  this  design  is  between  250  and 
300  sec.  Since  this  is  a  fighter  aircraft  and  the  pilot 
would  be  expected  :o  give  the  aircraft  a  new  input  at  least 
every  minute,  one  can  say  that  the  time  to  double  amplitude 
restriction  is  sufficient.  However,  to  have  Level  I  flying 
qualities,  the  phugoid  mode  must  be  stable  and  have  a  damp¬ 
ing  ratio  of  at  least  0.04.  To  meet  this  standard,  some 
kind  of  filter  must  be  added. 

Using  the  Mach  =  0.3,  sealevei  flight  condition  for  an 
example,  a  filter  has  been  designed  to  show  how  the  phugoid 
mode  can  be  stabilized.  The  gain,  Kq,  which,  according  to 
the  gain  schedule,  will  be  fed  back  at  this  flight  condition 
is  -3.483.  Looking  at  the  listing  of  the  root  locus  of  the 
uncompensated  system  with  the  Nz  loop  closed  and  the  q  loop 
open,  one  can  see  that  the  phugoid  mode  is  presently  unstable 
at  gains  greater  than  -13.0657.  To  stabilize  it,  a  lag  net¬ 
work  has  been  introduced  in  cascade  with  the  plant.  The  pole 
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of  the  filter  has  beer,  placed  close  to  the  on::r.  at  s  = 
-0.002  +  jO.  The  zero  of  the  filter  was  placed  to  the  left 
of  the  zero  at  -0.12822  +  jO  in  order  to  bring  the  locus 
left.  After  a  few  trials,  the  lag  filter  chosen  is  then: 


!s) 


=  s  +  1/T 
s  +  1/aT 


(31) 


where 

T  =  16.67  sec 
a  =  30 


so 


Gc(s) 


s  ■+■  0.06 
s  +  0.002 


(32) 


Now  one  can  see  from  Table  IX  that  the  phugoid  mode  is  stable 
below  a  gain,  Kq  =  -1.9427.  This  means  that  for  the  gain 
Kq  =  -3.483  which  will  be  fed  back  to  meet  short  period  mode 
requirement  for  x,  =  0.7,  the  phugoid  mode  will  also  be 
stable . 
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Table  IX 

Root  Locus  for  Phugoid  Mode  With  and 
Without  Lag  Filter 


Without  Lag  Filter 
Real  I nag  Gain 

0 . 43183e-2  0.60757e-l  0 

0.17287e-2  0.51098e-l  -6.2058 
0  0. 44441e-l  -13.066 

-0 . 22891e-2  0.34706e-l  -30.938 
-0 . 41647e-2  0.24884e-l  -72.668 
-0. 54938e-2  0.14973e-l  -205.25 
-0 . 61827e-2  0.49968e-2  -873.63 
0  0  0 


With  Lag  Filter 
Real  I  nag  Gam 

0 . 43183e-2  0.60757e-l  0 

0  0 . 57632e-l  -1.9417 

-O.14470e-1  0.43869e-l  -14.744 
-0 . 25280e-l  0.27091e-l  -4". 912 
-0 . 31051e-l  0 . 80022e-2  -122.40 
-0 . 2 1600e- 1  0  -133.89 

-0 . 51600e-l  0  -329.21 

-0 . 60000e- 1  0  0 


Option  42  in  TOTAL  was  used  again  to  determine  the  roots  of 
q/5c  with  the  lag  filter  at  the  desired  Kq  =  -3.483  for  this 
flight  condition.  These  roots  are: 

Filter  s  =  -0.0019251 

Short  Period  s  =  -3.0092  -  j3.3502 

Phugoid  s  =  -0.0027732  t  jO. 055429 

This  phugoid  mode  has  a  damping  ratio  of  0.04997,  which 

means  that  it  has  Level  I  flying  qualities.  The  time  re¬ 
sponse  of  the  longitudinal  mode  with  the  lag  filter  is  shown 
in  Fig  14.  The  root  loci  that  shew  the  phugoid  mode  with 
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and  without  the  filter  may  he  seen  in  Figs  15  and  16.  A 
design  for  only  one  flight  condition  is  shown  here.  If  it 
is  decided  that  the  phugoic  mode  must  be  stable  and  meet 
Level  I  flying  qualities  at  all  points  in  the  flight  envel¬ 
ope,  then  this  procedure  must  be  followed  for  all  flight 
conditions . 


nz  closed  phugoid  without  filler 


M-0.3.  SERLEVEL 

Ju.  t.it 


-4.W  •l.ll  <4.1*  4.1 


SiS-iS.  iJJ  I 


Fig  15.  Closed  Nz,  Open  q  Loop  Phugoid  Withou 
Filter  at  M=0.3,  Sealevel 


NZ 

CLOSEQ 

PHUGOID  WITH  FILTER  M=0.3. 

J—  «.  IS 

SERLEVEL 

a 

“  . . G 

•  '»  |.l( 

\-o-n 

0L  IK  (  S  I 


ms-o.qqti?o.p7  4iis-c.:si';?-c,7;'?s-8.^ci 


<S-0.  0 C2  1C  S'- 3.  C'?*5.<;..-Q4  >•  S'-?  S045-I*.  I 


Fig  16.  Closed  Nz,  Open  q  Loop  Phugoid  With 
Filter  at  M=0.3,  Sealevel 


M 


IV. 


Conclusions  ar.d  Recommendations 

This  study  has  shown  that  an  analog  backup  flight  con¬ 
trol  system  which  results  in  Level  I  handling  qualities  can 
be  designed.  Attention  must  be  paid  to  the  placement  of  the 
normal  acceleration  sensor  to  get  good  short  period  perfor¬ 
mance.  To  stabilize  and  improve  the  short  period  mode,  sim¬ 
ple  gain  feedback  of  the  Nz  and  q  states  were  used.  Normal 
acceleration  was  fed  back  first  to  stabilize  the  short  period 
mode  and  then  pitch  rate  feedback  was  added  to  improve  per- 
iormance  to  Level  I  Flying  Qualities  standards.  To  stabilize 
the  phugoid  mode  a  lag  filter  may  be  added. 

This  design  was  done  using  a  rigid  body  assumption. 
Further  study  could  be  done  by  looking  at  the  aeroelastic 
effects  of  the  longitudinal  axis.  Other  design  methods  may 
also  be  used  to  compare  complexity  and  performance.  A  design 
of  the  control  laws  for  the  lateral-directional  axis  should 
be  done  with  aeroelastic  effects  considered. 


41 


Bibliography 


1.  D'Azzo,  John  J.  and  Constantine  H.  Houpis.  Linear  Con¬ 
trol  System  Analysis  and  Design.  New  York:  McGraw-Hill 
Book  Company,  1981. 

2.  Grumman  Aerospace  Corporation.  Forward  Swept  Wing 
Demonstrator  Technology  Integration  and  Evaluation  Study. 
AFWAL-TR-80-3145 .  Vol  I;  Wright-Patterson  AFB,  Ohio: 
Flight  Dynamics  Laboratory,  June  1980. 

3.  Grumman  Aerospace  Corporation,  Forward  Swept  Wine 
Demonstrator  Technology  Integration  and  Evaluation  Studv. 
AFWAL-TR-80-3145.  Vol  II;  Wn ght-Pattersor.  AFB,  Ohio: 
Flight  Dynamics  Laboratory,  June  1980. 

4.  Heffley,  Bob  and  Wayne  Jewwell.  Aircraft  Handling  Quali¬ 
ties  Data.  NASA-CR-2144 .  Hawthorne,  California: 

Systems  Technology  Incorporated,  December  1972. 

5.  Hogg,  Robert  V.  and  Allen  T.  Craig.  Introduction  to 
Mathematical  Statistics.  New  York:  Macmillan  Publishing 
Co. ,  Inc. ,  1978 . 

6.  Hoh,  et  al.  Proposed  Mil  Standard  and  Handbook  -  Flvinc 

-  Qualities  of  Air  Vehicles  -  Vol  II:  Proposed  MIL  Hand- 

'  book .  AFWAL-TR-82- 3081 .  Wricht-Patterson  AFB,  Ohio: 

Flight  Dynamics  Laboratory,  July  82. 

7.  House,  Eric  B.  II.  Investigation  of  the  YF-16  in  High 
Angle  of  Attack  Asymmetric  Flight.  MS  Thesis.  AFIT/ 
GAE/AA/78M-6 .  Wright-Patterson  AFB,  Ohio:  School  of 
Engineering,  March  1978. 

8.  Larimer,  Stanley  J.  TOTAL  -  An  Interactive  Linear  Con¬ 
trol  Design  and  Analysis  Program.  ATIT/ENE,  Wright- 
Patterson  AFB,  OH. 

9.  McGlynn,  Henry  J.  TTYLON  -  A  Longitudinal  Aircraft 
Transfer  Function  Program.  ASD/ENFTC,  Wright-Patterson 
AFB,  OH. 

10.  McRuer,  Duane,  Irving  Ashkenas,  and  Dunstan  Graham. 
Aircraft  Dynamics  and  Automatic  Control.  Princeton: 
Princeton  University  Press,  1973. 

11.  Roskam,  Jan.  Airplane  Flight  Dynamics  and  Automatic 
Flight  Controls  -  Part  I.  Lawrence,  Kansas:  Roskam 
Aviation  and  Engineering  Corp.,  1979. 

V. 


42 


Appendix  A 
Wind  Tunnel  Data 

Wind  tunnel  data  in  the  form  of  weights,  geometry, 
moments  of  inertia,  stability  derivatives,  and  drag  data 
were  used  in  the  TTYLON  program  to  get  the  X-29A  longitu¬ 
dinal  transfer  functions.  This  data  follows. 


Appendix  B 

TTYLON  Flow  Diagrams 

The  computer  program  TTYLON,  which  is  maintained  by 
ASD/ENFTC  at  Wright-Patterson  AFB,  was  used  to  transform 
wind  tunnel  data  into  transfer  functions  to  model  the  X-29A. 
Several  flow  diagrams  follow  which  show  the  way  the  program 
works.  First,  an  overall,  top  level  diagram  is  shown.  Next, 
the  weight  and  moment  of  inertia  input  procedure  is  por¬ 
trayed.  The  last  two  flow  diagrams  show  the  method  by  which 
the  aerodynamic  data  is  input  and  the  way  in  which  a  trans¬ 
fer  function  for  each  desired  flight  condition  is  produced. 
The  dotted  lines  indicate  the  options  which  were  used  in 
this  design. 


Fig  B-l.  Overall  TTYLON  Flow  Diagram 
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Fig  B-2.  Parameter  Entry  Logic 


Fig  B-3.  Aero  Data  Entry  Sequence 


YES 


NO 


Aero  Data  Entry  Sequence  (Cont'a) 


ENTER:  | 
MACH 
ALTITUDE 
WEIGHT 
CG 

GAMMA 


LIST  TRIM 
DATA 


3LE  TC\N£_ 
TRIM?  / 


TO  E 


.1ST  STABILITY 
DERIVATIVES 


Short  period 

2PR0XIMATI0N2- 


ELEVATOR 

TF? 


YESJ  LIST  DEN 
> — •’CHARACTERISTICS 


SELECT 
TF  CODE 


TO  E 


THRUST  TF? 


STOP? 


LIST  NUM 
CHARACTERISTICS 


Fig  B-4.  Flight  Case  Data  Section 


Appendix  C 

TTYLON  Runs  for  Nz  Sensor  Location 


The  location  of  the  N2  accelerometer  will  determine  the 
position  of  the  open  loop  short  period  zeros  in  the  N2/£c 
transfer  function.  Lx,  the  distance  the  sensor  is  placed 
behind  the  nose  of  the  aircraft,  has  been  varied  from  0  ft 
to  50  ft.  (The  aircraft  is  48  ft  long) .  This  causes  the 
short  period  zeros  to  vary  from  a  frequency  of  2.8552  rad/ 
sec  to  9.1511  rad/sec.  The  actual  TTYLON  output  can  be  seen 
on  the  following  pages.  A  discussion  of  the  final  N2  accel¬ 
erometer  location  choice  may  be  found  in  the  section, 
"Placement  of  the  Normal  Acceleration  Sensor." 


DENOMINATOR  characteristics 

POLYNOMIAL  COEFr I  Cl  ENTS 
P<6)  -  -.68 100E-0 1 

P<  1)  «  15726E*00 

P( 2)  «  -.116  17E+02 

P( 3)  »  .  12057E*0  1 

P<4>  »  .  10006E+0  1 

POLYNOMIAL  FACTORS 

REAL  PART  I MAG  PARTS  OR  WN  PERIOD 

-.78548E-02  .76 164E-0  1  .092  .7649QE-81  ,62496E*62 

. 28675E+0  1 
-.40591E+8 1 

TF  CODES  <U,A0A,0,TH,N2 ,HDOT,STOP) 

TF  -  N2 

LX  -  9 

N2/DE  NUMERATOR  LX  *  8.000 

POLYNOMIAL  COEFFJ C1ENTS 

P<0)  «  -.31 792E-02 

PU)  -  . 2332OE+00 

P< 2)  «  .831 15£*62 

P< 3)  *  . 1 1336E+0 1 

P< 4)  *  .  9886 1E+-0  0 

POLYNOMIAL  FACTORS 


REAL  PART 

IMAG  PARTS 

DR 

WN 

PERIOD 

774S2E-02 

493S8E-02 

57 1?  l£«-88 

.9153 1E*B 1 

.062 

.9  1689E+0 1 

. 6S66 1E+08 

TF  *  NZ 
LX  »  ie 

NZ/DE  NUMERATOR  LX  «  10,000 

POLYNOMIAL  COEFFICIENTS 

P< 8)  *  -.31792E-02 

P<1)  -  . 23320E+00 

P<2)  »  .83 129E+02 

P(3>  *  .24923E+6 1 

P<4)  *  . 27930E+O 1 

polynomial  factgrs 

REAL  PART  I  MAG  PARTS  DR  t*N  PERIOD 

77450E-02 
•493S3E-02 

-.44476E+00  .S4372E+01  .882  .54554E+01  . 11556E+81 

TF  «  N2 

LX  «  20 

NZ/DE  NUMERATOR  LX  »  20.000 

POLYNOMIAL  COEFFICIENTS 

P<  8)  «  -.31 792E-8 2 

P(  I)  -  .2332OE+08 

P<2)  »  .83143E+02 

P<  3)  »  . 385 1OE*0 1 

P<4)  -  . 45974E+0 1 

POLYNOMIAL  FACTORS 

REAL  FART  IMAG  PARTS  DR  WN  PERIOD 

-.77448E-02 
. 49379E-02 

-.41 742E+00  ,42313E-r01  .098  .425 24£<-0  1  .  J4e48£*01 

TF  -  NZ 
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LX  -  36 

NZ/DE  NUMERATOR  LX  “  38.000 

polynomial  coefficients 

P<0>  “  - . 3 1792E-0 2 

P< 1)  =  .23320E-00 

P<2>  «  .82157£»02 

P( 3)  =  ,52697E»0 1 

P<  4)  =  . £48 1 SE»  0  1 


POLYNOMIAL  FACTOPS 

REAL  PART  IMAG  PARTS 

OR 

W4 

PERIOD 

-.77446E-02 

.49374E-02 

-.48549E+06  .35689E+01 

.113 

.36028E-81 

.  1754<SE*0  1 

TF  -  NZ 

JSJH ? 

NZ/OE  NUMERATOR  LX  =  40 

.060 

POLYNOMIAL  COEFFICIENTS 

P<0>  »  -.31 792E-02 

P(l)  «  . 23328E+0O 

P<2)  -  . 83 1 7 lE+62 

P<  3)  ■  .£5£84E*01 

P< 4)  ■  .8286IE+0  1 

* 

POLYNOMIAL  FACTORS 

REAL  PART  IMAG  PARTS 

DR 

LN 

PERIOD 

- . 77443E-02 
.49349E-82 

-.39681E+00  . 3 1382E+0 1 

.  125 

. 3 1832E+8  1 

. I9695E+0  1 

TF  *  NZ 


LX  -  58, 

NZ/DE  NUMERATOR  LX  -  50.368 


POLYNOMIAL  COEFFICIENTS 
P<0>  -  -.31792E-02 

P<  1)  =  .23320E+  CO 

P< 2)  “  .83165E+02 

P< 3)  -  .7927 1E*B I 

P< 4)  *  .  180  1 1E  +  02 

POLYNOMiAL  FACTORS 

REAL  PART  I MAG  PARTS 

-.7744  IE-82 
.49364E-02 

-.39454E+00  . 26552E+0 I 


DR 

PERIOD 

137 

. 28823E+0 1 

. 22087E+0 1 
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Appendix  D 

TTYLON  Runs  for  Flight  Envelope 

The  transfer  functions,  N 2/cc,  q/oc,  and  a/6c  were 
found  using  TTYLON.  Following  are  the  computer  listings 
showing  this  output  for  the  flight  conditions  given  in  Table 
I. 


I 
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ENTER  FLIGHT  CONDITION 
MACH  =  .3 


ALT  -  0 


WEIGHT  -  15000 
CG  -  451 
GAMMA  “  0 


FLIGHT  CASE  DATA 


MACH 

2 

.30060 

VFPS 

* 

335.47 

VSND 

* 

1118.22 

ALT 

x= 

0  . 

QDYN 

2 

133.74 

DEN 

2 

.8e237s* 

W6HT 

= 

15000  . 

1YY 

2 

47795. 

AREA 

m 

185.00 

CG 

s 

451.000 

REFCG 

= 

451.000 

CHORD 

= 

7 . 22 

GAMMA 

M 

0.000 

EPSTH 

2 

0.000 

ZT 

2 

6 . 0  O  O 

AOA 

» 

7.968 

ELEV 

m 

-1 1.796 

THRST 

1925.64 

NON-DIMENSI CNAL 
PERTURBATION 


0 

M 

A 

AD 

Q 

OE 


(PER 

(PER 

(PER 

(PER 

(PER 

(PER 


RAD) 

RAD) 

DEG) 

RAD) 

RAD) 

DEG) 


DIMENSIONAL 

PERTURBATION 


DERI  DAT  I  VE 
CD 

.07 7074S0 
.0  1481208 
.01720939 
0.00000000 
0.80000000 
.00 1 19177 

DERIVATIVE  PARAMETER 


U 

UD 

A 

AD 

G 

OE 

USE 


(PER 

(PER 

(PER 

(PER 

(PER 

(PER 


RAD) 

RAD) 

RAO) 

RAO) 

RAO) 

RAO) 


SHORT  PERIOD 


■ .  3  144SE-0  2 
• .  I536o£-04 
. 1 5734E+02 
.36474E-6I 
. 53466E+06 
*.  1  1069E+0 1 
APPROXIMATION 


MATRIX 

CL 

.59544855 
.33172654 
. 08583333 
. 4o976667 
6.75383333 
.00583867 
MATRIX 

2 

- . 949?®£-0 1 
. 10979E-O3 
-.2674l£»03 
26059E*00 
- . 38200E+0  1 
18236E+02 
(YES/NO) ?N 


CM 

0 .80000000 
.018554IO 
.0290  1  ©c.7 
-.96550000 
-7 . I8366©67 
.0  1403667 

M 

* . 25062E-02 
.  16046E-04 
. 60975E+0 1 
«. 3S0S8E-0 1 
-.28894E+0O 
.  30060E+0 1 


COMPUTE  ELEVATOR  TRANSFER  FUNCTIONS  (YES/NO) ?Y 

DENOMINATOR  CHARACTERISTICS 
POLYNOMIAL  COEFFICIENTS 
P<0)  -  -. 1200 1E+00 

P(l)  *  - . 14794E+80 


FERIOD 
, 4446  1 E  +  02 

- . 30539E+0  1 

TF  CODES  ( U , AOA , 0 , TH ,N2 , HDOT , STOP) 


P( 2)  =  -.5S964E*01 

P( 3)  »  . 1 1347E+0  l 
P( 4)  *  .  10080E+0  1 
POLYNOMIAL  FACTORS 

REAL  PART  I MAG  PARTS  DR  WN 

-. 1431 1E-0 I  .14132E+00  .101  . 14204E*00 

, 19478E+0 1 
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z 


LX  •  20 


NZ/OE  NUMERATOR  lx  *  20.800 

POLYNOMIAL  COEFFICIENTS 

p(0)  .  -.36576E-62 

p(  i)  m  -.34411E+00 

P<2)  m  .2105 6E*82 

p< 3)  m  . 184S1E+01 

p(4)  *  .24342E+01 


POLYNOMIAL  FACTORS 
REAL  PART 
-.73371E-02 
.23640E-0 1 
- . 39 156E+80 
TF  -  Q 


1MA6  PARTS 


.27 159E+0  1 


DR 


.  133 


Q/OE  numerator 

POLYNOMIAL  COEFFICIENTS 


P<0)  *  0. 

P( l)  m  .263 B6E~S 1 

P<2)  *  . 20955E+0 1 

p(3>  jr  .3O081E+01 


POLYNOMIAL  FACTORS 
REAL  PART 


IMAG  PARTS 


0. 

12828E-0  1 
-.(S8378E+00 
Tp  »  aOA 


DR 


.2942IE+01 


AOA/DE  NUMERATOR 
POLYNOMIAL  COEFFICIENTS 

p<0)  »  .sissze^i 

Pd)  »  .  17535E-0  1 

p<2)  «  . 29533E+0 1 

P<3)  -  -.54846E-01 


POLYNOMIAL  FACTORS 
REAL  PART 
-.30677E-02 
. 53653E+02 


IMAG  PARTS 
.  10336E+00 


DR 

.030 


WN 

. 18341E+00 


PERIOD 
. 2 1548E*0  1 

PERIOD 


PERIOD 
. 60790E+02 
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Y 


enter  flight  condition 

MACH  «  .6 
ALT  -  e 

WEIGHT  *  15600 


CG  -  451 
GAMMA  -  0 


FLIGHT  CASE  DATA 


MACH 

*S 

.60000 

VFPS 

670 .93 

VSND 

= 

1116.22 

ALT 

SS 

0  . 

QDYN 

= 

534.98 

DEN 

ts 

.0023769 

WGHT 

a 

15000  . 

IYY 

= 

47795. 

AREA 

= 

185.00 

CG 

a 

451.000 

REFCG 

= 

451.000 

CHORD 

s 

7 . 22 

GAWA 

a 

0.000 

EPSTH 

S 

0.000 

2T 

= 

0.000 

AOA 

« 

1.895 

ELEV 

.631 

THRST 

m 

2943.94 

NON-DIMENSIONAL 

PERTURBATION 


DERIVATIVE  MATRIX 


0 

M 

A 

AD 

Q 

DE 


(PER 

(PER 

(PER 

(PER 

(PER 

(PER 


RAD) 

RAD) 

DEG) 

RAD) 

RAD) 

DEG) 


DIMENSIONAL 


CD 

.02972913 
.0  1310585 
.00 189691 
0.00000000 
0.00000000 
.000 13091 

DERIVATIVE  PARAMETER 


PERTURBATION 
U  (PER  RAD) 
UD  (PER 
(PER 
(PER 
(PER 
(PER 


A 

AD 

0 

DE 

USE 


X 

- . 206S2E-0 5 

~.87344£-0«s 
.5063 1E+02 
.  17700E-81 
.  27359E+08 
. 13O67E+0  1 
SHORT  PERIOD  APPROXIMAT I  ON 


RAD) 

RAD) 

RAD) 

RAD) 

RAD) 


CL 

. 15057535 
- ^  00563595 
. 10439333 
.46906667 
7 . 24233333 
.00720467 
MATRIX 

2 

- . 3 1965E-0  1 
. 26396E-04 
I275SE+  04 
- . 53490E+00 
- . 82679E+0  I 
- . 87  o36E»  02 
(YES/NO) ?N 


CM 

-.00000000 
.  00589660 
. 03293©67 
-1  .08040000 
“7 . 6656666 / 
.0  1564oo7 

M 

-.  131 19E-02 
. 4264  IE-05 
.2S1S5E-02 
-.86SI6E-0 1 
-.61 666E*  00 
.  1 340  3E*  o 2 


COMPUTE  ELEVATOR  TRANSFER  FUNCTIONS  (YES/NO) ?Y 

DENOMINATOR  CHARACTERISTICS 
POLYNOMIAL  COEFFICIENTS 
P(0)  >=  -.  12255E+00 

P(l)  =  -  .  565  1 4E+  00 

P( 2)  *  - . 266  15E+02 

P( 3)  *  .2624  1E+0  I 

P( 4)  =  . 10000E+0  1 

POLYNOMIAL  FACTORS 

REAL  PART  I MAG  PARTS  DR  WN  PERIOD 

-.10S17E-01  . 669  1 OE-0  1  .160  .677 78E-0 1  .93906E  +  02 

.  4025 1E+0  1 
-.66276E+0  1 

TF  CODES  (U.AOA.Q.TH.NZ , HOOT, STOP) 
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NZ 


LX  -  20 

NZ/DE  NUMERATOR  LX  *  20.000 

POLYNOMIAL  COEFFICIENTS 

P<0>  -  -.U101E-01 

P< 1)  -  . 88329E+0 1 

P( 2)  *  . 45438E+0  3 

P(3)  »  . 19359E+02 

P<  4)  ■  . 11060E+02 

POLYNOMIAL  FACTORS 


REAL  PART 

IMAG  PARTS 

DR 

1*N 

-.21237E-0  1 

. 16499E-02 

86537E+00 

. 63482E+  0 1 

.  135 

.64069E+0 1 

TF  -  0 

O/DE  NUMERATOR 
POLYNOMIAL  COEFFICIENTS 
P<0)  -  0. 

P<  1)  -  .48683E+00 

P<  2)  =  .22075E+02 

P<  3)  »  . 134 13E+02 

POLYNOMIAL  FACTORS 

REAL  PART  IMAG  PARTS  OR  WN 

0. 

-.22357E-0  1 
16232E+6 1 
TF  -  AOA 

AOA/DE  NUMERATOR 
POLYNOMIAL  COEFFICIENTS 
P<0)  -  . 25*0  1E-0  1 

P<1)  -  . 26663E+00 

P( 2)  -  . 13  144E+02 

P< 3)  *  13O59E+00 

POLYNOMIAL  FACTORS 

REAL  PART  IMAG  PARTS  DR  WN 

- . 10 150E-0 1  . 42945E-0  1  .230  .44129E-01 

.  10068E+03 


PERIOD 

.  96976E+00 

PERIOD 


PERIOD 
. 1463 1E+03 
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ENTER  FLIGHT  CONDITION 
MACH  -  .8 


ALT  -  6 
WEIGHT  ■  15000 
CG  «  451 
GAMMA  «  0 


FLIGHT  CASE  DATA 


MACH 

e 

.88090 

VFPS 

894.53 

VSNO 

V 

1118.22 

ALT 

* 

0. 

QDYN 

* 

951.08 

DEN 

m 

.0023769 

WGHT 

= 

15008  . 

IYY 

= 

47795. 

AREA 

ss 

185 .90 

CG 

a 

451.000 

REFCG 

= 

451 .000 

CHORD 

= 

7.22 

GAMMA 

= 

0.000 

EPSTH 

S 

0.000 

ZT 

as 

0.000 

AOA 

c 

1.253 

ELEV 

=5 

1.926 

THRST 

as 

5741.71 

NON-DIMENSIONAL  DERIVATIVE  MATRIX 

PERTURBATION  CO  CL  CM 

0  (PER  RAD)  .03242504  .08453863  “.00000000 

M  (PER  RAD)  .01566414  -.04142369  .00445585 

A  (PER  DEG)  -.00232772  .11730000  .03555000 

AD  (PER  RAD)  0.00000000  .46860000  -1.15700080 

Q  (PER  RAD)  0.00000000  7.56800O00  -7.98700900 

DE  (PER  OEG)  -.00016033  .00808200  .01672000 

DIMENSIONAL  DERIVATIVE  PARAMETER  MATRIX 

PERTURBATION  X  Z  M 

U  (PER  RAD)  ~ . 349 16E-0 1  .41841E-02  -.12170E-82 

UD  (PER  RAO)  - . 38 I22E-06  . 17435E-04  .30325E-O5 

A  (PER  RAD)  .  13725E+03  - .  25462E-04  .54U4E-08 

AD  (PER  RAD)  .  15593E-8 1  - . 7 13 15E+00  12404E»00 

Q  (PER  RAD)  . 25 196E+00  - . 1 1523E+02  -.85667E»00 

DE  (PER  RAD)  . 72874E+0 1  - . 1 7464E+03  . 25462E+02 

USE  SHORT  PERIOD  APPROXIMATI ON  (YES/NO)  >4 

COMPUTE  ELEVATOR  TRANSFER  FUNCTIONS  (YES/NO) ?Y 

DENOMINATOR  CHARACTERISTICS 
POLYNOMIAL  COEFFICIENTS 
P(0)  «  -.10 16SE+00 

P(l)  »  16726E+0 1 

P(2>  *  -. 50828E+02 

P< 3)  -  .3S537E+0 1 

P(  4)  b  .  1080SE  +  0  1 
POLYNOMIAL  FACTORS 

REAL  PART  2 MAG  PARTS  DR  WN  PERIOD 

- .  16487E-0  1  . 4  1562E-0  1  .369  .44713E-0  1  .  15ll8E->-03 

•  54774E+0  l 
-.9303IE+0  1 

TF  CODES  <U,AQA,Q,TH,N2 , HOOT, STOP) 


c 
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N  Z 


LX  -  20 

NZ/OE  NUMERATOR  LX  *  20. 

POLYNOMIAL  COEFFICIENTS 
P<0)  -  -.47071E-01 

p<  1)  -  .58398E+02 

P<  2)  -  . 172 18E+04 

P( 3)  »  .54 195E+02 

P< 4)  -  . 2 1247E+02 

POLYNOMIAL  FACTORS 

REAL  PART  IMAG  PARTE 

- . 3474 1E-0 1 
.78775E-03 

-  - .  12572E+0 I  . 89048E*0 1 

TF  «  Q 

Q/OE  NUMERATOR 
POLYNOMIAL  COEFFICIENTS 
P<0)  *  0. 

PCI)  «  .21533E+0 1 

P< 2)  -  .A2754E+02 

P(3)  »  . 25487E+02 

POLYNOMIAL  FACTORS 

REAL  PART  IMAG  PARTS 

0. 

- . 34804E-0  1 
-.24275E+0  1 
TF  -  AOA 

AOA/DE  NUMERATOR 
POLYNOMIAL  COEFFICIENTS 
P<0)  -  . 3 1 18 IE-02 

P< X)  -  . 84457E+00 

P< 2)  -  . 24940E+02 

P( 3)  -  1951 1E+08 

POLYNOMIAL  FACTORS 

REAL  PART  IMAG  PARTS 

-.29638E-0 1 
-.42172E-02 
. 12786E+03 


DR 

WN 

PERIOD 

140 

.8993IE+0 1 

.  7B560E  +  00 

OR  WN  PERIOD 


OR  WN  PERIOD 
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Y 


ENTER  FLIGHT  CONDITION 
MACH  =  .4 

ALT  -  20609 

WEIGHT  =  15000 

CG  -  451 

GAMMA  *  0 

FLIGHT  CASE  DATA 

MACH  *  .40600  VFPS  =  415.40  VSND  *  103S.50 

ALT  -  20000.  QDYN  -  109.18  DEN  -  .0012655 

WGHT  =  15000.  IYY  =  47795.  AREA  «  185. 00 

CG  *=  451.008  REFCG  =  451.008  CHORD  -  7.22 

GAMMA  =  0.000  EPSTH  =  0.000  2T  -  0.000 

AOA  *=  9.582  ELEV  »  -14.768  THRST  «  2222.59 

NON-D I MENS  I ONAL  DERIVATIVE  MATRIX 

PERTURBATION  CD  CL  CM 

0  (PER  RAD)  .10855090  .72430230  -.00000800 

M  (PER  RAD)  .02146685  .28583536  .00047312 

A  (PER  DEG)  .02214529  .08547333  .028806*7 

AD  (PER  RAD)  0.00000000  .46953333  -1.00380003 

Q  (PER  RAD)  0.00000000  6.91666667  -7.34433333 

DE  (PER  OEG)  .00140297  .00541500  .01401667 

DIMENSIONAL  DERIVATIVE  PARAMETER  MATRIX 

PERTURBATION  X  Z  M 

U  (PER  RAD)  - . 12750E-02  -.75087E-01  -.20167E-02 

UD  (PER  RAD)  - . 1 1792E-04  .69B54E-04  . 10666E-84 

A  (PER  RAD)  . 12751E+02  -.222P7E+03  .48966E+01 

AD  (PER  RAD)  .2S612E-01  -.16949E+80  -.258S0E-01 

Q  (PER  RAO)  . 43350E+00  -.25679E+01  -.19475E+80 

DE  (PER  RAD)  -. 1 1965E+0 1  -. 13834E+02  .24504E+01 

USE  SHORT  PERIOD  APPROXIMATION  (YES/NO) ?N 

COMPUTE  ELEVATOR  TRANSFER  FUNCTI  ONS  (YES/NO)  ?Y 

DENOMINATOR  CHARACTERISTICS 
POLYNOMIAL  COEFFICIENTS 
P(0>  -  - . 63533E-0 1 

P(  1)  *  - .  1 1674E+00 

P( 2)  ■  - . 43925E+0 1 

P( 3)  «  . 76659E+00 

P( 4)  =  . 18000E+0 1 

POLYNOMIAL  FACTORS 

REAL  PART  IMAG  PARTS  DR  WN  PERIOD 

- . 12824E-0 1  . 1 1286E+00  .113  .  11359E+00  .55673E+02 

.  18793E  +  0 1 
-.26203E*0 1 

TF  CODES  < U, AOA, 0, TH ,NZ , HOOT, STOP) 
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NZ 


LX  ■  20 


NZ/DE  NUMERATOR  LX  ■  20.000 

POLYNOMIAL  COEFFICIENTS 

P<0>  *  -.17427E-02 

PC  1)  -  -.27060E+00 

P<2>  *  . 14424E+02 

P<3>  »  . 103 12E+0  1 

P< 4)  «  . 19536E+0 I 

POLYNOMIAL  FACTORS 


REAL  PART 

IMAG  PARTS 

DR 

WN 

PERIOD 

-.50776E-02 
.23761E-0 1 
-  .27326E+00 

. 27054E+0 1 

.100 

.27 192E+8 1 

. 23225E+0  1 

TF  ■  D 

Q/DE  NUMERATOR 
POLYNOMIAL  COEFFICIENTS 
P<0>  »  0. 

P<  1)  *  . 10469E-0 1 

P<2>  »  . 1 1736E+8  1 

PC  3)  »  . 245 13E+0 1 

POLYNOMIAL  FACTORS 

REAL  PART  IMAG  PARTS  DR  WN  PERIOD 

0. 

-.9093SE-02 
-.46966E+00 
TF  -  AOA 

AOA/DE  NUMERATOR 
POLYNOMIAL  COEFFICIENTS 
PC  9)  -  . 16333E-0 1 

P(l)  ■  . 90 149E-02 • 

P< 2)  *  .24276E+0  1 

P< 3)  -  - . 33760E-0 1 

POLYNOMIAL  FACTORS 

REAL  PART  IMAG  PARTS  DR  WN  PERIOD 

-.  19034E-02  . 8 1999E-0  1  .023  .82021E-01  .7*6 25E+02 

,7191 1E+02 
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ENTER  FLIGHT  CONDITION 
MACH  »  .6 

ALT  ■  20808 

WEIGHT  •  15000 

CG  »  451 

GAMMA  =  0 


FLIGHT  CASE  DATA 


MACH 

=* 

.60000 

VFPS 

* 

623.  10 

VSND 

r= 

1038.50 

ALT 

m 

28000  . 

QDYN 

*= 

245 . 66 

DEN 

M 

.86 12o55 

WGHT 

= 

15080  . 

IYY 

47795. 

AREA 

= 

135.00 

CG 

451.000 

REFCG 

* 

45 1.000 

CHORD 

m 

•7  ->  r> 

GAMMA 

= 

0.000 

EPSTH 

= 

0.000 

2T 

= 

0.000 

AOA 

4.035 

ELEV 

* 

-3.333 

THRST 

= 

1814.74 

NON-DIMENSIONAL  DERIVATIVE  MATRIX 

PERTURBATION  CD  CL  CM 

8  (PER  RAD)  .03983219  .32724656  -.08008008 

M  (PER  RAD)  .00983481  .07557654  -.002519^2 

A  (PER  DEG)  .80900309  .09358667  .03043223 

AD  (PER  RAD)  0.00000000  .46906667  -1.06040068 

Q  (FER  RAD)  0.00000000  7.24233333  -7 . 665<s66©7 

DE  (PER  DEG)  .00054575  .00567300  .01475333 

DIMENSIONAL  DERIVATIVE  PARAMETER  MATRIX 

PERTURBATION  X  2  M 

U  (PER  RAD)  - . 77483E-02  -.50411E-01  -.13636E-02 

UD  (PER  PAD)  - . 2 1 056E-05  .29846E-04  .4S534E-05 

A  (PER  RAD)  . 18380E+02  -.52378E+83  .11911E+02 

AD  (PER  RAD)  . 1855 1E-0 1  -.26295E+08  -.42760E-01 

Q  (PER  RAD)  . 26784E+00  - . 4080  1E+0  1  -.304®0E*60 

DE  (PER  RAD)  -.81075E+00  -.31821E+02  .58032E+0I 

USE  SHORT  PERIOD  APPROXIMATION  (YES/NO) ?N 

COMPUTE  ELEVATOR  TRANSFER  FUNCTIONS  (YES/NO) ?Y 

DENOMINATOR  CHARACTERISTICS 
POLYNOMIAL  COEFFICIENTS 
P(0)  *  - . 68  108E-0 1 

P( 1)  =  15726E+00 

P( 2)  =  -.116  17E+02 

P( 3)  *  . 12057E+0 1 

P( 4)  *  . 10008E+O 1 

POLYNOMIAL  FACTORS 

REAL  PART  IMAG  PARTS  DR  WN  PERIOD 

- . 70548E-02  .76164E-01  .092  .76490E-01  .S2496E+02 

. 23675E+8 1 
-. 4059 1E+0 1 

TF  CODES  (U.AOA.Q.TH.NZ, HOOT, STOP) 
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Y 


ENTER  FLIGHT  CONDITION 
MACH  -  1 

ALT  -  20000 

WEIGHT  -  15000 

CG  -  451 

GAMMA  -  0 


FLIGHT  CASE  DATA 


MACH 

* 

1.00000 

VFPS 

■ 

1038.50 

VSND 

S 

1838.50 

ALT 

20000  . 

QDYN 

* 

682.38 

DEN 

S 

.08 12655 

WGHT 

15000  . 

IYY 

=* 

47795. 

AREA 

= 

185.00 

CG 

* 

451.000 

REFCG 

** 

451.000 

CHORD 

= 

7.22 

GAMMA 

e 

0.000 

EPSTH 

0.000 

2T 

= 

0.000 

AOA 

1.719 

ELEV 

= 

.748 

THRST 

* 

4448.03 

NON-DIMENSIONAL  DERIVATIVE  MATRIX 

PERTURBATION  CD  CL  CM 

0  (PER  RAD)  .83521861  . 11776344  0.00000000 

M  (PER  RAD)  .01495357  -.06609 252  .05778858 

A  (PER  DEG)  -.00092374  .09711750  .02259650 

AD  (PER  RAD)  0.00000000  .45172500  -1.09200000 

Q  (PER  RAD)  0.00000000  6.83175000  -7.62075000 

DE  (PER  DEG)  -.00009914  .01042325  .01307525 

DIMENSIONAL  DERIVATIVE  PARAMETER  MATRIX 

PERTURBATION  X  Z  M 

U  (PER  RAD)  - . 23569E-0 1  -.18673E-02  .34752E-03 

UD  (PER  RAD)  -.36844E-06  .  12277E-04  .2091  1E-O5 

A  (PER  RAD)  . 90949E+82  -.15142E+04  .24701E+02 

AD  (PER  RAD)  . 12744E-0  1  -.42462E+00  -.72325E-0  1 

Q  (PER  RAO)  . 19290E+O8  -.64277E+01  -.50519E+0O 

DE  (PER  RAD)  .63884E+0  1  -.16159E+03  .  14287E+02 

USE  SHORT  PERIOD  APPROXIMATION  (YES/NO) ?N 

COMPUTE  ELEVATOR  TRANSFER  FACTIONS  (YES/NO)  ?Y 

DENOMINATOR  CHARACTERISTICS 
POLYNOMIAL  COEFFICIENTS 
P(0)  -  . 16050E-0 1 

P(  1)  -  - . 543 16E+00 

P( 2)  -  - . 23742E+02 

P( 3)  «  .20588E+0  1 

P( 4)  ■  . 10800E+0 1 

POLYNOMIAL  factors 

REAL  PART  IMAG  PARTS  DR  1«M  PERIOD 

- . 39750E-0 1 
. 16974E-0 1 
. 39644E+0 1 
- . 60004E+0 1 

TF  CODES  <U,AOA,Q,TH,N2 , HOOT, STOP) 
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NZ 


LX  *  20 


20.000 


( 


nz/de  numerator  IX  ' 

POLYNOMIAL  COEFFI Cl  ENTS 
P(g)  *  - . 1 1 996E-0 i 

p<  j>  m  .  l2433E<-02 

p<2)  »  .54352E+03 

p(3)  m  .leS3l£+02 

p<4)  »  .13903E+02 

P0LYN^L4£T0RS'»«  parts  OR 
-.23609E-01 

TF  -  Q 

Q/DE  NUMERATOR 
POLYNOMIAL  COEFFICIENTS 

P<0)  *  0  • 

pcu  .  ,  39983E+00 

p<2>  =  .17328E+02 

p<3>  .  .  1429SE+-02 

POLYNOMIAL  FACTORS 

REAL  PART  I MAG  PAK l 3 

0. 

-.23535E-01 
1 1S84E+0 1 
TF  “  AOA 

AOA/DE  NUMERATOR 
POLYNOMIAL  COErFI Cl  ENT 3 
P(0)  =  - . 158 19E-02 

p< d  «  .32033E+00 

p( 2)  «  .  14  1 10E+02 

p<  3)  «  - . I55o 1E+00 

POLYNOMIAL  FACTORS  paRTS  OR 

REAL  PART  I MAG  PAR i = 

- .26867E-0 1 
.41717E-02 
.  90699E+02 


WN 


.62779E+0 1 


WN 


WN 


c 
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L 


PERIOD 
. 10054E+0  1 

I 

PERIOD 


PERIOD 


ENTER  FLIGHT  CONDITION 
(V»CH  =  1.2 

ALT  =  20000 

WEIGHT  -  15000 

CG  ■  451 

GAMMA  *  0 


FLIGHT  CASE  DATA 


MACH 

= 

1.20000 

VFPS 

* 

1246.20 

VSND 

S 

1038.50 

ALT 

S 

20060. 

QDYN 

= 

982.63 

DEN 

= 

.00  12655 

WGHT 

= 

15000  . 

1YY 

= 

47795. 

AREA 

= 

185.00 

CG 

451 .000 

REFCG 

451.000 

CHORD 

m 

7.22 

GAMMA 

*: 

0.000 

EPSTH 

= 

0.000 

ZT 

= 

0.000 

AOA 

= 

1.625 

ELEV 

= 

-.220 

THRST 

s 

7126.20 

NON-DIMENSIONAL  DERIVATIVE  MATRIX 

PERTUP.SATI  ON  CD  CL  CM 

0  ( PER  RAD)  .03928507  .08140251  0.00000000 

M  (PER  RAD)  .01623864  -.08568897  .07392841 

A  (PER  DEG)  -.00358194  .09253560  .01313300 

AD  (PER  RAD)  0.00000000  .43485000  -1.02706000 

Q  (PER  RAD)  0.00000000  6.09550000  -7.25456006 

DE  (PER  DEG)  -.00057736  .01491550  .01006056 

DIMENSIONAL  DERIVATIVE  PARAMETER  MATRIX 

PERTURBATION  X  2  M 

U  (PER  RAD)  - . 3394SE-0  1  .2768OE-01  .  14839E-62 

UD  (PER  RAD)  -.31696E-66  . 11173E-04  . 1S591E-65 

A  (PER  RAD)  .1695SE+03  - .20776E+04  .20716E+02 

AD  (PER  RAD)  . 13918E-01  -.49653E+60  -.81632E-61 

Q  (PER  RAO)  . 19525E+00  -.68S23E+01  -.57709E*00 

DE  (PER  RAD)  .22343E+02  -.33272E+03  . I6772E+02 

USE  SHORT  PERIOD  APPROX I MAT I  ON  (YES/NO) ?N 

COMPUTE  ELEVATOR  TRANSFER  FUNCTIONS  (YES/NO) ?Y 

DENOMINATOR  CHARACTERISTICS 
POLYNOMIAL  COEFFICIENTS 
P(0)  =  . 9506 1E-0 1 

P( 1)  =  - . 7527 1E+00 

P( 2)  =  1 9504E+02 

P( 3)  =  . 23594E+0  1 

P(4)  =  . 10000E*  0 1 

POLYNOMIAL  factors 

REAL  PART  IMAG  PARTS  DR  WN  PERIOD 

-.91 1 17E-01 
. 53264E-0 1 
. 34 147E+0  1 
- . 57362E*0  1 

TF  CODES  (U.AOA,Q,TH,NZ , HOOT, STOP) 


NZ 


c 


LX  -  20 

NZ/OE  NUMERATOR  *-*•  1 

POLYNOMIAL  COEFFICIENTS 
P(0> 

P<  D 
P<  2) 

P<  3) 

P<4) 


20.000 


parts 


1975SE-0  1 
.27297E+02 
.8703SE+03 
. 25304E+02 
.20776E+02 
POLYNOMIAL  FACTORS 

REAL  PART  I MAO 
-.32099E-0 1 
.70784E-03  „  _ 

_ -59323£'t-00  44^1E+0  1 

TF  *  Q 

Q/DE  numerator 

POLYNOMIAL  COEFFICIENT a 

P<0>  =0- 
p(D  s  .69674E+00 

p(2)  **  .23034E+02 

pc3)  =  .  16795E+02 

POLYNOMIAL  FACTORS  . 

REAL  PART  IMAG  PARTS 

0. 

-.30944E-0  1 
13406E+0  1 
TF  *  AOA 

aoa/pe  numerator 

POLYNOMIAL  COEFFICIENTS 
P(0)  —  -.25165E-0I 

p<  1)  «s  .55474E+00 
P<2)  *  • 1^51 1E+02 

p<3)  »  -.2<4<S99E+00 

POLYNOMIAL  FACTORS  p<sRTS 

REAL  PART  IMAG  PARTS 

-.59260E-01 
.2570SE-0 1 
. 6 1875E+02 


OR 


.092 


WN 


.64694E+0 1 


PERIOD 


?"533E+00 


I  i 


DR 


l*IN 


PERIOD 


DR 


WN 


PERIOD 


75 


Y 


enter  TLIGHT  CONDITION 
MACH  »  .5 


ALT  -  40080 
WEIGHT  »  J5000 
CG  -  451 
GAMMA  =  0 


PLIGHT  CASE  DATA 


MACH  * 

.50000 

VPPS  « 

484.81  VSND 

53 

969.62 

ALT  * 

40000. 

DDYN  - 

6S.70  DEN 

S3  . 

0  0  0  5  c  -  4 .3 

WGHT  « 

15000. 

IYY  - 

47795.  AREA 

9 

185.00 

CG 

45! 

.000 

REFCG  » 

451.000  CHORD 

=3 

7.22 

GAMMA-  * 

0 

.000 

EPSTH  « 

8.000  2T 

S3 

0 .000 

AOA  * 

14 

.537 

ELEV  = 

-24.827  THRST 

=s 

3246. 79 

NON-OIMENSIONAL 

DERIVATIVE  MATRIX 

PERTURBATION 

CD 

CL 

CM 

0  C  PER 

RAO) 

. 24729548 

1 . 1 1615331 

0.00006080 

M  C  PER 

RAO) 

. 05S70757 

. 37327587 

-.019 16480 

A  C  PER 

DEG) 

.03742975 

.83667000 

.02378008 

AO  CPER 

RAO) 

0.00000000 

.46938000 

1.04 

21088 8 

Q  CPER 

RAD) 

0.00008880 

7 . 07950800 

7.50 

50  0  6  0  0 

DE  CPER 

DEG) 

.00218394 

.00505780 

.01 

393000 

DIMENSIONAL  1 

DERIVATIVE  PARAMETER  MATRIX 

PERTURBATION 

X 

2 

M 

U  CPER 

RAD) 

.  10245E-02 

- . 59739E-8 1 

.  167 

57E-02 

UD  CPER 

RAD) 

1237PE-04 

.47 740E-04 

.77  1 

23E-05 

A  CPER 

RAD) 

.86572E+0 1 

-.15 136E+03 

.28*13E*0 1 

AO  CPER 

RAD) 

.22 404E-0! 

- .  8640 1E-0 1 

.  13959E-0 1 

Q  CPER 

RAO) 

.  36069E+08 

-. 13910E+3  1 

.  107 

29E-0O 

OE  CPER 

RAD) 

'  13193E+0 I 

858 17E+0  1 

.  153 

23E*0  1 

USE  SHORT  PERIOD  APPROXIMATI ON  CYES/NO) ’N 

COMPUTE  ELEVATOR  TRANSFER  FUNCTIONS  C YES/NO) ?Y 

DENOMINATOR  CHARACTERISTICS 
POLYNOMIAL  COEFFICIENTS 
P<0>  »  - . 2S927E-6 1 

P< 1)  =  - , 95033E-0 I 

PC  2)  «  -.3U9?E*01 

PC 3)  *  .44466E+08 

PC  4)  =  .  18000E+0  1 

POLYNOMIAL  FACTORS 

REAL  PART  IMAG  PARTS  DR  WN  PERIOD 

- .  1572SE-0 1  .94651E-01  .  164  .95949E-01  .66333E+02 

.  I578OE  +  0 1 
199 12E+0  I 

TF  CODES  CU,AOA,Q,TH,N2 , HOOT .STOP) 
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NZ 


LX  -  20 

NZ/DE  NUMERATOR  LX  *  20.060 

POLYNOMIAL  COEFFICIENTS 

PC0>  ■  -.68040E-03 

PCI)  *  -.  17381E  +  00 

P( 2)  «  . 605 12E+0 1 

PC  3)  =  . 3754  1E+00 

P<  4)  =  . 12 168E+0  1 

POLYNOMIAL  FACTORS 


REAL  PART 

IMAG  PARTS 

DR 

WN 

PERIOD 

- . 34906E-02 
.321 49E-0  1 
1685SE+00 

.2225S£f 0  l 

.0  76 

.22322E+0 1 

. 28229E+0 1 

TF  «=  Q 

Q/DE  NUMERATOR 
POLYNOMIAL  COEFFICIENTS 
P(0>  =  0. 

P< 1)  =  . 27 10SE-02 

P<2)  =  .44275E+00 

PC  3)  “  . 15325E*0  1 

POLYNOMIAL  FACTORS 

REAL  PART  IMAG  PARTS  DR  WN  PERIOD 

0  . 

-.62581E-02 

- . 28265E+00 
TF  *  AOA 


AOA/DE  NUMERATOR 
POLYNOMIAL  COEFFICIENTS 
PC  0 )  =  . 70078E-02 

PC  1)  =  . 16300E-02 

PC  2)  *  .  1525oE+0 1 

PCS)  ■  -.181 13E-0 1 

POLYNOMIAL  FACTORS 

REAL  PART  IMAG  PARTS  DR  WN 

-.56145E-03  .67771E-01  .008  .67774E-01 

. 8423 1E+02 


PERIOD 

.F2712E+02 
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ENTER  FLIGHT  CONDITION 
MACH  -  .6 

ALT  *  40600 

WEIGHT  «  15000 

CG  -  451 

GAMMA  »  0 


FLIGHT  CASE  DATA 


MACH 

m 

.60000 

VFPS 

= 

581.77 

VSND 

= 

ALT 

*S 

40000  . 

GDYN 

= 

98.92 

DEN 

= 

.0005840 

WGHT 

15006  . 

1YY 

= 

47795. 

AREA 

S 

185.00 

CG 

zs 

451.000 

REFCG 

2 

45 1 .000 

CHORD 

= 

T . 

GAMMA 

= 

0.000 

EPSTH 

= 

0.000 

2T 

s 

0.000 

AOA 

* 

10.0  19 

ELEV 

= 

-15.382 

THRST 

= 

2480 . 39 

NON-DIMENSIONAL  DERIVATIVE  MATRIX 

PERTUR8AT I  ON  CD  CL  CM 

0  (PER  RAD)  .13346745  .79605*87  0.00000000 

M  (PER  RAD)  .02637223  .23666727  -.015535*5 

A  (PER  DEG)  .02722351  .03977006  .02939**7 

AD  (PER  RAD)  0.00000000  .46906667  -1.0S0400OO 

Q  (PER  RAO)  0.00000000  7.24233333  -7.665666*7 

DE  (PER  DEG)  .00154460  .00509333  .01426333 

DIMENSIONAL  DERIVATIVE  PARAMETER  MATRIX 

PERTURBATION  X  2  M 

U  (PER  RAD)  - . 40334E-03  -.54172E-01  -.14361E-02 

UO  (PER  RAO)  - , 59446E-05  .33647E-04  .55426E-05 

A  (PER  RAD)  . 64750E+8 1  -.2146SE+03  .45ll0E*0t 

AD  (PER  RAD)  .19276E-0I  -.1O911E+O0  -.17973E-01 

Q  (PER  RAD)  , 30692E+00  -.17372E+01  -.13 156E+C0 

DE  (PER  RAD)  -.14279E+01  - . 1  1685E+02  .22593E+01 

USE  SHORT  PERIOD  APPROXIMATION  (YE3/N0) ?N 

COMPUTE  ELEVATOR  TRANSFER  FUNCTIONS  (YES/NO) ?Y 

DENOMINATOR  CHARACTERISTICS 
POLYNOMIAL  COEFFICIENTS 
P ( 0 )  =  - . 30632E-0 1 

P< 1)  =  - . 9 1079E-0 1 

P( 2)  *  - . 45909E  +  0 1 

P( 3)  =  . 525O3E+00 

P(4)  *  .  10OOOE  +  0  1 

POLYNOMIAL  FACTORS 

REAL  PART  IMAG  PARTS  DR  WN  PERIOD 

-. 10247E-0 1  . 8B888E-0  1  .  126  .81534E-01  .7?678E +02 

.  19O91E+0 1 
-.241 36E+0 1 

TF  COOES  (U.AOA^.TH.NZ.HDOT.STOP) 
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N2 


LX  -  20 

N2/DE  NUMERATOR  LX  »  20.000 

POLYNOMIAL  COEFFICIENTS 

PC0)  -  -.56<S62£-03 

P<  1)  -  -.21030E+O0 

PC2)  *  .  1297 IE+02 

PC 3)  -  ,6S2B6C*ee 

PC 4)  »  .  17740E+0  1 

polynomial  factors 


REAL  PART 

IMAG  PARTS 

DR 

WN 

PERIOD 

23S29E-02 
18550E-0 1 
192  1 1E+00 

.26983E+0 1 

.071 

.27052E+0 1 

.  22286E+0 1 

TF  ■  Q 

Q/DE  NUMERATOR 
POLYNOMIAL  COEFFICIENTS 
PC 0)  «  0. 

PCI)  *  . 32543E-02 

PC  2)  «  . 75584E+60 

PC  3)  =  . 2259SE+0 1 

POLYNOMIAL  FACTORS 

REAL  PART  IMAG  PARTS  DR  WN  PERIOD 

0. 

- . 43659E-02 
-.330 13E+00 
TF  »  AOA 

AOA/DE  NUMERATOR 
POLYNOMIAL  COEFFICIENTS 
PCO)  »  . 7682  IE-02 

PCI)  -  . 54500E-02 

PC  2)  =  . 22494E+0  1 

PC 3)  *  -.20742E-01 

POLYNOMIAL  FACTORS 

REAL  PART  IMAG  PARTS  DR  WN  PERIOD 

-.12272E-02  .38427E-01  .021  .5S439E-01  .  I0754E+03 

. 10845E+03 
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Y 


ENTER  FLIGHT  CONDITION 
MACH  *  1 


ALT  -  40000 
WEIGHT  -  15000 
CG  -  451 
GAMMA  -  0 


FLIGHT  CASE  DATA 


MACH 

= 

1 .00000 

VFPS 

= 

969.62 

VSND 

= 

969.62 

ALT 

40000  . 

QDYN 

=5 

274.79 

DEN 

■ 

.0005846 

WGHT 

15003. 

IYY 

*= 

47795. 

AREA 

= 

185.00 

CG 

sr 

451.000 

REFCG 

SS 

451.000 

CHORD 

* 

7.22 

GAMMA 

* 

0.000 

EPSTH 

3 

0.000 

ZT 

0.000 

AOA 

* 

3.873 

ELEV 

X 

-1.969 

THRST 

■ 

2147.24 

NON-DIMENSIONAL  DERI  VAT I WE  MATRIX 


PERTURBAT I  ON 

CD 

CL 

CM 

0 

(PER 

RAD) 

.04214259 

.29221647 

-.00000000 

M 

(PER 

RAD) 

.0 1697672 

-. 15295772 

-.04446168 

A 

(PER 

DEG) 

.00S19174 

.09215500 

.02202400 

AD 

(PER 

RAD) 

0.00000000 

.45172500 

-1.09200000 

Q 

(PER 

RAD) 

0.00000000 

6.83175000 

-7.62075000 

DE 

(PER 

DEG) 

.00072353 

. 00S 13950 

.01293*50 

D I  MENS  I  ONAL 

DERIVATIVE  PARAMETER  MATRIX 

PERTURBATION 

X 

Z 

M 

U 

(PER 

RAD) 

-.9453SE-02 

-. 86439E-02 

-. 10264E-02 

UD 

(PER 

RAD) 

- . 86282E-06 

. 12745E-04 

. 2 1749E-05 

A 

(PER 

RAD) 

. 19307E+02 

- . 58088E+03 

. 96232E*  O  1 

AD 

(PER 

RAD) 

. 12330E-0  1 

-. 18213E+O0 

- . 3 10  79E-0  1 

0 

(PER 

RAD) 

. 16733E+00 

- . 2767 1E+0  1 

-.21 783E»00 

DE 

(PER 

RAO) 

-. 10753E+0  1 

- . 5 1040E+02 

. 57 140E+0 1 

USE  SHORT  PERIOD  APPROXIMATION  (YES/NO) ?N 

COMPUTE  ELEVATOR  TRANSFER  FUNCTIONS  (YES/NO) ?Y 

DENOMINATOR  CHARACTERISTICS 
POLYNOMIAL  COEFFICIENTS 
P<0>  *  - . 22379E-0 1 

P(l>  -  126S0E+00 

P<  2)  =  - . 952 18E+0  1 

P<  3)  «  . 85880E+00 

P(  4)  *  .  10O00E+0  1 

POLYNOMIAL  FACTORS 

REAL  PART  IMAG  PARTS  DR  WN  PERIOD 

-.67527E-02  .47974E-01  .139  .48445E-0  1  . 13098E+03 

. 2694 1E+0 1 
-.35394E+0 l 

TF  CODES  <U,A0A,0,TH,NZ , HOOT, STOP) 
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LX  -  20 


N2 


N2/DE  NUMERATOR  LX  *  20.000 

POLYNOMIAL  COEFFICIENTS 

P<0>  -  -.20544E-02 

P<  1)  »  .6  1S<41E»00 

P<  2)  -  . 87699E+02 

P<  3)  -  . 27847E+0 1 

P< 4)  -  .51390E+0 1 

POLYNOMIAL  FACTORS 


REAL  PART 

IMAG  PARTS 

DR 

- . 94906E-02 

.24638E-02 

-.26742E+00 

.4 1219E+0  1 

.065 

. 4 1306E+0 I 

TF  -  Q 

Q/DE  NUMERATOR 
POLYNOMIAL  COEFFICIENTS 
P<  0)  -  0. 

PCI)  =»  . 304 16E-0 1 

PC  2)  =  . 2977SE+0  1 

PC 3)  -  . 57 156E+0  1 

POLYNOMIAL  FACTORS 

REAL  PART  IMAG  PARTS  DR  UM 

0  . 

10423E-0 1 
-.51057E+00 
TF  -  AOA 

AOA/DE  NUMERATOR 
POLYNOMIAL  COEFFICIENTS 
PC0>  -  . 32527E-02 

PCI)  «  . 48877E-0  1 

PC  2)  -  .56845E+01 

PC  3)  «  52750E-0  I 

POLYNOMIAL  FACTORS 

REAL  PART  IMAG  PARTS  DR  WN 

- . 430 I4E-02  . 23530E-0  1  .180  .23920E-01 

. 10777E+03 


PERIOD 

. 15242E+0  1 

PERIOD 


PERIOD 
.  26703E+03 
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ENTER  FLIGHT  CONDITION 
MACH  »  1.4 

ALT  -  40000 

WEIGHT  ■  15000 

CG  *  451 

GAMMA  *=  0 


FLIGHT  CASE  DATA 

MACH  * 

1.40000 

VFPS  * 

1357, 46  VSND 

=  9 69.  £-2 

ALT 

40000  . 

QDYN  - 

538.59  DEN 

=  .8085843 

WGHT  * 

15003 . 

IYY  = 

47795.  AREA 

=  1S5.O0 

CG 

451.000 

REFCG  = 

451.000  CHORD 

*  7.22 

GAMMA  = 

0.000 

EPSTH  = 

0.000  2T 

=  0.000 

AOA  «= 

2.449 

ELEV  = 

.739  THRST 

=  40S1.30 

NON-DIMENSIONAL 

DERIVATIVE  MATRIX 

PERTURBATION 

CD 

CL 

CM 

0  (PER 

RAD) 

.04092404 

. 14879555 

- . 00800000 

M  (PER 

RAO) 

.01432249 

- .08553634 

-.00949051 

A  (PER 

DEG) 

.00022160 

.08452500 

.00481200 

AD  (PER 

RAD) 

0.00000000 

.41797500 

- .  9620000  0 

Q  (PER 

RAD) 

0.00000000 

5.35925000 

6 . 88825000 

DE  (PER 

DIG) 

.00003.593 

.0  1408.350 

. 0090  9950 

DIMENSIONAL  DERIVATIVE  PARAMETER  MATRIX 

PERTURBATION 

X 

2 

M 

U  (PER 

RAD) 

-. 17137E-0  1 

. 4 1495E-62 

. 27780E-O  2 

UD  (PER 

RAD) 

-.31944E-0O 

. 74700E-05 

.121 19E-85 

A  (PER 

RAD) 

.  72<5<5OE+02 

-. 104  13E  +  04 

. 4 1 337E+  O 1 

AD  (PER 

RAD) 

.  10131E-0 l 

-.2369  1E+00 

.  38436E-0 1 

Q  (PER 

RAD) 

. 130 14E+00 

-.30431E+O  1 

•27572E+00 

DE  (PER 

RAD) 

.69 177E+0  1 

- . 1 7235E+  O  3 

. 73473E+0  1 

USE  SHORT  PERIOD 

APPROXIMATION 

(YES/NO) 15 

COMPUTE  ELEVATOR  TRANSFER  FUNCTIONS  (YES/NO) ?Y 

DENOMINATOR  CHARACTERISTICS 
POLYNOMIAL  COEFFICIENTS 
P(0)  =  -.63973E-02 

P(l)  *  - . 63237E-0  1 

P< 2)  *  39098E+0 1 

P<  3)  =  . 10989E+O 1 

P(4)  -  .  10000E+0  1 

POLYNOMIAL  FACTORS 

REAL  PART  IMAG  PARTS  DR  WN  PERIOD 

-.82773E-02  .39491E-01  .205  .40349E-61  . 15910E*03 

. 15136E+0 1 
- . 25960E+0  1 

TF  CODES  <U,AOA,Q,TH,N2 , HOOT, STOP) 
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NZ 


LX  ■  20 

NZ/DE  NUMERATOR  LX  *  20.000 

POLYNOMIAL  COEFFICIENTS 

P<0>  -  - . 400 18E-02 

PCI)  -  . 36980E+0 l 

P<2)  •  .23167E+03 

P<3)  «  . 60745E+0 l 

PC  4>  «  .  10237E+02 

POLYNOMIAL  FACTORS 

REAL  PART  I MAG  PARTS  C 

-.  16973E-01 
. 10 173E-02 

- . 28872E+00  . 47495E+0  1  ,6< 

TF  «  Q 

Q/DE  NUMERATOR 
POLYNOMIAL  COEFFICIENTS 
PC0>  ■  0. 

P<1>  «  .93665E-01 

P<  2)  -  . 563 15E+0 1 

P( 3)  «  .78 522E+01 

POLYNOMIAL  FACTORS 


.47582E+0 1 


PERIOD 


1  3229E+0 1 


I MAG  PARTS 


REAL  PART 

0. 

-.  17837E-0 1 
- . 700 18E+00 
TF  »  AOA 


AOA/DE  NUMERATOR 
POLYNOMIAL  COEFFICIENTS 
PC0>  *  . 22863E-03 

P<  1)  -  .  12433E+00 

P( 2)  ■  . 779 1 1E+0 1 

PC  3)  *  - . 12706E+00 

POLYNOMIAL  FACTORS 

REAL  PART  I MAG  PARTS 

13834E-0 1 
-.21207E-02 
.61334E+02 


PERIOD 


PERIOD 
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Appendix  E 


Root  Loci 

The  root  locus  design  method  was  the  frequency  domain 
technique  used.  The  computer  program  TOTAL,  which  is  main¬ 
tained  by  AFIT/ENE  at  Wright-Patterson  AEB,  was  used  to  get 
listings  and  plots  of  root  loci.  First  the  open  loop  NZ/6C 
root  loci  were  listed  for  each  flight  condition.  Following 
are  these  listings  with  only  the  short  period  branch  listed. 
One  can  see  that  the  damping  ratio  (zeta)  and  gain  at  each 
point  in  the  locus  are  listed.  As  is  explained  in  the  sec¬ 
tion,  "Selection  of  Gains,"  Knz  =  -1.2  was  chosen.  Follow¬ 
ing  the  open  loop  Nz  root  loci  are  the  listings  of  the  roots 
which  result  from  the  K$z  =  -1.2.  Option  42  of  TOTAL  was 
used  to  get  these  results. 

Once  the  Nz  loop  was  closed,  the  roots  of  q/6c  were 
found  with  the  desired  damping  ratio,  0.7.  Option  43  of 
TOTAL  was  used  to  get  these  results.  Finally,  Option  42  was 
used  to  get  the  roots  of  q/<5c  which  resulted  from  the  gain 
schedules  for  Kq  in  Methods  I,  II,  and  III. 
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GAIN™- 1,41 


M  =  0.3 

h  =  0  ft 


OPEN-LOOP  < OLTF)  ROOT  LOCUS  USING  OPTION  4  1 
CODE 

8 -LOCUS  PT.  2-2 ERO  4- 1 MAG  INARY  AXIS 

1-POLE  3-BREAK  PT .  3-SENSITIVITY  PT. 

REGION  OF  CALCULATION-REAL:  CC=  -15.0  TO  AA*  18.8 

IMAJ:  DD=>  -.188E-0  1  TO  BB*  15.0 

BRANCH  STARTING  AT  <-.814311)  ♦  J< . 14132) 

TYPE  L  TO  LIST,  S  Tu  SKIP,  OR  S  TO  ABORT  >S 

BRANCH  STARTING  AT  <-.014311)  ♦  J<-. 14132) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  *  TO  ABORT  >S 

BRANCH  STARTING  AT  <1,9478)  ♦  J<0.) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  *  TO  ABORT  >S 

BRANCH  STARTING  AT  <-3.0539)  ♦  J<0.) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  *  TO  ABORT  >L 


BRANCH  NUMEER  4 


CALCULATION  STEP  SIZE  »  .1250 

PRINTING  STEP  SIZE  -  .4250 


LOCUS  REAL 

LOCUS  I MAG 

DIST  TO  ORIGIN 

GAIN 

z  eta 

CD 

-3.8539088 

0  . 

3.8539008 

0. 

1.00000 

1 

-3.4789088 

0. 

3.4789800 

. 753 159E-0 1 

1 .00000 

0 

-4.3039800 

0  . 

4.3039000 

.  135557 

1.00000 

0 

-4.9289000 

0. 

4.9289000 

.  IS2892 

1.00080 

0 

-5.5539000 

0. 

5.5539000 

.219992 

1 .00000 

0 

-4. 1789000 

0  . 

4.1789000 

.249224 

1.00000 

8 

-4.803900O 

0  . 

6.8039000 

.272465 

1  .00000 

0 

-7.4289008 

0. 

7.4289800 

.291140 

1.00000 

0 

-8.83390OO 

0. 

8.0539000 

.38o311 

1.00000 

0 

-8.4789080 

0  . 

8.4789000 

.318770 

1.00000 

0 

-9.3039000 

e . 

9.3039000 

.329107 

1 .00000 

0 

-9.9289800 

0 . 

9.9289000 

.337766 

1 .00000 

0 

-10.553900 

0 . 

10.553900 

.345084 

1  .00000 

0 

-11. 178908 

0. 

1 1 . 178900 

.351326 

1 . 00000 

0 

-11.803980 

0 . 

1  1.803900 

. 356685 

1 .00000 

0 

-12.428900 

0. 

12.428908 

. 36 1 32 1 

1.00000 

0 

-13.053900 

8. 

13.053900 

.365359 

1.00000 

0 

-13.478960 

0 . 

13.678900 

.368896 

1 . 00000 

0 

-14.303900 

8. 

14.303900 

.3720  12 

1.00000 

0 

-14.928900 

0  . 

14.928900 

.374771 

1 .00000 

0 

-15.053900 

e . 

15.053908 

.375285 

1  .00000 

0 

BOUNDARY 


SUB-BRANCH  STARTING  AT  X  »  -.22949E-01  Y  *>  -.  10880E-01 

CALCULATION  STEP  SIZE  *  .1258 

PRINTING  STEP  RI ZE  =  -258 
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LOCUS  real 

LOCUS  }h*G  DIST  TO  ORIGIN 

GAIN 

zeta.  c: 

-.23481923E-01  0. 

.234S1923E-0 1 

2.42176 

1.00000  3 

0.  0. 

.  I00O0600E-02 

3.28130 

1.00600  4 

TROUBLES 

NEAR  A  POLE. 

REDUCING 

THE  STEP  SIZE  MAY  HELP. 

SUB-BRANCH  STARTING 

AT  X  -  1.8765 

Y  ■ 

15.000 

CALCULATION  STEP  SIZE  * 

.  1250 

PRINTING 

STEP  SIZE 

.6250 

LOCUS  REAL 

LOCUS  1MAG 

DIST  TO  ORIGIN 

GAIN 

2  ETA 

CD 

1 .6882509 

14.4040  18 

14.502618 

.440 195 

.11641 

0 

1.5878167 

13.805634 

13.887738 

.4430  1  1 

.  10857 

0 

1.3350548 

13.2049 90 

13.2723 07 

.446258 

.  10059 

0 

1.1699938 

12.662184 

12.656379 

.450030 

.09244 

0 

1.0126614 

1 1.997316 

12.039978 

.454447 

.08411 

0 

.86308430 

1 1.390482 

1 1 .423134 

.459671 

.07556 

0 

.72120787 

10.781784 

10.805884 

.465914 

.06675 

0 

.58729672 

10. 171320 

10 . 188261 

.473470 

.05764 

0 

.46113464 

9.559190 1 

9.5703062 

.482744 

.04818 

0 

.34282477 

8.9454942 

8.9520610 

.494321 

.03830 

0 

. 23 238997 

8.3303324 

8.3335733 

.509066 

.02789 

0 

.12985336 

7.7:33050 

7.7148979 

. 5283  1  1 

.01683 

0 

•35230873E-01 

7.0960 122 

7.0960997 

.554211 

.06497 

0 

0. 

6.8512476 

6.8512476 

.567033 

'.00000 

4 

-.8351554 3E-0 1 

6.2318568 

6.2324164 

.609024 

.01340 

0 

-.1590322 9 

5.61 14400 

5.6136931 

.672690 

.02833 

0 

-.22650526 

4.9900970 

4.9952350 

.778263 

.04534 

0 

-.28585856 

4 . 3o79260 

4.3772703 

.981052 

.06531 

0 

-.33695534 

3 . 7453227 

3.7601508 

1 ■ 50  366 

.08961 

0 

-.37950859 

3. 1214788 

3. 1444637 

5.27637 

. 12069 

0 

-.39156080 

2.9159000 

2.9420727 

0. 

. 13309 

2 

c 
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41 


M  =  0.4 
h  =  2 OK  ft 


OPEN-LOOP  ( OLTF)  ROOT  LOCUS  USING  OPTION  41 
CODE 

0-LOCUS  PT.  2-ZERO  4- I MAG I NARY  AXIS 

1-POLE  3-BREAK  PT .  5-SENSITIUITY  PT. 

REGION  OF  CALCULATION-REAL:  CC=  -15.0  TO  AA=  10.0 

IMAJ:  00=  —.100 E— 0 1  TO  BB=  15.0 

BRANCH  STARTING  AT  <-.012824)  ♦  J(  . 11286) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  *  TO  ABORT  >S 

BRANCH  STARTING  AT  (-.0  12824)  ♦  JC-.  11286) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  *  TO  ABORT  >S 

BRANCH  STARTING  AT  (1.8793)  ♦  J<0.) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  *  TO  ABORT  >S 

BRANCH  STARTING  AT  (-2.6203)  ♦  J(0.) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  *  TO  ABORT  >L 


BRANCH  NUMBER  4 


CALCULATION  STEP  SI2E  ■  .1250 

PRINTING  STEP  SIZE  «  .6250 


LOCUS  REAL 

LOCUS  IMAG 

DIST  TO  ORIGIN 

GAIN 

zeta 

CD 

-2.6203000 

0. 

2.6203000 

0. 

1.00000 

1 

-3.2453000 

0. 

3.2453000 

.  100248 

1.00000 

0 

-3.8703000 

0. 

3.8703000 

.  179688 

1.00000 

0 

-4.4953000 

0  . 

4.4953000 

.241082 

1.00000 

0 

-5. 1203000 

0  . 

5.1203000 

.288335 

1.00000 

0 

-5.7453000 

0. 

5.7453000 

.324922 

1.00000 

0 

-6.3703000 

0. 

6.3703000 

.353551 

1.00000 

0 

-6.9953000 

0  . 

6.9953000 

.376229 

1.00000 

0 

-7.6203000 

0. 

7.6203000 

.394423 

1.00000 

0 

-8.2453000 

0. 

8.2453000 

.409198 

1  .00000 

0 

-8.8703030 

0 . 

8.876 7000 

.421338 

1.00000 

0 

-9.4953000 

0  . 

9.495J800 

.431419 

1.00000 

0 

-10. 120300 

0  . 

10  .  120300 

.439S75 

1.00000 

0 

-10.745300 

0. 

10.745380 

.447033 

1.00000 

0 

-11.370300 

0. 

11.370300 

.453144 

1.00000 

0 

-11.995300 

0. 

11.995300 

.45840  1 

1.00000 

0 

-12.620300 

0 . 

12.620300 

.462955 

1.00000 

0 

-13.245300 

0 . 

13.245300 

.466926 

1.08000 

0 

-13.870300 

0. 

13.870300 

.470411 

1 .00000 

0 

-14.495300 

0  . 

14.495300 

.473484 

1.00000 

0 

-15.  120300 

0. 

15. 128300 

.476210 

1.00000 

0 

60UN0ARY 


SUB-BRANCH  STARTING  AT  X  *  1.8867  Y  -  15.000 
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''AD-A137  S14 
UNCLASSIFIED 


DESIGN  OF  LONGITUDINAL  CONTROL  LAWS  FOR  THE  X-29A  2/2 

BACKUP  ANALOG  FLIGHT  CO.  .  <U)  AIR  FORCE  INST  OF  TECH 
WRIGHT-PRTTERSON  AFB  OH  SCHOOL  OF  ENG I.  .  H  L  EMRICK 
SEP  83  AFIT/GAE/AA/83S-3  F/G  1/3 


NL 


CALCULATION  STEP  SIZE 
PRINTING  STEP  SIZE 


i  a 

-N .*•  •  .V 


LOCUS  REAL 
1.7088453 
1.5384859 
1.3753989 
1.2194073 
1.07 11328 
.92999572 
.79421494 
.44980774 
.55078943 
.43917385 
.33497109 
.23818841 
. 14882895 
.448878 10 E-0 1 
0. 

-.47958427E-01 

-. 12859854 

-.18204014 

-.22852090 

-.24854131 

-.27324000 


LOCUS  IMAG 
14.400851 
13.799521 
13. 194177 
12.5 90909 
11.983805 
11.374953 
10.744442 
10. 152343 
9.5388032 
8.9238541 
8.3074054 
7.4901482 
7.0715729 
4.4519713 
5.6980823 
5.2747914 
4.4547437 
4.0320340 
3.4087497 
2.7850559 
2.7054000 


.  1250 

.4250 

01  ST  TO  ORIGIN 
14.501885 
13.8850  18 
13.247441 
12.449840 
12.031580 
1 1.412907 
10.793849 
10.174434 
9.5544918 
8.9344542 
8.3143540 
7.4938341 
7.0731389 
4.4523180 
5.8980823 
5.2772291 
4.4545198 
4.0341433 
3.4144210 
2.7979744 
2.7191453 


GAIN 
.542840 
.545744 
.549133 
.553041 
.557414 
.543013 
.549455 
.577233 
.584754 
.598405 
.413434 
.433142 
.459204 
.495323 
.741144 
.818872 
.951248 
1.21742 
1 . 9785 1 
14.9320 
0. 


ZETA  CD 
.11784  0 

.11080  0 
.  10347  e 
.09441  0 

.08903  0 

.08149  0 

.07377  0 

.04583  0 

.05745  0 

.04915  0 

.04029  0 

.03094  0 

.02104  0 

.01037  6 

.00000  4 

.01288  0 
.02742  0 

.04510  0 

.04489  0 

.09598  0 

. 10049  2 


8  c® 


GAIN*- 1,41 


M  =  0.5 

h  *  40K  ft 


OPEN-LOOP  < OLTF)  ROOT  LOCUS  USING  OPTION  41 
CODE 

0 -LOCUS  PT.  .  -2-ZERO  4-IMA6INARY  AXIS 

1-POLE  3-BREAK  PT.  3-SENSITIVITY  PT. 

REGION  OF  CALCULATION-REAL!  CC-  -13.0  TO  AA->  10.0 

IMAJi  DD»  -.100E-01  TO  B0»  15.0 

BRANCH  STARTING  AT  (-.013728)  ♦  J(. 094631) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  *  TO  ABORT  >S 

BRANCH  STARTING  AT  (-.015728)  ♦  J(-. 094651) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  *  TO  ABORT  >S 

BRANCH  STARTING  AT  (1.578)  +  J(0.) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  •  TO  ABORT  >S 

BRANCH  STARTING  AT  (-1.9912)  +  J(0.) 

.  TYPE  L  TO  LIST,  S  TO  SKIP,  OR  *  TO  ABORT  >L 


BRANCH  NUMBER  4 


CALCULATION  STEP  SIZE  »  .1250 

PRINTING  STEP  SIZE  ■  .6250 


LOCUS  REAL 

LOCUS  I  NAG 

DIST  TO  ORIGIN 

GAIN 

ZETA 

CD 

-1.9912000 

0. 

1.9912000 

0. 

1.00000 

1 

-2.6162000 

0. 

2.6162000 

. 192622 

1.00000 

0 

-3.2412000 

0. 

3.2412000 

.337893 

1.00000 

0 

-3.8662000 

0. 

3.8662000 

.443711 

1.00000 

0 

-4.4912000 

0. 

4.4912000 

.520722 

1.00000 

0 

-3. 1 162000 

0. 

5.  1 162800 

.577490 

1.00000 

0 

-5.7412080 

0. 

5.7412000 

.620881 

1.00000 

0 

-6.3662080 

0. 

6.3662000 

.652636 

1.00000 

0 

-6.9912000 

0. 

6.9912000 

.677973 

1.00000 

0 

-7.6162000 

0. 

7.6162000 

.698822 

1.00000 

0 

-8.2412800 

0. 

8.2412000 

.714131 

1.00000 

0 

-8.8662080 

0. 

8.8662000 

.727253 

1.00000 

0 

-9.4912800 

0. 

9.4912000 

.738074 

1.00000 

0 

-10. 1 16200 

0 

10. 1 16200 

.747100 

1.00000 

0 

-10.741208 

0. 

10.741200 

.754704 

1.00000 

0 

-1 1 .366200 

0. 

11.366200 

.761 169 

1.00000 

0 

-11.991200 

0. 

1 1.991200 

.76671 1 

1.00000 

0 

-12.616200 

0. 

12.616200 

.771498 

1.00000 

0 

-13.241200 

0. 

13.241200 

.775663 

1.00000 

0 

-13.866200 

0. 

13.866200 

.779308 

1.00000 

0 

-14.491200 

0. 

14.491200 

.782517 

1.00000 

0 

-15. 1 16200 

0  . 

15.  1  16200 

.785358 

1.00000 

0 

BOUNDARY 


SUB-BRANCH  STARTING  AT  X  -  1.7188  Y  ■  15.000 


89 


(CD 


CALCULATION  STEP  SI2E 
PRINTING  STEP  SIZE 


LOCUS  REAL 
1.5444749 
1.4147427 
1.2754201 
1 .  1404418 
1.0119017 
.88975279 
.77402487 
.44473443 
.54188535 
.44548440 
.37554382 
.29205953 
.21503218 
.  14445407 
.80307429E-01 
.22558543E-01 
0. 

-.48749098E-01 
-.91415851E-01 
-.12830949 
-.  14033831 
-.  14858000 


LOCUS  IMAG 
14.394344 
13.787053 
13. 178247 
12.547995 
11.954343 
11.343418 
10.729228 
10. 113841 
9.4973844 
8.8798441 
8.2413745 
7.4419780 
7.0217455 
4.4007440 
5.7790492 
5. 1547255 
4.8927453 
4.2494735 
3.4441325 
3.0222243 
2.3980444 
2.2258000 


.  1250 

.4250 

DIST  TO  ORIGIN 
14.479113 
13.859455 
13.239822 
12.419434 
1 1 .99910  7 
11.378240 
10.757112 
10.135483 
9.5139911 
8.8920582 
8.2499058 
7.4475549 
7.0250372 
4.4023758 
5.7794071 
5. 1547748 
4.8927453 
4.2499520 
3 . 44727S3 
3.024 9447 
2.4034008’ 
2.2321749 


GAIN 
.854319 
.857347 
.840872 
.844931 
.849448 
.875245 
.881877 
.889852 
.899544 
.91 1575 
.924485 
.944894 
.971472 
1.00445 
1.05572 
1.12941 
l. 17244 
1.32334 
1.42031 
2.41511 
9. 18034 
0. 


ZETA 
.  10805 
.  10222 
.09433 
.09037 
.08433 
.07820 
.07195 
.04558 
.05984 
.05235 
.04541 
.03819 
.03041 
.02254 
.01389 
.00437 
.00000 
.01142 
.02504 
.04242 
.04471 
.07552 


Cl 


GAIN— l 
41 


M  *  0.6 

h  *  0  ft 


OPTION  > 


OPEN-LOOP  (OLTF)  ROOT  LOCUS  USING  OPTION  41 
CODE 

0 -LOCUS  PT.  2-ZERO  4- I MAG I NARY  AXIS 

1-POLE  3-BREAK  PT.  5-SENSITIVITY  PT. 

REGION  OF  CALCULATION-REAL:  CC-  -15.0  TO  AA-  10.0 

IMAJ:  DD-  -. 100E-8 1  TO  BB-  15.0 

BRANCH  STARTING  AT  (-.010817)  ♦  JC . 04491) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  *  TO  ABORT  >S 

BRANCH  STARTING  AT  (-.010817)  ♦  J(-. 04491) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  *  TO  ABORT  >S 

BRANCH  STARTING  AT  (4.0251)  ♦  J(0.) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  *  TO  ABORT  >S 

BRANCH  STARTING  AT  (-4.4274)  ♦  J(0.) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  *  TO  ABORT  >L 


BRANCH  NUMBER  4 


CALCULATION  STEP  SIZE  -  .1250 

PRINTING  STEP  SIZE  «  .4250 


LOCUS  REAL 

LOCUS  IMAG 

DIST  TO  ORIGIN 

GAIN 

ZETA 

-4.6274000 

0. 

6.6276000 

0. 

1.00000 

-7.2526000 

0. 

7.2526000 

. 785703E-02 

1.00000 

-7.8776000 

0. 

7.8776000 

.  15O327E-0 1 

1.00000 

-8.5026000 

0. 

8.5026000 

.2 15300E-0 1 

1.00000 

-9. 1276000 

0. 

9. 1276000 

.273S04E-0 1 

1.00000 

-9.7526000 

0. 

9.7526000 

. 326305E-0 1 

1.00000 

-10.377600 

0. 

10.377600 

. 373339E-0 1 

1.00000 

-11.002600 

0. 

11.002600 

. 4  15454E-0 1 

1.00000 

-11.627600 

0. 

11.627600 

.453178E— 0 1 

1.00000 

-12.252600 

0. 

12.252600 

.487003E-0 1 

1 .00000 

-12.877600 

0. 

12.877600 

. 5 17377E-0 1 

1.00000 

-13.502600 

0. 

13.502600 

. 544702E-0 1 

1.00000 

-14. 127600 

0. 

14.  127600 

. 569332E-0 1 

1.00000 

-14.752600 

0. 

14.752600 

. 59 1 582E-0 1 

1 .00000 

-15.002600 

0. 

15.002600 

. 59987SE-0 1 

1.00000 

BOUNDARY 


SUB-BRANCH  STARTING  AT  X  -  -.97454E-02  Y  -  -.100600-01 

CALCULATION  STEP  SIZE  -  .1250 

PRINTING  STEP  SIZE  -  .4250 

LOCUS  REAL  LOCUS  IMAG  DIST  TO  ORIGIN  GAIN  ZETA 

TROUBLE  IN  FINDING  THE  NEXT  POINT  TO  WITHIN  10E-5  ACCU'ACY 
REDUCING  THE  STEP  SIZE  MAY  HELP. 


SUB-BRANCH  STARTING  AT  X  -  .32273  Y  -  15.000 

CALCULATION  STEP  SIZE  *  .1250 

PRINTING  STEP  SIZE  -  .6230 


LOCUS  REAL 

LOCUS  IMAG 

DIST  TO  ORIGIN 

GAIN 

ZETA 

20565840 

14.386870 

14.387540 

.  127309 

.0 1429 

9372450 8E-01 

13.771177 

13.771493 

.  131615 

.00681 

13.233617 

13.233 617 

. 136079 

.00000 

10277008 

12.617126 

12.617545 

.  142234 

.00815 

20063426 

11.999837 

12.001514 

.  149864 

.01672 

29338946 

11.381790 

11.385576 

.  159331 

.62579 

38163014 

10.763024 

10.769787 

. 172114 

.03544 

46475117 

10.143577 

10.154218 

.  189069 

.04577 

54294786 

9.5234897 

9.5389543 

.213003 

.65692 

61621598 

8.9028007 

8.9241012 

.249063 

.06905 

68455185 

8.2815494 

8.3097937 

.309037 

.06238 

74795239 

7.6597750 

7.6962059 

.4271 17 

.09718 

80641522 

7.0375169 

7.0835690 

.761090 

. 11384 

85993874 

6.4148145 

6.4721974 

7.37146 

. 13287 

86537000 

6.3482060 

6.4069110 

0. 

. 13507 

GAIN*- 1 ,41 


M  =  0.6 
h  «  20K  ft 


OPEN-LOOP  <OLTF>  ROOT  LOCUS  USING  OPTION  41 
CODE 

0-LOCUS  PT.  2-ZERO  4- I MAG I NARY  AXIS 

l-POLE  3-BREAK  PT.  ,  5-SENSITIUITY  PT. 

.  REGION  OF  CALCULATION-REAL i  CC-  -13.0  TO  AA-  10.0 

IMAJt  DD»  -.  100E-0 1  TO  BB>  15.0 

BRANCH  STARTING  AT  <-.0070548)  +  J< . 074164) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  *  TO  ABORT  >S 

BRANCH  STARTING  AT  (-.0070548)  ♦  J<-. 076144) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  *  TO  ABORT  >S 

BRANCH  STARTING  AT  (2.8673)  ♦  J<0.) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  *  TO  ABORT  >S 

BRANCH  STARTING  AT  (-4.0371)  ♦  J<0.) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  *  TO  ABORT  >L 


BRANCH  NUMBER  4 


CALCULATION  STEP  SIZE  -  .1230 

PRINTING  STEP  SIZE  -  .6250 


LOCUS  REAL 

LOCUS  IMAG 

DIST  TO  ORIGIN 

GAIN 

ZETA 

CD 

—4 .0571000 

0. 

4.0591000 

0. 

1.00000 

1 

-4.6841000 

0. 

4.6841000 

. 283670E-0 1 

1 .00000 

0 

-5.3091000 

0. 

5.3091000 

. 530369E-0 1 

1.00000 

0 

-5.9341000 

0. 

5.9341000 

. 74 126 1E-0  1 

1.00000 

0 

-6.5391000 

0. 

6.5591000 

. 9200 19E-0 1 

1.00000 

0 

-7. 1841000 

0. 

7. 1841000 

. 107109 

1.00000 

0 

-7.8091000 

0. 

7.8091000 

.  1 19883 

1.00000 

0 

-8.4341000 

0. 

8.4341000 

.  130713 

1.00000 

0 

-9.0591000 

0. 

9.0591000 

.  139931 

1 .00000 

0 

-9.6841000 

0. 

9.6841000 

.  147313 

1.00000 

0 

-10.309100 

0. 

10.309100 

.  154591 

1 .00000 

0 

-10.934100 

0. 

10 .934100 

.160446 

1.00000 

0 

-11.559100 

0. 

11.559100 

. 165532 

1.00000 

0 

-12. 184100 

0. 

12. 184100 

. 16997 1 

1 .00000 

0 

-12.809100 

0. 

12.809100 

. 173865 

1.00000 

0 

-13.434100 

0. 

13.434100 

. 177297 

1.00000 

0 

-14.039100 

0. 

14.059100 

.  130336 

1.00000 

0 

-14.684100 

0. 

14.684100 

.  183038 

1 .00000 

0 

-15.039100 

0. 

15.059100 

.  184518 

1.00000 

0 

BOUNDARY 


SUB-BRANCH  STARTING  AT  X  -  -.14648E-02  Y  -  -.1O000E-01 

CALCULATION  STEP  SIZE  -  .1250 

PRINTING  STEP  SIZE  -  .6230 

LOCUS  REAL  LOCUS  IMAG  DIST  TO  ORIGIN  GAIN  ZETA  CD 

TROUP’  r  IN  FINDING  THE  NEXT  POINT  TO  WITHIN  10E-5  ACCURACY. 


93 


REDUCING  THE  STEP' SI2E  MAT  HELP 


SUB-BRANCH  STARTING  AT  X  -  .87138  Y  -  15.600 


CALCULATION  STEP  SIZE  -  .1230 

PRINTING  STEP  SI2E  -  .4250 


LOCUS  REAL 

LOCUS  imG 

DIST  TO  ORIGIN 

GAIN 

2  ETA 

CD 

75814272 

14.385380 

14.405345 

.251421 

.05243 

0 

44982834 

13.749843 

13.785147. 

.255088 

.04714 

0 

54423875 

13. 153490 

13.144827 

.259135 

.04149 

0 

44737558 

12.534359 

12.544339 

.243902 

.03544 

0 

35324783 

11.918484 

11.923720 

.249578 

.02943 

0 

24391248 

11.299908 

11.302989 

.274422 

.02335 

0 

17931345 

10.480441 

10 .482144 

.284795 

.01479 

0 

99474938E-0 1 

10.040783 

10.041275 

.295214 

.00989 

0 

24399 187E-01 

9.440310 1 

9.4403417 

.308440 

.00258 

6 

9.2294142 

9.2294142 

.313792 

.00000 

4 

48489423E-0 1 

8.4084038 

8.4084778 

.332824 

.00798 

0 

13241883 

7 . 9844833 

7.9877843 

.358818 

.0 1440 

0 

19178307 

7.3444915 

7.3449883 

.394023 

.02403 

0 

24418948 

4.7418455 

4.7443590 

.453018 

.03449 

0 

29584704 

4.1188428 

4. 1259908 

.549722 

.04829 

0 

34077238 

5.4954410 

5.5046144 

.745527 

.04169 

6 

38099543 

4.8717580 

4.8844332 

1.33179 

.07797 

0 

41454983 

4.2477724 

4.2481498 

48. 1548 

.09740 

0 

41742000 

4.2318000 

4.2523371 

0. 

.09814 

2 

'  .v.  ir.%  :r. 


GAIN— 1,41 


M  «-0.6 

h  =  40K  ft 


OPEN-LOOP  < OLTF)  ROOT  LOCUS  USING  OPTION  41 
CODE 

0-LOCUS  PT.  2-ZERO  4- I MAG I NARY  AXIS 

1-POLE  3-BREAK  PT..  5-SENSITIUITY  PT. 

REGION  OF  CALCULATION-REAL i  CC-  -15.0  TO  AA»  10.0 

IMAJj  DO  -.  100E-0 1  TO  BB»  15.0 


BRANCH 

STARTING 

AT 

<-.010247) 

♦ 

J< .080888) 

TYPE  L 

TO  LIST, 

S 

TO  SKIP,  OR 

% 

TO  ABORT 

>S 

BRANCH 

STARTING 

AT 

<-.010247) 

♦ 

J<-. 080888) 

TYPE  L 

TO  LIST, 

S 

TO  SKIP,  OR 

s 

TO  ABORT 

>s 

BRANCH 

STARTING 

AT 

<1.9091)  ♦ 

J<0.) 

TYPE  L 

TO  LIST, 

S 

TO  SKIP,  OR 

* 

TO  ABORT 

>s 

BRANCH 

STARTING 

AT 

<-2.4136)  •• 

-  J<0.) 

TYPE  L 

TO  LIST, 

S 

TO  SKIP,  OR 

* 

TO  ABORT 

>L 

BRANCH  NUMBER  4 


CALCULATION  STEP  SIZE  «  .1250 

PRINTING  STEP  SIZE  -  .6250 


LOCUS  REAL 

LOCUS  IMAG 

DIST  TO  ORIGIN 

GAIN 

ZETA 

CD 

-2.4136000 

0. 

2.4136000 

0. 

1.00000 

I 

-3.0386000 

0. 

3.0386000 

. 1 12034 

1.00000 

0 

-3.6636000 

0. 

3.6636000 

.281191 

1.00000 

0 

-4.2886000 

0. 

4.2886808 

.2700 14 

1.00000 

0 

-4.9136000 

0. 

4.9136000 

.322790 

1.00000 

6 

-5.5386000 

0. 

5.5386080 

.363459 

1.00000 

0 

-6. 1636000 

0. 

6. 1636000 

.395121 

1.00000 

0 

-6.7886000 

0. 

6.788600O 

.420076 

1 .00000 

0 

-7.4136000 

0. 

7.4136000 

" .439997 

1.00000 

0 

-8.0386000 

0. 

8.0386000 

.45610  1 

1.00000 

0 

-8.6636000 

0. 

8.6636000 

.469274 

1.00000 

0 

-9.2886000 

0. 

9.2886080 

.488  168 

1.00000 

0 

-9.9136000 

0. 

9.913 6000 

.489270 

1.00000 

0 

-10.538600 

0. 

10.538600 

.496947 

1.00000 

0 

-11. 163600 

0. 

1  1.  163600 

.503477 

1.00000 

0 

-11.788600 

0. 

1 1.788608 

.509076 

1 .00000 

0 

-12.413600 

0. 

12.413680 

.513911 

1.00000 

0 

-13.038600 

0. 

13.038600 

.513115 

1 .00000 

0 

-13.663600 

0. 

13.663600 

.521792 

1.00000 

0 

-14.288600 

0. 

14.288608 

.525028 

1 .00000 

c> 

-14.913600 

0. 

14.913608 

.527889 

1 .80000 

0 

-15.038600 

0. 

15.038600 

.528421 

1 .00000 

0 

BOUNDARY 


SUB-BRANCH  STARTING  AT  X  -  1.2607  Y  -  15.OO0 

CALCULATION  STEP  SIZE  -  .1250 

PRINTING  STEP  SIZE  -  .6250 

95 


LOCUS  REAL 

LOCUS  IMAG 

DIST  TO  ORIGIN 

GAIN 

2  ETA 

1 . 1398766 

14.386800 

14.43  1886 

.597081 

.07898 

1.0241511 

13.772609 

13.810636 

.600275 

.07416 

.91347338 

13. 157489 

13. 189160 

.603958 

.06926 

.80785628 

12.541480 

12.567472 

.608237 

.06428 

.70728805 

11.924625 

11.945582 

.613252 

.05921 

. 6 1 179240 

11.306965 

11.323505 

.619185 

.05403 

.52136835 

10.688543 

10.781251 

.626278 

.04872 

.43602014 

10.069399 

10.078835 

.634864 

.04326 

.35575101 

9.4495771 

9.4562712 

.645404 

.03762 

.28056295 

8.8291179 

8.8335745 

.658562 

.03176 

.21045628 

8.2080639 

8.2107616 

.675317 

.02563 

. 14542894 

7.5864577 

7.5878514 

.697170 

.01917 

•85475344E-01 

6.9643415 

6.9648660 

.726542 

.01227 

.38584498E-0 1 

6.3417582 

6.3418320 

.767567 

.00482 

• 

5.9670882 

5.9670882 

.800803 

.00000 

.  46830595E-0 1 

5.3438468 

5.3440520 

.879447 

.00876 

.88687694E-0 1 

4.7202515 

4.7210846 

1.01165 

.01879 

.  12566535 

4.8963478 

4.0982749 

1.27115 

.03066 

. 15795046 

3.4721836 

3.4757743 

1.97110 

.04544 

. 18595203 

2.8478122 

2.8538767 

8.66802 

.06516 

.  1921 1000 

2.6983000 

2.7051302 

0. 

.07102 

96 


y  j 


GAIN—  1 ,41 


M  *  0.8 

h  s  0  ft 


OPEN-LOOP  (OLTF)  ROOT  LOCUS  USING  OPTION  41 
CODE 

0 -LOCUS  PT.  2-ZERO  4- IMAGINARY  AXIS 

1-POLE  3-BREAK  PT.  5-SENSITIL'ITY  PT. 


REGION  OF  CALCULATION-REAL!  CC-  -10.0  TO  AA«  6.** 

IMAJ:  DD=  -12.0  TO  BB=  12. *3 

BRANCH  STARTING  AT  (-.014487)  ♦  J(. 04 1542) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  *  TO  ABORT  >S 

BRANCH  STARTING  AT  (-.014487)  ♦  J(-. 04 1542) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  *  TO  ABORT  >S 

BRANCH  STARTING  AT  (5.4774)  ♦  J<0.) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  *  TO  ABORT  >S 

BRANCH  STARTING  AT  (-9.3031)  +  J(0.) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  *  TO  ABORT  >L 


BRANCH  NUMBER  4 

CALCULATION  STEP  SI2E  »  . 12S0 

PRINTING  STEP  SIZE  -  .4400 


LOCUS  REAL 
-9.3031000  0. 

-9.9431000  0. 

-10.071100  0. 

BOUNDARY 


LOCUS  I MAG 


DIST  TO  ORIGIN 
9.3031000 
9.9431000 
10.071 100 


GAIN 

0. 

. 299934E-02 
. 357729E-02 


:eta 
1.00000 
•  .06000 
1.00000 


SUB-BRANCH  STARTING  AT  X  -  -.93900  Y  »  -12.000 


CALCULATION  STEP  SIZE  »  .1280 

PRINTING  STEP  SIZE  »  .4400 


LOCUS  REAL 
-1.0123030 
-1.0814220 
-1 .  1445542 
-1.2074941 
-1.2572000 


LOCUS  I MAG 
-11.344213 
- 10 . 72795S 
-10.091282 
-9.4542099 
-8.9047000 


DIST  TO  ORIGIN 
1 1.409211 
1 0 . 7S2324 
10 . 154208 
9.5310343 
8.9930 103 


GAIN 
. 172771 
.222210 
.325321 
.449008 
0. 


ZETA 
.08873 
. 1O030 
. 1 1289 
.12471 
. 13980 


SUB-BRANCH  STARTING  AT  X  »  -.93900  Y  ■ 

CALCULATION  STEP  SIZE  »  .12S0 

PRINTING  STEP  SIZE  *  .4400 

LOCUS  REAL  LOCUS  I MAG  DIST  TO  ORIGIN  GAIN 

_  <  ii  li.nu  it  ,10  9*7  1  t  i  i 


12.0  •)  O 


ZETA 

0  O  Q  —  7 
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GAIN— 1,41 


M  -  1.0 
h  »  20K  ft 


OPEN-LOOP  < OLTF)  ROOT  LOCUS  USING  OPTION  41 
COOE 

0-LOCUS  PT.  2-2ER0  4- I MAG I NARY  AXIS 

1-POLE  3-BREAK  PT.  5-SENSITIUITY  PT. 

REGION  OF  CALCULATION-REAL:  CC-  -15.0  TO  AA-  10.0 

IMAJ:  DO  — .  100E-0 1  TO  BB-  15.0 

BRANCH  STARTING  AT  (-.03975)  ♦  J<0.) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  •  TO  ABORT  >S 

BRANCH  STARTING  AT  (.816974)  ♦  J<0.) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  *  TO  ABORT  >S 

BRANCH  STARTING  AT  (3.9644)  ♦  J(0.) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  *  TO  ABORT  >S 

BRANCH  STARTING  AT  (-6.0004)  ♦  J(0.) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  *  TO  ABORT  >L 


BRANCH  NUMBER  4 
CALCULATION  STEP  SIZE  «  .1250 


PRINTING 

STEP  SIZE 

.6250 

LOCUS  REAL 

LOCUS  I MAG 

DIST  TO  ORIGIN 

GAIN 

ETTA 

-6.0004000 

0. 

6.0004000 

0. 

1.00000 

-6.6254000 

0. 

6.6254000 

•630787E-02 

1.00000 

-7.2504000 

0. 

7.2504000 

.  120895E-0 1 

1  .00000 

-7.8754000 

0. 

7.8754000 

. 173338E-0 1 

1.00000 

-8.5004000 

0. 

B. 5004000 

•220579E-0 1 

1.00000 

-9. 1254000 

0. 

9. 1254000 

. 262952E-0 1 

1 .00000 

-9.7504000 

0. 

9.7504000 

. 30087 1E-0 1 

1.00000 

-10.375400 

0. 

10 .375400 

. 334773E-0 1 

1.00000 

-11.000400 

0. 

11.000400 

. 365085E— 0 1 

1.00000 

-11.625400 

0. 

11.625408 

. 39221 2E-0 1 

1.00000 

-12.250400 

0. 

12.250400 

. 4 16522E-0 1 

1.00000 

-12.875400 

0. 

12.875400 

.43S345E-0 1 

1.00000 

-13.500400 

0. 

13.500400 

. 457977E-0 1 

1.00000 

-14. 125400 

0. 

14. 125400 

. 475675E-0 1 

1.00000 

-14.750400 

0. 

14.750400 

. 49  1667E-0 1 

1.00000 

-15.000400 

0. 

15.000400 

. 49763 1E-0 1 

1.00000 

BOUNDARY 

SUB-BRANCH  STARTING  AT 

X  »  .65527 

Y  - 

15.000 

CALCULATION  STEP  SIZE  ■ 

.  1250 

PRINTING 

STEP  SIZE 

.6250 

LOCUS  REAL 

LOCUS  I  MAG 

DIST  TO  ORIGIN 

GAIN 

ZETA 

.53325417 

14.387029 

14.396908 

.  989728E-0 l 

.03704 

.41642920 

'  13.773046 

13.779348 

.  182104 

.03022 

.30470242 

13. 158115 

13. 161643 

. 105860 

.02315 

.  19808153 

12.542279 

12.543843 

. 1 10433 

.81579 

zc 


."9337384 1E-0  1 
. 1863 1843E-03 
6. 

-  •  9 1250569E-0  1 

-.  17736810 

-.25834673 

-.33418096 

-.40486569 

-.47039622 

-.53076832 

-.58597825 

-.59373000 


T1T925578 

11.308057 

11.306830 

10.688529 

10.069492 

9.4497622 

8.8293817 

8.2083933 

7.5868400 

6.9647644 

6.3422094 

6.2498000 


11.925969 
1 1.308057 
1 1.306830 
10.688919 
10.071054 
9.4532930 
8.8357036 
8.2183719 
7.60 14086 
6.9849594 
6.3692221 
6.2779388 


. 1 16099 
. 123272 
.  123288 
.  132620 
.  145179 
.  162873 
.  189453 
.233450 
.319288 
.556509 
4.00809  . 

0. 


.00810  0 
.00002  0 
.00000  4 

.00854  0 

.01761  0 

.02733  0 

.03782  6 

.04926  0 

.06188  0 
.07599  0 

.09200  0 

.09457  2 


Ic® 


»  «  *»  -  •  -  *  -  »  » 


GAIN— 1,41 


M  «  1.0 

h  «  40K  ft 


OPEN-LOOP  <OLTF>  ROOT  LOCUS  USING  OPTION  41 
CODE 

•-LOCUS  PT.  2-ZERO  4- I MAG I NARY  AXIS 

1-POLE  3-BREAK  PT.  3-SENSITIVITY  PT. 

REGION  OF  CALCULATION-REAL*  CC«  -13.0  TO  AA-  10.0 

IMAJt  DD*  -. 100E-0 1  TO  0B-  13.0 

BRANCH  STARTING  AT  <-.0047527)  ♦  J< . 047972) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  *  TO  ABORT  >S 

BRANCH  STARTING  AT  <-.0047327)  J<-. 047972) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  •  TO  ABORT  >S 

BRANCH  STARTING  AT  <2.4941)  ♦  J<0.) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  *  TO  ABORT  >S 

BRANCH  STARTING  AT  < -3. 5394)  +  J<0.) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  *  TO  ABORT ,>L 


BRANCH'  NUMBER  4 

CALCULATION  STEP  SIZE  -  .1230 

PRINTING  STEP  SIZE  -  .4250 


LOCUS  REAL 

LOCUS  IMAG 

DIST  TO  ORIGIN 

GAIN 

ZETA 

CD 

-3.3394008 

0. 

•  3.5394000 

0. 

1.00000 

1 

-4. 1444000 

0. 

4.1444000 

. 258845E— 0 1 

1.00000 

0 

-4.7894000 

0. 

4.7894000 

•483594E-0 1 

1.00000 

0 

-3.4144000 

0. 

5.4144000 

. 47943  IE-01 

1.00000 

e 

-4.0394000 

0. 

4.0394000 

•843709E— 0 1 

1.00000 

0 

-4.4444000 

0. 

4.4444000 

.981779E-0 1 

1.00000 

0 

-7.2894000 

0. 

7.2894000 

.  109791 

1.00000 

0 

-7.9144000 

0. 

7.9144000 

. 1 19580 

1.00000 

0 

-8.3394000 

0. 

8.5394000 

. 127843 

1.00000 

0 

-9. 1444000 

0. 

9.  1444080 

. 134910 

1.00000 

0 

-9.7894000 

0. 

9.7894000 

. 140933 

1.00000 

0 

-10.414400 

0. 

10.414400 

. 144116 

1.00000 

0 

-11.039400 

0. 

11.039400 

.150585 

1.00000 

0 

-11.444400 

0. 

11.444400 

. 134473 

1.00000 

0 

-12.289400 

0. 

12.289400 

. 157870 

1.00000 

0 

-12.914400 

0. 

12.914400 

. 140852 

1.00000 

0 

-13.539400 

0. 

13.539400 

. 143482 

1.00000 

0 

-14.144400 

0. 

14.144400 

.  145812 

1.00000 

0 

-14.789400 

0. 

14.789400 

.  147884 

1  .00000 

0 

-15.039400 

0. 

13.039400 

.  148452 

1.00000 

0 

BOUNDARY 


SUB-BRANCH  STARTING  AT  X  -  -.34180E-02  Y-  -.10000E-01 

CALCULATION  STEP  SIZE  -  .1250 

PRINTING  STEP  SIZE  -  .4250 

LOCUS  REAL  LOCUS  I MAT  njpr  TQ  ORIGIN  GAIN  ZETA  CD 


101 


SUB-BRANCH  STARTING  AT  X  -  .97900  Y  -  15.000 

CALCULATION  STEP  SIZE  -  .1250 

PRINTING  STEP  SIZE  -  .6250 


LOCUS  REAL  LOCUS  I MAG 
.86964982  14.3B4644 

.76523054  13.768438 

. 66538350  13.151458 

.57011393  12.533763 

.47942675  11.915379 

.39332653  11.296340 

.31181746  18.676679 

.23490332  10.056431 

.16258738  9.4356299 

.94872308E-0 1  8.8143104 

.31760032E-01  8.1925064 

•  7.8609146 

.5605161 1E-0 1  7.2384344 

.10750463  6.6155573 

.15435995  5.9923175 

.  19662648  5.3687497 

.23432109  4.7448888 

.26742000  4 . 1219000 


DIST  TO  ORIGIN 

GAIN 

14.410908 

.221735 

13.789679 

.224482 

13.168280 

.227685 

12.546723 

.231454 

11.925020 

.235935 

11.363185 

.241330 

10.681231 

.247916 

10.059174 

.256094 

9.4370306 

.266456 

8.8148209 

.279914 

8.1925680 

.297953 

7.8609146 

.310255 

7.2386515 

.341009 

6.6164308 

.388255 

5.9943053  • 

.468770 

5.3723491 

.633019 

4.750671 1 

1.13399 

4.  1305657 

0. 

ZETA  CD 
.06035  0 
.05549  0 
.05053  0 
.64544  0 
.04020  0 
.03480  0 
.02919  6 
.02335  6 
.01723  8 
.01076  0 
.00388  6 
.00000  4 
.00774  6 
.01625  0 
.02575  0 
.03660  0 
.04932  0 
.06474  2 
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GAIN— 1,41 


M  «  1.2 
h  «  20K  ft 


OPEN-LOOP  < OLTF)  ROOT  LOCUS  USING  OPTION  41 
CODE 

8-LOCUS  PT.  2-ZERO  4- I MAG I NARY  «XIS 

1-POLE  3-BREAK  PT.  3-SENSITIVITY  PT. 

REGION  OF  CALCULATION-REAL!  CC«  -18.0  TO  AA-  6. 00 

IMAJs  DD»  -12.0  TO  BB-  12.0 

BRANCH  STARTING  AT  (-.8? 1117)  ♦  J<0:> 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  *  TO  ABORT  >S 

BRANCH  STARTING  AT  <.033264)  +  J<8.) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  *  TO  ABORT  >S 

BRANCH  STARTING  AT  <3.4147)  ♦  J<8.) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  *  TO  ABORT  >S 

BRANCH  STARTING  AT  < -3. 7362)  ♦  J<0.> 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  *  TO  ABORT  >L 


BRANCH 

NUMBER  4 

CALCULATION  STEP  SIZE  » 

.  1288 

PRINTING 

STEP  SIZE  « 

.6400 

LOCUS  REAL 

LOCUS  IMAG 

DI ST  TO  ORIGIN 

GAIN 

ZETA 

5.7362800 

0. 

5.7362000 

0. 

1.00000 

6.3762000 

0. 

6.3762000 

.40 1956E-02 

1.00000 

7.8162000 

0. 

7.0162000 

. 775702E-02 

1 . 00000 

7.6562000 

0. 

7.6562000 

.11184 1E-0 1 

1  .00000 

8.2962000 

0. 

8.2962000 

.  142970E-0 1 

1.00000 

8.9362008 

0. 

8.9362000 

.  171069E-0 1 

1.00000 

9.3762000 

0. 

9.5762000 

.  196337E-0 1 

1 .00000 

10.088200 

0. 

10.088200 

. 2 14674E-0 1 

1.00000 

BOUNDARY 

SUB-BRANCH  STARTING  AT 

X  -  .35800 

Y  -  -12.000 

CALCULATION  STEP  SIZE  - 

.  1280 

PRINTING 

STEP  SIZE  - 

.6400 

LOCUS  REAL 

LOCUS  IMAG 

DIST  TO  ORIGIN 

GAIN  ZETA 

.22127242 

-1 1.374780 

11.376932 

. 8 18842E-0 1  .01945 

.920  18295E-0 l 

-10.747972 

10 .748365 

.SS2560E-01  .O0356 

• 

-10.277548 

10.277543 

. 944845E-0 1  .00000 

. 1 1643212 

-9.6482320 

9.6489345 

.105814  .01207 

.22352421 

-9.0176018 

9.0204215 

.122943  .02500 

.32725949 

-8.3857431 

8.3921264 

.151585  .03900 

.42163005 

-7.7527425  ' 

7.7641991 

. 20S528  .05430 

.50863180 

-7. 1 186871 

7.  1368349 

.373609  .07127 

.38826625 

■*■6 .4836643 

6.5102964 

5.59816  .09036 

.59323080 

-6.4420000 

6.4692570 

0.  .09170 

SUB-BRANCH  STARTING  AT  X 


35800 


Y  - 


12.000 


CALCULATION  STEP  SIZE  -  .1280 

PRINTING  STEP  SIZE  -  .6400 


LOCUS  REAL  1  LOCUS  I MAG 
.22127242  11.374780 

.92018295E-01  10.747972 

18.277548 
.11643212  9.6482320 

.22552421  9.0176018 

.32725949  8.3857431 

.42163005  7.7527425 

.50863188  7.1186871 

.58826625  6.4836643 

.59323000  6.4420000 


01  ST  TO  ORIGIN  GAIN 

11.376932  •  .818842E-01 

10.748365  .082560E-01 

18.277548  .944845E-81 

9.6489345  .105814 

9.0204215  .122943 

8.3921264  .151585 

7.7641991  .208528 

7.1368349  .373609 

6.5102964  5.59816 

6.4692570  0. 


:eta  co 
.01943  0 

.00856  «. 

.00000 
.01207  i 
.02500  0 

.03900  0 

.05430  0 

.07127  0 

.0903.*  0 

.091-0  2 
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GAlNM-1,41 


M  ■  1.4 
h  ■  40K  ft 


OPEN-LOOP  <OLTF)  ROOT  LOCUS  US I NO  OPTION  41 
COOC 

•-LOCUS  PT.  2-2ERO  4- I MAG I NARY  AXIS 

1-POLE  3-BREAK  PT.  5-SENSITIVITY  PT. 

REGION  OP  CALCULAT I ON-REAL i  CO  -15.0  TO  AA-  10.0 

IMAJ :  00-  -.  100E-0 1  TO  BO  15.0 

•RANCH  STARTING  AT  <-.0082773)  ♦  J<  . 039491) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  «  TO  ABORT  >S 

•RANCH  STARTING  AT  (-.0082773)  ♦  J<-. 039491) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  *  TO  ABORT  >S 

•ranch  starting  at  < 1.5134)  *  j<0.) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  4  TO  ABORT  >S 

BRANCH  STARTING  AT  (-2.594)  ♦  J<0.) 

TYPE  L  TO  LIST,  S  TO  SKIP,  OR  »  TO  ABORT  >L 


BRANCH  NUMBER  4 


CALCULATION  STEP  SIZE  -  ,1250 

PRINTING  STEP  SIZE  •  .4250 


LOCUS  REAL 

LOCUS  I MAG 

01  ST  TO  ORIGIN 

GAIN 

ZETA- 

CD 

-2.5960000 

0. 

2.5960000 

0. 

1.00000 

l 

-3.2210000 

0. 

3.2210000 

. 927762E-02 

1.00000 

0 

-3.8460000 

0. 

3.8460080 

.  185849E-01 

1.00000 

0 

-4.4710000 

0. 

4.4710000 

•273683E-0 1 

1.00000 

0 

-5.0960000 

0. 

5.0960000 

.353435E-0  1 

1.00000 

0 

-5.7210000 

0. 

5.7210000 

. 424 133E-0  1 

1.00000 

0 

-6.3460000 

0. 

6.3460000 

. 4859 14E-0 1 

1.00000 

0 

-6.9710000 

0. 

*  6.9710000 

. 53948 IE-01 

1.00000 

0 

-7.5960000 

0. 

7.5960000 

. 585765E-0 1 

1.00000 

0 

-8.2210000 

0. 

8.2210000 

.625732E-0 1 

1.00000 

0 

-8.6440000 

0. 

8.8460000 

.660288E-0 1 

1.00000 

0 

-9.4710000 

0. 

9.4710000 

. 69024 1E-0 1 

1.00000 

0 

-10.094000 

0. 

10.096000 

. 7 16289E-0 1 

1.00000 

0 

-10.721000 

0. 

10.721000 

. 739025E-0 1 

1.00000 

0 

-11.346000 

0. 

11.346000 

. 758950E-0 1 

1.00000 

0 

-11.971000 

0. 

11.971000 

. 776482E-0 1 

1.00000 

0 

-12.596000 

0. 

12.596000 

. 79 1971E-0 1 

1.00000 

0 

-13.221000 

0. 

13.221000 

. 8057 1 1E-0 1 

1.00000 

0 

-13.846000 

0. 

13.646000 

. 8 17945E-0 1 

1.00000 

0 

-14.471000 

0. 

14.471000 

. 82888 1E-0  1 

1.00000 

0 

-15.096000 

0. 

15.096000 

. 83869 1E-0 1 

1 .00000 

0 

BOUNDARY 


SUB-BRANCH  STARTING  AT  X  »  -.78125E-02  Y  -  -.10G00E-01 

CALCULATION  STEP  SIZE  »  .1250 

PRINTING  STEP  SIZE  «  .6250 
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TROUBLE  IN  FINDING  THE  NEXT  POINT  TO  WITHIN  10E-5  ACCURACY. 
REDUCING  THE  STEP  SIZE  MAY  HELP. 


SUB-BRANCH  STARTING  AT  X  -  1.6631  Y  ■  15.000 


CALCULATION  STEP  SIZE  ■  .1250 

PRINTING  STEP  SIZE  «  .6250 


LOCUS  REAL 

LOCUS  I MAG 

DIST  TO  ORIGIN 

GAIN 

ZETA 

1.4897629 

14.399518 

14.476378 

.111569 

. 10291 

1.3236763 

13.796993 

13.860344 

.112999 

.09550 

1.1646991 

13.192554 

13.243867 

.  1 14674 

.08794 

1.6128532 

12.586284 

12.626972 

. 116656 

.88021 

.86815994 

11.978267 

12.009687 

.119029 

.07229 

.73063920 

11.368588 

11.392042 

.  121909 

.06414 

.60031010 

10.757331 

10.774068 

. 125460 

.05572 

.47719073 

10. 144581 

10. 155798 

.  129925 

.04699 

.36129814 

9.5304234 

9.5372694 

.  135673 

.03788 

.25264841 

8.9149432 

8.9185225 

. 143297 

.02833 

. 15125657 

8.2982258 

8.2996042 

.  153809 

.01822 

.57136678E-0 1 

7.6803568 

7.6805693 

. 169082 

.00744 

• 

7.2791211 

7.2791211  . 

.  183212 

.00000 

.82107510E-0 1 

6.6595414 

6.6600475 

.217200 

.81233 

.  15690786 

6.0390371 

6.0410752 

.285023 

.02597 

.22439375 

5.4176948 

5.4223398 

.481565 

.84138 

.28455585 

4.7956006 

4.8040355 

6.02775 

.85923 

.28872000 

4.7495000 

4.7582675 

0. 

.06068 

42 


ENTER  GAIN  OF  INTEREST  (GAIN) .TOLERANCE  (GTOL) 0-1 .2, . t 


OPEN-LOOP  < OLTF)  ROOT  LOCUS  USING  OPTION  42  M  «  0.3 

h  ■  0  ft 

CODE 

0 -LOCUS  PT.  2-2 ERO  4- I MAG I NARY  AXIS 

1-POLE  3-BREAK  PT.  5-SENSITIVITY  PT. 

GAIN  OF  INTEREST  (GAIN)  -  -1.260  ♦OR-  GTOL  -  .1006 


SENSITIVITY  OF  INTEREST  <OLK>  «  -2.923  +0R-  KTOL  «  .2436 


ROOTS  OF  INTEREST 
X  -  .43183E-02 

X  -  . 43183E-02 

X  -  -.29209 

X  -  -.29209 


Y  -  -.607 57E-01 

Y  -  . 60757E-0 1 

Y  -  -4.0147 

Y  ■  4.0147 


REGION  OF  CALCULATION-REAL:  CC-  -5.00  TO  AA-  1.03 

IMAJ:  DO-  -3.00  TO  BB-  3.03 

BRANCH  STARTING  AT  (-.014311)  ♦  J(.14132> 

42 

ENTER  GAIN  OF  INTEREST  (GAIN) .TOLERANCE  (GTOL) 0-1 .2, . 1 


OPEN-LOOP  (OLTF)  ROOT  LOCUS  USING  OPTION  42 
CODE 

0 -LOCUS  PT.  2-ZERO 

1-POLE  3-BREAK  PT. 


M  *  0.4 
h  *  20K  ft 

4—  IMAGINARY  AXIS 

5- SENSITIVITY  PT. 


GAIN  OF  INTEREST  (GAIN)  -  -1.280  +  0R-  GTOL  -  .1060 


SENSITIVITY  OF  INTEREST  (OLK)  -  -2.344  ♦OR-  KTOL  -  .1954 

ROOTS  OF  INTEREST 

X  -  . 46977E-02  Y  -  -.52395E-01 

X  -  .46977E-02  Y  -  .52395E-01 

X  -  -. 17979  Y  -  4.0600 

X-  -.17979  Y-  -4.0600 


REGION  OF  CALCULATION-REAL;  CC«  -5.80  TO  AA-  1.03 

IMAJ:  DD-  -3.00  TO  BB-  3.00 

BRANCH  STARTING  AT  (-.012824)  ♦  J( . 1 1286) 

42 

ENTER  GAIN  OF  INTEREST  (GAIN)  .TOLERANCE  ( GTOL)  s >- 1 . 2 , . 1 


OPEN-LOOP  (OLTF)  ROOT  LOCUS  USING  OPTION  42 
CODE 

0-LOCUS  PT.  2-ZERO 


M  «  0.5 
h  *  40K  ft 

4- I MAG I NARY  AXIS 
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’  v'  »■ .  ■ .  ■ . .  J.  r»  ^rrrryrr  T.  ;■ 


.1.  f.rrr 


GAIN  OF  INTEREST  <GAIN>  -  -1.200 

SENSITIVITY  OF  INTEREST  (OLK)  -  -1.460 

ROOTS  OF  INTEREST 


♦OR-  GTOL  - 
♦OR-  KTOL 


X  ■ 

.54698E-02 

Y  - 

-.51732E 

X  « 

.54698E-02 

Y  - 

.51732E 

X  - 

-. 1 1778E-0  1 

Y  - 

-4 . 7494 

X  - 

-. 11778E-0I 

Y  - 

4.7494 

REGION  OF 

CALCULATION-REAL! 

CC« 

-5.00 

IMAJ I 

DO- 

-3.00 

TO  AA- 
TO  BB- 


BRANCH  STARTING  AT  (-.013728)  ♦  JC . 094651) 
42 


ENTER  GAIN  OF  INTEREST  < GAIN) .TOLERANCE  (GT0L)s>-I.2 


1.00 

3.00 


1000 
.  1217 


OPEN-LOOP  < OLTF)  ROOT  LOCUS  USING  OPTION  42 
CODE 

0-LOCUS  PT.  2-ZERO 

1-POLE  3-BREAK  PT. 


M  *  0.6 

h  *  o  ft 

4-  I MAG I NARY  AXIS 

5- SENSITIVITY  PT. 


GAIN  OF  INTEREST  (GAIN)  -  -1.200  +0R-  GTOL  -  .1000 

SENSITIVITY  OF  INTEREST  (OLK)  *  -13.27  +08-  KTOL  -  1.106 

ROOTS  OF  INTEREST 


X  « 

-.98176E-02 

Y  « 

-.91793E-02 

X  » 

-.98176E-02 

Y  » 

.91793E-02 

X  « 

-.82973 

Y  - 

6.7733 

X  - 

-.82973 

Y  ■ 

-6.7733 

REGION  OF 

• 

CALCULATION-REAL! 

CO 

-5.00  TO 

AA» 

1.00 

IMAJ : 

DD» 

-3.00  TO 

BB* 

3.00 

BRANCH  STARTING  AT  (-.010817)  ♦  J( .06691) 

42 

ENTER  GAIN  OF  INTEREST  (GAIN) .TOLERANCE  (GTOL) :>-l .2, . 1 


OPEN-LOOP  (OLTF)  ROOT  LOCUS  USING  OPTION  42 
CODE 

0-LOCUS  PT.  2-ZERO 

1-POLE  3-BREAK  PT. 


M  =  0.6 

h  ='  20K  ft 

4-  I MAG I NARY  AXIS 

5- SENSITIVITY  PT. 


GAIN  OF  INTEREST  (GAIN)  *  -1.200 


♦  OR-  GTOL  «  .  1000 


SENSITIVITY  OF  INTEREST  (OLK)  -  -5.409  +  0R-  KTOL  *  .4507 


ROOTS  OF  INTEREST 
X  *  19634E-02 

X  ■  - .  19634E-02 

X  •  -.37512 

X  ■  -.37312 


Y  -  - . 24 173E-0 1 

Y  -  . 24  173E-0  1 

Y  -  -4.9676 

Y  ■  4.9676 
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TO  BB- 


3.00 


IMAJj  DO-  -3.00 
BRANCH  STARTING  AT  (-.0070548)  ♦  J(  . 076164) 

42 

ENTER  GAIN  OF  INTEREST  (GAIN) .TOLERANCE  (GTOL) :>-l .2, . 1 


OPEN-LOOP  ( OLTF)  ROOT  LOCUS  USING  OPTION  42  M  «  0.6 

h  *  40K 

COOE 

0 -LOCUS  PT.  2-2 ERO  4- I MAG I NARY  AXIS 

1-POLE  3-BREAK  PT.  3-SENSITIVITY  PT . 

GAIN  OF  INTEREST  (GAIN)  -  -1.200  ♦OR-  GTOL  -  .1000 

SENSITIVITY  OF  INTEREST  (OLK)  -  -2.12?  +OR-  KTOL  »  .1774 


ROOTS  OF  INTEREST 

X  •  .40 106E-02  Y  -  -.38340E-01 

X  -  .40 106E-02  Y  -  . 38340E-0 1 


X  -. 

-.11849 

Y  - 

4.2240 

X  - 

-.11849 

Y  - 

-4.2240 

REGION  OF  CALCULATION-REAL  I  CC-  -5.00  TO  AA-  1.00 

IMAJj  DO-  -3.00  TO  BB*  3.00 

BRANCH  STARTING  AT  (-.010247)  ♦  J( .080888) 

4 2 

ENTER  GAIN  OF  INTEREST  (GAIN) .TOLERANCE  (GTOL) i>-1 .2, . 1 


OPEN-LOOP  (OLTF)  ROOT  LOCUS  USING  OPTION  42 


M  *  0.8 
h  *  0  ft 


COOE 

0-LOCUS  PT.  2-2 ERO  4- I MAG I NARY  AXIS 

1-POLE  3-BREAK  PT.  5-SENSITIL'ITY  PT. 

GAIN  OF  INTEREST  (GAIN)  -  -1.200  +OR-  GTOL  -  .1000 

SENSITIVITY  OF  INTEREST  (OLK)  -  -25.52  ♦OR-  KTOL  -  2.127 

ROOTS  OF  INTEREST 


X 

- 

-.33280E-0  1 

Y  - 

0. 

X 

m 

-.64486E-03 

Y  - 

0. 

X 

- 

-1.2305 

Y  » 

9.2053 

X 

- 

-1.2305 

Y  - 

-9.2053 

REGION  OF  CALCULATION-REAL j  CC-  -5.00  TO  AA-  1.00 

IMAJ :  DO-  -3.00  TO  BB-  3.00 

BRANCH  STARTING  AT  (-.016487)  ♦  J(. 04 1362) 

42 


ENTER  GAIN  OF  INTEREST  (GAIN) .TOLERANCE  (GTOL) »>-l .2. . 1 


OPEN-LOOP  (OLTF)  ROOT  LOCUS  USING  OPTION  42 


M  =  1.0 
h  =  20K  ft 


CODE 

0-LOCUS  PT.  2-ZERO  4-JMAGINARY  AXIS 

1-POLE  3-BREAK  PT.  5-SENSITIVITY  PT . 
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Wv'var 


.  • . 


GAIN  OF  INTEREST  (GAIN)  -  -1.280 

SENSITIVITY  OF  INTEREST  (OLK)  »  -16.6? 


ROOTS  OF  INTEREST 
X  -  -.24510E-01 

X  -  . I8225E-02 

X  -  -.56668 

X  -  -.56668 


Y  ■ 

Y  - 

Y  - 

Y  - 


0. 

0. 

-6.5666 

6.5666 


♦OR-  GTOL  - 
♦OR-  KTOL 


.  1000 
11391 


Lf. 

t'. 


i 

i 


/■V-- '  AV 


REGION  OF  CALCULATI ON-REAL i  CC-  -5.00 

IMAJ :  DD=  -3.00 

BR«4CH  STARTING  AT  (-.03975)  ♦  J(0.) 

42 


TO  AA— 
TO  BB= 


1.00 

3.00 


ENTER  GAIN  OF  INTEREST  (GAIN) .TOLERANCE  (GTOL) :>-l .2, . 1 


OPEN-LOOP  ( OLTF)  ROOT  LOCUS  USING  OPTION  42 


M  =  1.0 

h  =  40K  ft 


CODE 

0 -LOCUS  PT. 
1-POLE 


2- 2 ERO 

3- BREAK  PT. 


GAIN  OF  INTEREST  (GAIN)  -  -1.208 

SENSITIVITY  OF  INTEREST  .  (OLK)  -  -6.167 

ROOTS  OF  INTEREST 


4-  IMAGINARY  AXIS 

5- SENSITIVITY  PT. 

♦OR-  GTOL  *  .1080 

♦OR-  KTOL  «  .5139 


X  » 

-.37698E-02 

Y  » 

-. 12623E 

X  » 

-.37698E-02 

Y  - 

. 12623E 

X  - 

-.23649 

Y  - 

-4.7865 

X  « 

-.23649 

Y  « 

4.7065 

ON  OF 

CALCULATI  ON-REAL : 

CC- 

-5.00 

IMAJ: 

DD* 

-3.00 

TO  AA- 
TO  BB=* 


1.00 

3.00 


BRANCH  STARTING  AT  (-.0067527)  ♦  J( .047972) 

42 


ENTER  GAIN  OF  INTEREST  ( GAIN) .TOLERANCE  (GTOL) :>-l .2, . 1 

OPEN-LOOP  (OLTF)  ROOT  LOCUS  USING  OPTION  42 
CODE 


M  =  1.2 
h  =  20K  ft 


0-LOCUS  PT.  2- 

ZERO 

4- IMAGINARY  AXIS 

1-POLE  3- 

BREAK 

PT. 

5-SENSITIVITY  PT. 

GAIN  OF  INTEREST  (GAIN)  - 

-1. 

200 

♦  OR-  GTOL  =  .  1000 

SENSITIVITY  OF  INTEREST  (OLK)  - 

-24.93 

♦OR-  KTOL  -  2.078 

ROOTS  OF  INTEREST 

X  -  - . 34739E-0  1 

Y  - 

0. 

X  -  .32160E-02 

Y  - 

8. 

X  -  -.56932 

Y  - 

6.6404 

X  -  -.56932 

Y  - 

-6.6404 

REGION  OF  CALCULATION-REAL: 

CC- 

-5.00 

TO  AA-  1 .00 

IMAJ: 

DD- 

-3.00 

TO  BB-  3.00 
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42 


c® 


ENTER  GAIN  OF  INTEREST  (GAIN) .TOLERANCE  (GTOL) :>-i .2, . I 


OPEN-LOOP  (OLTF)  ROOT  LOCUS  USING  OPTION  42 


M  *  1.4 
h  *  40K  ft 


CODE 

•-LOCUS  PT. 

1-POLE 


2- 2ERO 

3- BREAK  PT. 


4- IMAGINARY  AXIS 

5—  SENS I T I V I TY  PT. 


GAIN  OF  INTEREST  (GAIN)  »  -1.200 

SENSITIVITY  OF  INTEREST  <OLK>  - 


ROOTS  OF  INTEREST 
X-  — . 15595E-0 1 
X  «  -.36254E-03 

X  -  —.26432 

X  -  -.26632 


-12.28 


0. 

0. 

4.9925 

-4.9925 


♦OR-  GTOL 


♦OR-  KTOL 


.  1000 


1.024 


REGION  OF  CALCULATION-REAL:  CC-  -5.00 

IMAJ:  DO-  -3.00 
BRANCH  STARTING  AT  (-.0082773)  ♦  JC . 039491) 


TO  AA«  1.00 
TO  BB»  3.00 


I  c 


GAIN— 1 ,43 

ENTER  ZETA, RAD, GTOL  >.7,1,.  I 


OPEN-LOOP  ( OLTP)  ROOT  LOCUS  USING  OPTION  43 
CODE 

0-LOCUS  PT.  2-ZERO 

1-POLE  3-BREAK  PT. 

DAMPING  FACTOR  OF  INTEREST  ( ZETA)  ■  .7000  . 

GAIN  OF  INTEREST  <GAIN)  -  3.632 

SENSITIVITY  OF  INTEREST  (OLK)  -  -3.713 


M  -  0.3 
h  *  o  ft 


4-  IMAGINARY  AXIS 

5- SENSITIVITY  PT. 

♦OR-  GTOL  =  .365 

♦OR-  KTOL  =  .5 


ROOTS  OF  INTEREST 
X  -  .263 15E-02 

X  «  .26313E-02 

X  -  -3.1471 

X  -  -3.1471 


Y  -  . 54473E— 0 1 

Y  ■  • -.54473E-0 1 

Y  «  3.2105 

Y  -  -3.2105 


REGION  OF  CALCULATION-REAL  I  CC«  -10.0  TO  AA»  1.60 

IMAJj  DD-  -10.0  TO  BB»  10.0 

BRANCH  STARTING  AT  (.0043183)  ♦  J< .060757) 

GA 114— 1,43 

ENTER  ZETA, RAD, GTOL  >.7,l,.l 


OPEN-LOOP  < OLTF)  ROOT  LOCUS  USING  OPTION  43 


M  *  0.4 
h  =  20K  ft 


CODE 

0-LOCUS  PT.  2-ZERO  4- I MAG I NARY  AXIS 

1-POLE  3-BREAK  PT.  3-SENSITIVIT/  PT. 


DAMPING  FACTOR  OF  INTEREST  (ZETA)  -  .7000 

GAIN  OF  INTEREST  (GAIN)  »  3.179  +0R-  GTOL  *  .317' 


SENSITIVITY  OF  INTEREST  (OLK)  -  -5.799  +OR-  KTOL  *  .5: 


ROOTS  OF  INTEREST 
X  -  . 37230E-02 

X  -  .37230E-02 

X  «  -3.0781 

X  «  -3.0781 


REGION  OF  CALCULATION-REAL: 

IMAJ : 

BRANCH  STARTING  AT  (.0046977)  ♦ 
GAIN— 1 ,43 


Y  =  . 48477E-0 1 

Y  =  - . 48477E-0 I 

Y  -  -3.  140  1 

Y  ■  3. 1401 


CC-  -10.0  TO  AA»  1.0O 
DD*  -10.0  TO  BB*  10.0 
J( .052395) 


ENTER  ZETA,  RAD,  GTOL  >.7,1,'.  1 


OPEN-LOOP  (OLTF)  ROOT  LOCUS  USING  OPTION  43 


CODE 

0-LOCUS  PT.  2-ZERO  4-IMAGINAPY  A.cJS 

1-POLE  3-BREAK  PT.  •  5-SENSITIVITY  PT . 

DAMPING  FACTOR  OF  INTEREST  (2ETA)  -  .7000 

GAIN  OF  INTEREST  (GAIN)  «  2.079  +  OR-  GTOL  -  .2079 

SENSITIVITY  OF  INTEREST  (OLK)  -  -6.925  *OR-  KTOL  «  .6925 

ROOTS  OF  INTEREST  M  =  0.5 

X  «  .5 1096E-02  Y  »  . 4'954 IE-0 1  h  *  40K  ft 

X  ■  .51096E-02  Y  »  - . 4954 1 E-0  1 

X  -  -3.4739  Y  -  -3.5441 

X  ■  -3.4739  Y  -  3.5441 


REGION  OF  CALCULATION-REAL:  CC=  -10.0  TO  AA»  1.00 

IMAJ :  DO*  -10.8  TO  BB=*  10.0 

BRANCH  STARTING  AT  (.0054693)  +  J(. 051752) 

GAIN— 1,43 

ENTER  2 ETA , RAD , GTOL  >.7,l,.l 


OPEN-LOOP  (OLTF)  ROOT  LOCUS  USING  OPTION  43 
CODE 

0-LOCUS  PT.  2-ZERO 

1-POLE  3-BREAK  PT. 

DAMPING  FACTOR  OF  INTEREST  (ZETA)  -  .7000 

GAIN  OF  INTEREST  (GAIN)  -  8.545  +  OR-  GTOL  *  .S545 

SENSITIVITY  OF  INTEREST  (OLK)  •  -9.343  +OR-  KTOL  =  .9343 

ROOTS  OF  INTEREST 

X«  - . 10 139E-0 1  Y-  . 57812E-02 

X  ■  —.10 139E-0 1  Y  -  - . 578 12E-02 

X  »  -5.5008  Y  **  -5.6117 

X  ■  -5.5008  Y  -  5.6117 


REGION  OF  CALCULATION-REAL:  CC*  -18.0  TO  AA»  1.00 

IMAJ:  DD»  -10.0  TO  BB=  .  10.0 

BRANCH  STARTING  AT  (-.0098176)  +  J< . 0091793) 

GAIN—  1 ,43 

ENTER  ZETA, RAD, GTOL  >.7,i,.i 


M  «  0.6 
h  e  0  ft 

4- IMAGINARY  AXIS 

5- SENSITIVITY  PT. 


OPEN-LOOP  (OLTF)  ROOT  LOCUS  USING  OPTION  43 
CODE 

0-LOCUS  PT.  2-ZERO 

1-POLE  3-BREAK  PT . 


M  =  0.6 
h  =  20K  ft 

4- IMAGINARY  «<IS 

5- SENSITIVITY  PT . 


DAMPING  FACTOR  OF  INTEREST  (ZETA)  =  .7000 

GAIN  OF  INTEREST  (GAIN)  »  5.254  ♦  OR-  GTOL  *  .5254 


SENSITIVITY  OF  INTEREST  (OLK)  - 


-6.918 


♦OR-  KTOL 


.6918 


ROOTS  OF  INTEREST 
X  *  -.24673E-02 

X  *  -.24673E-02 

X  *  -3.8334 

X  *  -3.8334 


REGION  OF  CALCULATION-REAL; 

IMAJi 

BRANCH  STARTING  AT  (-.0019634) 
GAIN-- 1 ,43 


Y  -  . 21923E-0 1 

Y  -  -.21923E-0  1 

Y  *  3.9110 

Y  -  -3.9110 


CC-  -10.0  TO  AA—  1.00 

DO-  -10.0  TO  88*  10.0 

♦  J(. 024173) 


ENTER  ZETA.RAD.GTOL  >.7,l,.l 


OPEN-LOOP  ( OLTF)  ROOT  LOCUS  USING  OPTION  43 
CODE 

0-LOCUS  PT.  2-ZERO 

1-POLE  3-BREAK  PT. 

DAMPING  FACTOR  OF  INTEREST  (ZETA)  «  .7000 

GAIN  OF  INTEREST  (GAIN)  -  3.003  +OR-  6T0L  *  .3003 


SENSITIVITY  OF  INTEREST  (OLK)  - 

-4.011  fOR- 

KTOL 

• 

ROOTS 

OF  INTEREST 

X  * 

.35794E-02 

Y  * 

. 3433 1E-0 1 

X  » 

•35796E-02 

Y  * 

-.34331E-0 1 

X  * 

-3. 1234 

Y  * 

3.  1845 

X  * 

-3. 1234 

Y  * 

-3. 1845 

REGION  OF  CALCULATION-REAL; 

CC* 

-10.0  TO  AA- 

1 . 

00 

IMAJi 

DD« 

-10.0  TO  BB» 

10 

.0 

BRANCH  STARTING  AT  (*.0040  104)  ♦ 

J( .03834) 

GAIN— 1 ,43 

ENTER  2 ETA , RAD , GTOL  >.7,l,.l 


M  -  0.6 
h  *  40K  ft 


4-  I MAG I NARY  AXIS 

5- SENSITIVITY  PT. 


OPEN-LOOP  (OLTF)  ROOT  LOCUS  USING  OPTION  43 


M  =  0.8 
h  =  0  ft 


CODE 

0-LOCUS  PT.  2-ZERO  4- I MAG I NARY  AXIS 

l-POLE  3-BREAK  PT.  5-SENSITIVITY  PT. 


DAMPING  FACTOR  OF  INTEREST  (ZETA)  -  .7000 


GAIN 

OF  INTEREST  (GAIN) 

as 

12 

.21 

♦OR-  GTOL  * 

SENSITIVITY  OF  INTEREST 

(OLK) 

- 

-12.49 

♦OR-  KTOL  - 

ROOTS 

OF  INTEREST 

X  * 

- . 33683E-0 1 

Y 

■ 

0. 

X  * 

-.46?48£-63 

Y 

m 

0. 

X  » 

-7.5734 

Y 

m 

7.7264 

X  ■ 

-7.5734 

Y 

m 

-7.7264 

REGION  OF  CALCULATION-REAL; 

CC* 

-15.0 

TO 

AA* 

1.00 

IMAJs 

00* 

-15.0 

TO 

BB* 

15.0 

1.221 

1.269 
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gain*-  i  , 43 

ENTER  ZETA,RAD,GTOL  >.7,l,.l 


OPEN-LOOP  <OLTF>  ROOT  LOCUS  USING  OPTION  43 
CODE 

0-LOCUS  PT.  2-ZERO 

1-POLE  3-BREAK  PT. 


M  -  1.0 
h  «  20K  ft 

4-  I HAG I NARY  AXIS 

5- SENSITIVITt  PT. 


DAMPING  FACTOR  OF  INTEREST  (ZETA)  =  .7000 


GAIN  OF  INTEREST  (GAIN)  -  10.08  +0R-  GTOL  -  1.008 

SENSITIVITY  OF  INTEREST  (OLK)  -  -9.190  +  0R-  KTOL  *  .9190 

ROOTS  OF  INTEREST 

X  *  -.24328E-0 1  Y  «  0. 

X  *  . 1467 IE-02  Y  -  0. 

X*  -5.1616  Y«  -5.2656 

X  «  -5. 1616  Y  -  5.2656 


REGION  OF  CALCULATION-REAL!  CC»  -10.0  TO  AA*  1.00 

IMAJi  DD“  -10.0  TO  BB=*  10.0 

BRANCH  STARTING  AT  (-.02451)  «•  J(0.) 

GAIN— 1,43  . 

ENTER  ZETA,RAD,GTOL  >.7,1,. 1 


OPEN-LOOP  ( OLTF)  ROOT  LOCUS  USING  OPTION  43 
CODE 

0 -LOCUS  PT.  2-ZERO 

1-POLE  3-BREAK  PT. 

DAMPING  FACTOR  OF  INTEREST  (ZETA)  «  .7000 

GAIN  OF  INTEREST  (GAIN)  -  13.45  +OR-  GTOL  -  1.345 

SENSITIVITY  OF  INTEREST  (OLK)  *  -9.438  +OR-  KTOL  *  .9438 


M  *  1.2 
h  *  20K  ft 

4- I MAG I NARY  AXIS 
"5-SENSITIVITY  PT. 


ROOTS  OF  INTEREST 
X  -  -.33971E-01 

X  ■  . 25598E-02 

X  *  -5.2882 

X  «  -5.2882 


Y  «  8. 

Y  *  0 . 

Y  -  5.3946 

Y  -  -5.3946 


REGION  OF  CALCULATION-REAL; 

cc- 

-15. 

0 

TO 

AA* 

1  .00 

IMAJi 

DD* 

-15. 

0 

TO 

BB* 

15.0 

BRANCH  STARTING  AT  (-.034739)  ♦  J<0.) 
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GAIN*- l ,43 

ENTER  ZETA,RAD,GTOL  >.7,l,.l 


OPEN-LOOP  <OLTF>  ROOT  LOCUS  US I NO  OPTION  43 
CODE 

0-LOCUS  FT.  2-ZERO 

1-POLE  3-BREAK  PT. 


M  «  1.4 
h  *  40K  ft 

4-IMAGINARY  AXIS 
3-SENSITIVITY  PT. 


OAMP1NG  FACTOR  OF  INTEREST  ( ZETA)  ■  .7000 

GAIN  OF  INTEREST  (GAIN)  -  10.23  *OR-  GTOL  -  1.023 

SENSITIVITY  OF  INTEREST  (OLK)  -  -7.136  +  OR-  KTOL  *  .7136 


ROOTS  OF  INTEREST 
X  ■  — .  29750E-03 

X-  -. 15833E-0 1 
X  -  -3.8340 

X  -  -3.8340 


¥ 

Y 

Y 

Y 


0. 

0. 

3.9113 

-3.9113 


REGION  OF  CALCULATION-REAL:  CC-  -10.0 

IMAJt  DO-  -10.0 

BRANCH  STARTING  AT  (-.015595)  ♦  J(0.) 


TO 
TO  BB» 


1.00 

10.0 


N 


A 

1 
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42 


Method  I 


ENTER  CAIN  OF  INTEREST  < GAIN)  , TOLERANCE  ( GTOL)  : >3 . 6 1 , .  1 


OPEN-LOOP  < OLTF)  ROOT  LOCUS  USING  OPTION  42  M  ■  0.3 

h  *  0  ft 

CODE 

0-LOCUS  PT.  2-ZERO  4— IMAGINARY  AXIS 

1-POLE  3-BREAK  PT.  5-SENSITIVITY  PT. 

GAIN  OF  INTEREST  (GAIN)  -  3.610  +  OR-  6T0L  -  .1000 


SENSITIVITY  OF  INTEREST  <OLK>  -  -5.647  ^OR-  KTOL  «  .1564 


ROOTS  OF  INTEREST 
X  «  . 264S2E-02 

X  »  . 26482E-02 

X  -  -3.113? 

X  «  -3.1139 


Y  -  . 54535E-0 1 

Y  ■  -.54535E-0 1 

Y  •  -3.2356 

Y  -  3.2356 


REGION  OF  CALCULATION-REAL:  CD*.  -5.00  TO  AA»  1.00 

IMAJ:  DD-  -3.00  TO  3B=  3.00 

BRANCH  STARTING  AT  (.0043183)  ♦  J( .060757) 

42 

ENTER  GAIN  OF  INTEREST  (GAIN) .TOLERANCE  (6T0L) :>-3.33, . 1 


OPEN-LOOP  (OLTF)  ROOT  LOCUS  USING  OPTION  42 
COOE 

0-LOCUS  PT.  2-ZERO 

1-POLE  3-BREAK  PT. 


M  =  0.4 
h  «  20K  ft 

4—  IMAGINARY  AXIS 

5—  SENS I T I V I TV  PT. 


GAIN  OF  INTEREST  (GAIN)  «  -3.330  ♦OR-  GTOL  -  .1000 


SENSITIVITY  OF  INTEREST  (OLK)  ■  6.074  *0R-  KTOL  -  .1824 


ROOTS  OF  INTEREST 
X  »  . 368 14E-02 

X  ■  . 366 14E-02 

X  -  -3.2155 

X  -  -3.2155 


Y  -  . 483 13E-0 1 

Y  -  - . 483 1 3E-0 1 

Y  -  -3.0213 

Y  *  3.0213 


REGION  OF  CALCULATION-REAL:  CC-  -5.00  TO  AA=  1.00 

IMAJ:  DD=  -3.00  TO  BB=  3.00 

BRANCH  STARTING  AT  (.0046977)  ♦  J( . 052395) 

42 

ENTER  GAIN  OF  INTEREST  ( GAIN) .TOLERANCE  ( GTOL) : >-2 . 87 , . 1 


OPEN-LOOP  (OLTF)  ROOT  LOCUS  USING  OPTION  42 
CODE 

0-LOCUS  PT.  2-ZERO 

1-POLE  3-BREAK  PT. 


M  =  0.5 
h  =  40K  ft 

4-  IMAGINARY  AXIS 

5- SENSITIVITY  PT. 


GAIN  OF  INTEREST  (GAIN)  ■  -2.870  ♦OR-  GTOL  ■ 

SENSITIVITY  OF  INTEREST  (OLK)  •  9.561  ♦OR-  KTOL  » 


ROOTS  OF  INTEREST 
X  -  .49740E-02 

X  -  .49748E-02 

X  -  -4.7920 

X  -  -4.7920 


REGION  OF  CALCULAT I  ON-REAL t 

IMAJs 

BRANCH  STARTING  AT  (.0054698)  ♦ 


Y  »  .48774E-0 1 

Y-  - . 48774E-0  1 

Y  «  1.5657 

Y  »  -1.5657 


CO  -5.00  TO  AA-  1.00 
DO  -3.00  TO  BO  3.00 
J(. 051752) 


.  1000 
.3332 
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42 

ENTER  GAIN  OF  INTEREST  (GAIN)  , TOLERANCE  (GTOL) :>-8. 19, .  1 


OPEN-LOOP  (OLTF)  ROOT  LOCUS  USING  OPTION  42 
CODE 

0-LOCUS  PT.  2-ZERO 

I -POLE  3-BREAK  PT. 

GAIN  OF  INTEREST  (GAIN)  •  -8.190 

SENSITIVITY  OF  INTEREST  (OLK)  -  8.954 


M  *  0.6 
h  s  0  ft 

4-  IMAGINARY  AXIS 

5- SENSITIVITY  PT. 

♦OR-  GTOL  ■  . 1000 

♦OR-  KTOL  »  . 1093 


ROOTS  OF  INTEREST 
X«  -. 10 129E-0 1 
X*  -.101 29E-01 
X  «  -5.3044 

X  «  -5.3044 


Y  -  .59 194E-02 

Y  -  -. 59194E-02 

Y  ■  -5.7410 

Y  «  5.7410 


REGION  OF  CALCULATION-REAL!  CC«  -5.00  TO  AA*=  1.00 

IhAJs  DD»  -3 .‘00  TO  BB*  3.00 

BRANCH  STARTING  AT  (-.0098174)  ♦  J( .0891793) 

42 

ENTER  GAIN  OF  INTEREST  (GAIN) .TOLERANCE  ( GTOL)  : >-4 . 89 , .  1 


OPEN-LOOP  ( OLTF)  ROOT  LOCUS  USING  OPTION  42 


M  *  0.6 

h  «  20K  ft 


CODE 

0-LOCUS  PT.  2-ZERO  4- IMAGINARY  AXIS 

1-POLE  3-BREAK  PT.  5-SENSITIVIT .  PT. 


GAIN  OF  INTEREST  (GAIN)  »  -4.890  +OR-  GTOL  -  .1000 


SENSITIVITY  OF  INTEREST  (OLK)  -  4.439  ♦  OR-  KTOL  •  .1317 


ROOTS  OF  INTEREST 


X 

m 

-.24374E-02 

Y  ■ 

.22040E-0  1 

X 

m 

-.24374E-02 

Y  - 

-.22040E-0  1 

X 

m 

-3.5940 

Y  » 

4.0887 

X 

m 

-3.5940 

Y  - 

-4.0887 

REGION  OF  CALCULATION-REAL: 

CC- 

-5.00 

TO 

AA* 

ITWJi 

DD- 

-3.00 

TO 

BB« 

BRANCH  STARTING  AT  (-.0019434)  ♦  J( . 024173) 

42 

ENTER  GAIN  OF  INTEREST  (GAIN)  .TOLERANCE  (GTOL)  :>-3.21 , .  1 


OPEN-LOOP  (OLTF)  ROOT  LOCUS  USING  OPTION  42 


CODE 

0-LOCUS  PT. 
1-POLE 


2- ZERO 

3- BREAK  PT. 
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M  =  0.6 
h  =  4 OK  ft 

4-  I MAG I MARY  A* I S 

5- SENSIT1 VI  T'i  r.T 


GAIN  OF  INTEREST  (GAIN) 

-  -3.210 

♦OR-  GTOL  - 

.  1000 

SENSITIVITY  OF  INTEREST 

(OLK)  «  4.425 

♦OR-  KTOL  - 

.2001 

ROOTS  OF  INTEREST 

X  -  .35518E-02  Y  -  .36205E-01 

X  ■  . 35518E-02  Y  -  -.36205E-01 

X  -  -3.3303  Y  -  -2.9933 

X  -  -3.3303  Y  «  2.9933 


REG! ON  OF  CALCULATION-REAL!  CC-  -5.00  TO  AA-  1.00 

IMAJ!  DO-  -3.00  TO  BB-  3.00 

BRANCH  STARTING  AT  (.0040106)  ♦  J( .03834) 

42 

ENTER  GAIN  OF  INTEREST  (GAIN)  .TOLERANCE  (GTOL)  .*>-12.94, .  1 


OPEN-LOOP  ( OLTF)  ROOT  LOCUS  USING  OPTION  42  *  M  «  0.8 

h  e  0  ft 

CODE 

0-LOCUS  PT.  2-ZERO  4- I MAG I NARY  AXIS 

1-POLE  3-BREAK  PT.  5-SENSITIVITY  PT. 

GAIN  OF  INTEREST  (GAIN)  -  -12.94  +  OR-  GTOL  -  .1000 

SENSITIVITY  OF  INTEREST  (OLK)  «  13.45  +  0R-  KTOL  -  .1039 

ROOTS  OF  INTEREST 

X-  -.33700E-0 1  Y-  0. 

X  -  -.46191E-03  •  Y  -  0. 

X  -  -7.9355  Y  -  7.4579 

X  -  -7.9555  Y  «  -7.4579 


REGION  OF  CALCULATION-REAL!  CC-  -5.00  TO  AA-  1.06 

IMAJ :  00-  -3.00  TO  BB-  3.00 

BRANCH  STARTING  AT  (-.03328)  ♦  J(0.) 

42 

ENTER  GAIN  OF  INTEREST  (GAIN) .TOLERANCE  (GTOL) :>-9.87, . 1 


OPEN-LOOP  (OLTF)  ROOT  LOCUS  USING  OPTION  42 
CODE 

0-LOCUS  PT.  2-ZERO 

1-POLE  3-BREAK  PT. 

GAIN  OF  INTEREST  (GAIN)  -  -9.870 

SENSITIVITY  OF  INTEREST  (OLK)  -  8.994 


M  *  1 

h  *=  20K  ft 

4-  I MAG I NARY  AXIS 

5- SENSITIVITY  PT. 

♦OR-  GTOL  -  . 1000 

♦OR-  KTOL  -  .9113E-0 


ROOTS  OF  INTEREST 

X  -  -.24331E-01  Y  -  0. 

X  -  . 14732E-02  Y  -  0. 

X  -  -5.0438  Y  -  -5.3380 

X  -  -5.0438  Y  -  5.3380 


REGION  OF  CALCULATION-REAL:  PC' 


-5 .  AO 


TO  AA' 


l  .OO 


IMAJ:  DO-  -3.00 

BRANCH  STARTING  AT  (-.02451)  ♦  J(0.)  • 

42 


TO  B6»  3.00 


ENTER  GAIN  OF  INTEREST  (GAIN) .TOLERANCE  ( GTOL) : >-5 . 22 , . I 


OPEN-LOOP  (OLTF)  ROOT  LOCUS' USING  OPTION  42 


CODE 

•-LOCUS  PT. 
1-POLE 


2- ZERO 

3- BREAK  PT. 


M  »  1  .0 

h  «  4 OK  ft 

4-IMAGINARY  AXIS 
3-SENSITIVITY  PT. 


GAIN  OF  INTEREST  (GAIN)  -  -3.220 


♦OR-  GTOL  ■  . 1000 


SENSITIVITY  OF  INTEREST  (OLK) 


5.774  ♦OR-  KTOL  -  .1106 


ROOTS  OF  INTEREST 
X  «  -.39191E-02 

X  »  -.39191E-02 

X  «  -3.1235 

X  «  -3. 1235 


Y  -  .  1 1737E-0 l 

Y«  - . 1 1737E-0 1 

Y  »  -3.9261 

Y  -  3.9261 


REGION  OF  CALCULATION-REAL:  CC«  -5.00 

IMAJ:  DD-  -3.00 
BRANCH  STARTING  AT  (-.0037693)  ♦  J(. 012623) 
42 


TO  AA-  1.00 
TO  BB*  3.08 


ENTER  GAIN  OF  INTEREST  (GAIN) .TOLERANCE  ( GTOL) : >- 13 . 3 , . 1 


OPEN-LOOP  (OLTF)  ROOT  LOCUS  USING  OPTION  42 


CODE 

• -LOCUS  PT. 
1-POLE 


2- 2ER0 

3- BREAK  PT. 


M  *  1.2 

h  *  20K  ft 

4- IMAGINARY  AXIS 

5- SENSITIVITY  PT. 


GAIN  OF  INTEREST  (GAIN) 


-13.30 


♦OR-  GTOL 


SENSITIVITY  OF  INTEREST  (OLK) 


9.334 


♦OR-  KTOL 


. 70 18E- 


ROOTS  OF  INTEREST 
X  »  -.33977E-0 1 

X  -  . 25636E-02 

X  «  -5.2366 

X  -  -5.2366 


Y  -  0. 

Y  -  0. 

Y  -  5.4320 

Y  -  -3.4320 


REGION  OF  CALCULATION-REAL:  CC-  -3.00 

IMAJ:  DO-  -3.00 

BRANCH  STARTING  AT  (-.034739)  ♦  J(0.) 

42 


TO  AA*  1.00 
TO  BB»  3.00 


ENTER  GAIN  OF  INTEREST  (GAIN) .TOLERANCE  (GTOL) :>-8.24, . 1 


OPEN-LOOP  (OLTF)  ROOT  LOCUS  USING  OPTION  42 


COOE 

0-LOCUS  PT. 
1-POLE 


2- ZERO 

3- BP‘AK  PT. 


M  «  1.4 
h  ■  40K  ft 

4- IMAGINARY  AXIS 

5- SENSITIVITY  PT . 


GAIN  OF  INTEREST  <GAIN> 


8.240 


SENSITIVITY  OF  INTEREST  <OLK>  ■  5.734 

1 


ROOTS 

OF  INTEREST 

X  - 

-.30835E-03 

Y  ■ 

0. 

X  - 

15795E-0 1 

Y  ■ 

0. 

X  - 

-3.1342 

Y  - 

-4.3811 

X  « 

-3. 1342 

Y  - 

4.3811 

REGION  1 

DF  CALCULATION-REALs 

CC- 

-5.08 

IMAJs 

OD" 

-3.08 

BRANCH  STARTING  AT  <-.015595)  ♦ 

j(0.; 

> 

% 


42 


Method  .n 


Gb 


ENTER  GAIN  OF  INTEREST  < GAIN) .TOLERANCE  (GTOL) j>-3.77, . 1 


OPEN-LOOP  ( OLTF)  ROOT  LOCUS  USING  OPTION  42 


CODE 

0 -LOCUS  PT. 

1-POLE 

GAIN  OF  INTEREST  (GAIN) 


2- ZERO 

3- BREAK  PT. 


M  »  0.3 
h  «  0  ft 

4- IMAGINARY  AXIS 
3-SENS IT IVITY  PT. 


-3.770 


SENSITIVITY  OF  INTEREST  (OLK) 


♦OR-  GTOL  « 
3.897  +0R-  KTOL 


.  1080 


.  1564 


ROOTS  OF  INTEREST 
X  «  .25857E-02 

X  -  . 25857E-02 

X  -  -3.2390 

X  «  -3.2390 


Y  ■  . 5430 2E-0 1 

Y  -  -.54302E-0 1 

Y  «  3. 1384 

Y  »  -3.1384 


REGION  OF  CALCULATION-REAL:  CC«  -5.00 

IMAJ:  DO-  -3-.  00 
BRANCH  STARTING  AT  (.0043183)  ♦  J(. 060757) 
42 


TO  AA-  1.00 
TO  BB»  3.00 


ENTER  GAIN  OF  INTEREST  (GAIN) .TOLERANCE  ( GTOL) : >-3 . 35 , . 1 


OPEN-LOOP  (OLTF)  ROOT  LOCUS  USING  OPTION  42 


CODE 

0-LOCUS  PT. 

1-POLE 

GAIN  OF  INTEREST  (GAIN) 


2- ZERO 

3- BREAK  PT. 


M  *  0.4 

h  *  20K  ft 

4- IMAGINARY  AXIS 

5- SENSITIVITY  PT. 


-3.350 


SENSITIVITY  OF  INTEREST  (OLK) 


♦OR-  GTOL 


6.  1 10  +0R-  KTOL 


.  1000 


.  1824 


ROOTS  OF  INTEREST 
X  -  .36759E-02 

X  -  .  36759E-02 

X  «  -3.2338 

X  -  -3.233S 


Y  -  . 48291E-0 1 

Y  -  - . 4829 1E-0 1 

Y  ■  3.0047 

Y  »  -3.0047 


REGION  OF  CALCULATION-REAL:  CC-  -5.00 

IMAJ:  DD»  -3.00 
BRANCH  STARTING  AT  (.0046977)  ♦  J( .052395) 
42 


TO  AA* 
TO  B6= 


ENTER  GAIN  OF  INTEREST  (GAIN)  .TOLERANCE  < GTOL)  : >-2 . 65 , .  1 


OPEN-LOOP  (OLTF)  ROOT  LOCUS  USING  OPTION  42 


0.5 

4 OK  ft 


CODE 

0 -LOCUS  PT. 
1-POLE 


2- ZERO 

3- BREAK  PT. 


4- IMAGINARY  AXIS 

5- SENSITIVITY  PT. 


•  .  *  _  »  «  _  *  .  *  >  m  ■  1  ’  ■  r  ■  'ji'ji  V-  m  ■  J-  *  .  •  «.  -  -T  • 


GAIN  OF  INTEREST  (GAIN)  ■  -2.630  ♦  OR-  GTOL  - 

SENSITIVITY  OF  INTEREST  (OLK)  -  8.82?  ♦  OR-  KTOL 


.  100 
.3 


ROOTS 
X  ■ 

X  - 
X  - 
X  « 


OF  INTEREST 
. 30 1 15E-62 
.501 13E-02 
-4.4236 
-4.4256 


Y  - 

Y  « 

Y  ■ 

Y  - 


.  48984E-0 1 
- .  48984E-0 1 
-2.368? 
2.3687 


REGION  OF  CALCULATION-REAL i  CO  -15.0 

IMAJ:  DO  -13.0 
BRANCH  STARTING  AT  (.0054698)  ♦  J(. 05 1732) 
42 


TO 

TO  BO 


1.0d 

15.0 


ENTER  GAIN  OF  INTEREST  (GAIN) .TOLERANCE  (GTOL) :>-9.42, . 1 


OPEN-LOOP  ( OLTF)  ROOT  LOCUS  USING  OPTION  42 


M  =  0.6 
h  *  0  ft 


CODE 

0 -LOCUS  PT.  2-ZERO 

1-POLE  3-BREAK  PT. 


4- I MAG I NARY  *XIS 
3-SENSITIVITY  PT. 


GAIN  OF  INTEREST  (GAIN)  «  -9.420 


♦  OR-  GTOL  ■  .  100i 


SENSITIVITY  OF  INTEREST  (OLK)  «  10.30  +0R-  KTOL  *  .11 

ROOTS  OF  INTEREST 
X-  - , 10 164E— 0  1 
X-  -.10164E— 01 
X  -  -5.9789 

X  -  -5.9789 


REGION  OF  CALCULATION-REAL:  CO  -15.0  TO  AA»  1.00 

IMAJ:  DO  -15.0  TO  BB=  15.0 

BRANCH  STARTING  AT  (-.0098176)  ♦  J(. 0091793) 

42 


Y  »  . 54389E-02 

Y  -  — .54389E-02 

Y  -  3.2494 

Y  «  -5.2494 


ENTER  GAIN  OF  INTEREST  (GAIN) .TOLERANCE  (GTOL) :>-5.58, . 1 


OPEN-LOOP  (OLTF)  ROOT  LOCUS  USING  OPTION  42 
CODE 

0-LOCUS  PT.  2-ZERO 

1-POLE  3-BREAK  PT. 


M  =  0.6 

h  *  20K  ft 

4- IMAGINARY  A* IS 

5- SENSITI  VIT'i  PT. 


GAIN  OF  INTEREST  (GAIN)  -  -5.380 


♦OR-  GTOL  =  .1001 


SENSITIVITY  OF  INTEREST  (OLK)  -  7.347  +  0R-  KTOL  *  .11 


ROOTS  OF  INTEREST 
X  -  - . 2493SE-0 2 

X  «  -.24935E-02 

X  -  -4.0482 

X  -  -4.0482 


Y  ■  . 2 1 802E-0  1 

Y  -  - . 2 1802E-0 1 

Y  -  3.7316 

Y  -  -3.7316 


REGION  OF  CALCULATION-REAL:  CC* 


15.0 


TO  AA= 


1 .00 


GAIN  OF  INTEREST  (GAIN)  ■  -10.43 


♦OR-  GTQL 


.  1000 


$  <£■ 


SENSITIVITY  OF  INTEREST  <OLK>  -  9.523 


ROOTS  OF  INTEREST 

X  *  -.24323E-0 1  Y  ■  0. 

X  -  . 14368E-02  Y  -  0. 

X  -  -5.3280  Y  -  -3. 1358 

X-  -3.3280  Y-  5.1338 


REGION  OF  CALCULATION-REAL:  CC=  -5.00 

IMAJ:  DO"  -3.00 

BRANCH  STARTING  AT  (-.02451)  ♦  J(0.) 

42 


♦OR-  KTOL 


.91 13E-0 


TO  AA"  1.00 
TO  BB"  3.00 


ENTER  GAIN  OF  INTEREST  (GAIN) .TOLERANCE  (GTOL) :>-5.92, . 1 


<a 


OPEN-LOOP  (OLTF)  ROOT  LOCUS  USING  OPTION  42 


CODE 

0-LOCUS  PT. 
1-POLE 


2- ZERO 

3- BREAK  PT. 


GAIN  OF  INTEREST  (GAIN)  -  -5.920 

SENSITIVITY  OF  INTEREST  (OLK)  » 


6.549 


M  =  1.0 
h  =  40K  ft 

4- IMAGINARY  axis 

5- SENSITIVITY  PT. 

♦OR-  GTOL  »  . 1000 

♦OR-  KTOL  *  . 1 106 


ROOTS  OF  INTEREST 
X  -  -.39368E-02 

X  «  -.39368E-02 

X  »  -3.5186 

X  »  -3.5166 


Y  -  .  1 1629E-0 1 

Y  »  - .  1 1629E-0 1 

.Y  -  -3.6391 

Y  »  3.6391 


REGION  OF  CALCULATI ON-REAL:  CC"  -5.00  TO  AA*  1.00 

IMAJ:  DO*  -3.00  TO  BB*  3.00 

BRANCH  STARTING  AT  (-.0037698)  ♦  J( . 012623) 

42 

ENTER  GAIN  OF  INTEREST  (GAIN) .TOLERANCE  (GTOL) :>-12.92, . I 


OPEN-LOOP  (OLTF)  ROOT  LOCUS  USING  OPTION  42 


CODE 

0-LOCUS  PT. 

1-POLE 

GAIN  OF  INTEREST  (GAIN) 


2- ZERO 

3- BREAK  PT. 


-12.92 


SENSITIVITY  OF  INTEREST  (OLK) 


9.068 


M  =  1.2 
H  *  20K  ft 

4-IMAGINf-P'r  AXIS 
3-SENS I T I V I TV  PT, 


♦OR-  GTOL 


.  1000 


♦  OR-  KTOL  =  .  70 18E-0 


ROOTS  OF  INTEREST 
X  -  -. 33995E-0 1 

X  -  . 23804E-02 

X  -  -5. 1033 

X  ■  -5. 1033 


Y  -  0. 

Y  ■  0. 


5.5252 

-5.5252 


jV, 


BRANCH  STARTING  AT  (-.0019634)  ♦  J( . 024173) 

42  • 

ENTER  GAIN  OF  INTEREST  (GAIN)  , TOLERANCE  (  GT0L>  :  >-3  ..18 ,  .  1 


OPEN-LOOP  ( OLTF)  ROOT  LOCUS  USING  OPTION  42 


CODE 

0-LOCUS  PT. 

1-POLE 


2- ZERO 

3- BREAK  PT. 


M  =  0.6 
h  =  40K  ft 


4-  I MAG I NARY  AXIS 

5- SENSITIVITY  PT. 


GAIN  OF  INTEREST  (GAIN) 


-3.180 


♦OR-  GTOL 


SENSITIVITY  OF  INTEREST  (OLK) 


6.365 


♦OR-  KTOL 


.2001 


ROOTS  OF  INTEREST 
X  -  .35553E-02 

X  «  .355582-02 

X  «  -3.3003 

X  «  -3.3003 


Y  *  . 36223E-0 1 

Y  -  - . 36223E-8 l 

Y  -  -3.8230 

Y  -  3.0236 


REG I  ON  OF  CALCULATION-REAL:  CC=  -15.8  TO  AA«  1.00 

IMAJj  DD=  -15. '8  TO  BB=  15.0 

BRANCH  STARTING  AT  (.0040106)  ♦  J( .03834) 

42 

ENTER  GAIN  OF  INTEREST  (GAIN) .TOLERANCE  (GTOL) :>-12.95, . 1 


OPEN-LOOP  (OLTF)  ROOT  LOCUS  USING  OPTION  42 


M  =  0.8 
h  *=  0  ft 


CODE 

0 -LOCUS  PT.  2-ZERO  4- I MAG I NARY  AXIS 

1-POLE  3-BREAK  PT.  5-SENS I TI VI Tr  PT. 


GAIN  OF  INTEREST  (GAIN)  -  -12.95  +OR-  GTOL  -  .1080 


SENSITIVITY  OF  INTEREST  (OLK)  «  13.46  +0R-  KTOL  -  .1039 


ROOTS  OF  INTEREST 
X  «  -.3370  1E-0  1 

X  «  -.4618  IE-8 3 

X  »  -7.9607 

X  -  -7.9607 


Y  ■  0. 

Y  »  0. 

Y  -  7.4548 

Y  «=  -7.4540 


REGION  Or  CALCULATION-REAL:  C'«  -15.0  TO  AA«  1.00 

IMAJ :  DD=  -15.0  TO  BB~  15.0 

BRANCH  STARTING  AT  (-.03328)  ♦  J(0.) 

42 

ENTER  GAIN  OF  INTEREST  (GAIN) .TOLERANCE  ( GTOL)  : >- 10 . 64 , .  1 


OPEN-LOOP  (OLTF)  ROOT  LOCUS  USING  GFTION  42  K  =  1  .0 


h  =  20K  ft 

CODE 

0-LOCUS  PT.  2-ZERO  4- 1 MAG I  NARY  AXIS 

1-POLE  3-BREAK  PT.  5-SENSITIVITY  PT. 

GAIN  OF  INTEREST  (GAIN)  *  -10.«.4  +  0R-  GTOL  *  .1008 
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1 


SENSITIVITY  OF  INTEREST  <OLK> 


9.696 


♦OR-  KTOL  « 


9 1  13E-0 


ROOTS  OF  INTEREST 
X  -  - . 24320E-0 1 

X-  . 14515E-02 
X  -  -5.4146 

X  -  -5.4146 


Y  ■  0. 

Y  «  0. 

Y  -  -5.0643 

Y  «  5.0648 


REGION  OF  CALCULATION-REAL!  CC«  -5.00  TO  AA«  1.06 

IMAJs  DD»  -3.00  TO  BB*  3.00 

BRANCH  STARTING  AT  C-. 02451)  ♦  J(0.) 

42 

ENTER  GAIN  OF  INTEREST  (GAIN) .TOLERANCE  (GTOL) s>-6.0 1 , . 1 


OPEN-LOOP  < OLTF)  ROOT  LOCUS  USING  OPTION  42 


CODE 

0-LOCUS  PT. 
1-POLE 


2- ZERO 

3- BREAK  PT. 


M  =  1.0 

h  -  4 OK  ft 

4- IMAGINAFY  ~  I  S 

5- SENSITIVITY  PT. 


GAIN  OF  INTEREST  (GAIN) 


-6.010 


♦OR-  GTOL 


.  1000 


SENSITIVITY  OF  INTEREST  (OLK) 


6 .648 


♦OR-  KTOL  =  .1106 


ROOTS  OF  INTEREST 


X  - 

— . 39390E-02 

Y  - 

. 1 1616E-0  1 

X  - 

-.39390E-02 

Y  » 

-. 116I6E-01 

X  - 

-3.5604 

Y  » 

-3.5975 

X  - 

-3.5604 

Y  ** 

3.5975 

REGION  OF  CALCULATION-REAL;  CC«  -5.00 
.  IMAJ;  DC*  -3.08 
BRANCH  STARTING  AT  (-.0037698)  ♦  J(. 012623) 
42 


TO  AA=  1.00 
TO  BB~  3.00 


ENTER  GAIN  OF  INTEREST  (GAIN)  .TOLERANCE  (GTOL)  : >- 13 .  1 1 . .  I 


OPEN-LOOP  (OLTF)  ROOT 
CODE 

0-LOCUS  PT. 

1-POLE 


LOCUS  USING  OPTION  42 


2- ZERO 

3- BREAK  PT. 


M  -  1.2 
h  =  20K  ft 

4- IMAGINAF.  i  rv\IS 

5- SENSITIVI T i  FT. 


GAIN  OF.  INTEREST  (GAIN)  *  -13.11  +OR-  GTOL  »  .1000 

SENSITIVITY  OF  INTEREST  (OLK)  ■  9.201  ♦OR-  KTOL  =  .7018E-U 


ROOTS  OF  INTEREST 

X  -  - . 33936E-0 1  Y  »  9. 

X  -  . 25730E-02  Y  ®  0. 

X  -  -5. 1699  Y  -  5.4792 

X  -  -5. 1699  Y  -  -5.4792 


REGION  OF  CALCULATION-REAL;  CC=  -5.00  TO  AA»  1.00 

IMAJ;  DD*  -3.00  TO  9B=  .  .  0 
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BRANCH  STARTING  AT  C-. 034739)  ♦  J<0.) 

42 

ENTER  GAIN  OF  INTEREST  (GAIN)  .TOLERANCE  ( GTOL)  t >-9 . 46 , .  I 


OPEN-LOOP  < OLTF)  ROOT  LOCUS  USING  OPTION  42 


CODE 

0 -LOCUS  PT. 
I -POLE 


2- 2ERO 

3- BREAK  PT. 


M  *  1 .4 
h  *  40K  ft 

4-  IMAGINARY"  axis 

5- SENSITIVITY  PT. 


GAIN  OF  INTEREST  (GAIN) 


-9.460 


♦OR-  GTOL 


SENSITIVITY  OF  INTEREST  <OLK>  ■  6.335 


♦OR-  KTOL 


.696  IE -0 


ROOTS  OF  INTEREST 

X  «  -.30 155E-03  Y  «  0. 

X  ■  -. 13 790E-01  Y  »  0. 

X  ■  -3.5589  Y  -  4.1167 

X«  -3.5589  Y  «  -4.1167 


REGION  OF  CALCULATION-REAL:  CC»  -5.00  TO  AA=  1.00 

IMAJ :  OD-  -3.00  TO  BB*  3.00 

BRANCH  STARTING  AT  (-.0155395)  ♦  J(0.) 


42 


ENTER  GAIN  OF  INTEREST  (GAIN)  , TOLERANCE  (GTOL) :>-3.49, .  1 


OPEN-LOOP  (OLTP)  ROOT  LOCUS  USING  OPTION  42 


M  «  0.3 

h  *  0  ft 


CODE 

8-LOCUS  PT.  2-ZERO 

1-POLE  3-BREAK  PT. 


4- IMAGINA=>Y  AXIS 

5- SENSITIWT-r  PT. 


GAIN  OF  INTEREST  (GAIN)  »  -3.490 


♦  OR-  GTOL  »  .  1000 


SENSITIVITY  OF  INTEREST  (OLK)  -  3.439  ♦OR-  KTOL  »  .1564 


ROOTS  OF  INTEREST 
X  -  .26956E-02 

X  »  . 26956E-02 

X  -  -3.8201 

X  -  -3.0201 


REGION  OF  CALCULATION-REAL: 

IMAJ: 

BRANCH  STARTING  AT  (.8043133)  ♦ 
42 


Y  -  . 547 1 1E-0 1 

Y  -  -.5471  IE-01 

Y  -  3.3635 

Y  »  -3.3035 


CC-  -5.00  TO  AA=  1.00 
DD-  -3.80  TO  68=  3.30 
J(. 060757) 


ENTER  GAIN  OF  INTEREST  (GAIN) , TOLERANCE  ( GTGD  : >-3 . 06 ,  .  1 


OPEN-LOOP  (OLTF)  ROOT 
CODE 

0-LOCUS  PT. 

1-POLE 


LOCUS  USING  OPTION  42 


2- ZERO 

3- BREAK  PT. 


M  =  0.4 
h  «  20K  ft 


4-  1  MAG  I  NARY  h.-'IS 

5- SENSITIVIT.  FT. 


GAIN  OF  INTEREST  (GAIN)  «  -3.060 


♦OR-  GTOL  -  .1000 


SENSITIVITY  OF  INTEREST  (OLK)  -  5.581 


♦  OR-  KTOL  =  .  1S24 


ROOTS  OF  INTEREST 
X  »  . 37562E-02 

X  -  .37562E-82 

X  -  -2.9694 

X  ■  -2.9694 


Y  -  .48607E-81 

Y  -  - . 43o07E-0 1 

Y  -  -3.2263 

Y  -  3.2268 


REGION  OF  CALCULATION-REAL:  CC-  -5.00  TO  AA-  1.08 

IMAJ:  DD-  -3.00  TO  BB*  3.0C* 

BRANCH  STARTING  AT  (.0046977)  ♦  J(. 052395) 

42 


ENTER  GAIN  OF  INTEREST  ( GAIN) .TOLERANCE  (GTOL- :>-2.37, . 1 


OPEN-LOOP  (OLTF)  ROOT  LOCUS  USING  OPTION  42 
COOE 


129 


I -POLE 


3-BREAK  PT 


5-SENSITIVITY  PT 


GAIN  OF  INTEREST  (GAIN)  ■  -2.370  ♦OR-  GTOL  -  .1668 

SENSITIVITY  OF  INTEREST  <OLK>  »  7.896  ♦OR-  KTOL  -  .3332 

ROOTS  CF  INTEREST 
X  -  .50595E-02 

X  -  .56595E-02 

X  ■  -3.9592 

X  ■  -3.9592 


REGION  OF  CALCULATION-REALS  CC-  -5.00  TO  AA»  1.00 

IMAJs  DO*  -3.00  TO  BB»  3.00 

BRANCH  STARTING  AT  (.0054698)  ♦  J< . 051752) 

42 

ENTER  GAIN  OF  INTEREST  < GAIN) .TOLERANCE  ( GTOL) s >-9 . 24 , . 1 


OPEN-LOOP  (OLTF)  ROOT  LOCUS  USING  OPTION  42 


M  *  0.6 
h  *  0  ft 


CODE 

0 -LOCUS  PT.  2-2 ERO  4- I MAG I NARY  AXIS 

1-POLE  3-BREAK  PT.  5-SENSITIVITY  PT. 


GAIN  OF  INTEREST  (GAIN)  »  -9.246  +0R-  GTOL  «  .I960 


SENSITIVITY  OF  INTEREST  (OLK)  -  10.10  +0R-  KTOL  »  .1093 


ROOTS  OF  INTEREST 
X  »  -. 10 159E-0  1 

X-  - . 10 159E-0  1 
X  -  -5.8805 

X  «  -5.8805 


Y  ■  .55097E-02 

Y  -  - . 55097E-02 

Y  -  5.3295 

Y  -  -5.3295 


• 

REGION  OF  CALCULATION-REAL:  CC«  -5.00  TO  AA=  1.00 

IMAJs  DD*  -3.00  TO  BB»  3.00 

BRANCH  STARTING  AT  (-.0098176)  ♦  J(. 0091793) 

42 

ENTER  GAIN  OF  INTEREST  ( GAIN) .TOLERANCE  ( GTOL) s >-5 . 42 , . 1 


OPEN-LOOP  (OLTF)  ROOT 
CODE 

0-LOCUS  PT. 
1-POLE 

GAIN  OF  INTEREST 


LOCUS  USING  OPTION  42 


2- ZERO 

3- BREAK  PT. 

GAIN)  «  -5.420 


M  =  0.6 
h  *=  20K  ft 


'’-IMAGINARY  AXIS 
w  SENSITIVITY  PT. 

♦OR-  GTOL  =  . 1000 


SENSITIVITY  OF  INTEREST  (OLK)  -  7.137  +0R-  KTOL  *  .1317 


ROOTS  OF  INTEREST 
X  -  - . 248O7E-02 

X  -  -.248O7E-02 

X  «  -3.9429 


.21861E-0  1 
- . 2 186 1E-0 1 
3.8222 


Y  - 

Y  « 

Y  - 


TO  68*  3.06 


IMAJ:  DD-  -3.00 
BRANCH  STARTING  AT  <-.0019434)  +  J(. 024173) 
42 


ENTER  GAIN  OF  INTEREST  (GAIN) , TOLERANCE  < GTOL) : >-2 . 88 ,  .  1 


OPEN-LOOP  (OLTF)  ROOT 
COOE 

0-LOCUS  PT. 

1-POLE 


LOCUS  USING  OPTION  42 


2- ZERO 

3- BREAK  PT. 


M  =  0.6 
h  «  40K  ft 

4-  IMAGINAR  Y  I  S 

5- SENSITIVITY  PT. 


GAIN  OF  INTERFST  (GAIN)  »  -2.880 


♦  OR-  GTOL  *  .  1000 


SENSITIVITY  OF  INTEREST  (OLK)  -  5.744 


♦OR-  KTOL  =  .20 


ROOTS  OF  INTEREST 
X  -  .35943E-02 

X  *  . 35943E-02 

X  -  -3.0001 

X  -  -3.0001 


Y  -  . 34407E-0  1 

Y  *  - . 34407E-0  1 

Y  -  -3.2901 

Y  -  3.2901 


REGION  OF  CALCULATION-REAL:  CC*  -5.00  TO  AA=  1.00 

IMAJ:  DD*  -3.00  TO  BB*  3.66 

BRANCH  STARTING  AT  (.0040104)  ♦  J< . 03834) 

42 

ENTER  GAIN  OF  INTEREST  < GAIN)  .TOLERANCE  < GTOL)  : >- 12 . 44  ,  .  1 


OPEN-LOOP  (OLTF)  ROOT  LOCUS  USING  OPTION  42 


M  *  0.8 

h  *  0  ft 


CODE 

0-LOCUS  PT.  2-ZERO  4— IMAGINARY  —'IS 

1-POLE  3-BREAK  PT.  5-SENSITIVIT v  PT. 


GAIN  OF  INTEREST  (GAIN)  *  -12.44 


♦  OR-  GTOL  -  .  1000 


SENSITIVITY  OF  INTEREST  (OLK)  *  13.14  ♦OR-  KTOL  =  .10 


ROOTS  OF  INTEREST 
X  *  -.3070 1E-0  1 

X  -  -.51049E-03 

X  *  -7.8115 

X  ■  -7.8115 


REGION  OF  CALCULATION-REAL: 

IMAJ: 

BRANCH  STARTING  AT  (-.03328)  ♦ 
42 


Y  -  0. 

Y  -  0. 

Y  -  7.5555 

Y  -  -7.5555 


CC*  -5.00  TO  AA*  1.00 

DD*  -3.08  TO  BB*  3.00 

J<0  .) 


ENTER  GAIN  OF  INTEREST  (GAIN) .TOLERANCE  < GTOL)  : >- 10 . 45 
,  .  1 


OPEN-LOOP  (OLTF)  ROOT  LOCUS  USING  OPTION  42 
CODE 

0-LOCUS  PT. 


M  *  1 .0 

h  *  2 OK  ft 


2-ZERO 


4- I MAG l NARY  A ■ IS 


REGION  OF  CALCULATION-REALS  CC-  -5.00  TO  AA-  1.00 

IMAJ:  DO  -3.00  TO  BB«  3.00 

BRANCH  STARTING  AT  (-.034739)  ♦  J<0.) 

42 

ENTER  GAIN  OF  INTEREST  (GAIN) .TOLERANCE  ( GTOL) s >-9 . 27 , . 1 


OPEN-LOOP  < OLTF)  ROOT  LOCUS  USING  OPTION  42 
CODE 

0-LOCUS  PT.  2-ZERO 

1-POLE  3-BREAK  PT. 

GAIN  OF  INTEREST  (GAIN)  ■  -9.270 

SENSITIVITY  OF  INTEREST  <OLK>  -  6.453 

1 


K  «  1.4 
h  «  40K  ft 

4-  IMAGINARY  AXIS 

5- SENSITIVITY  PT. 

♦OR-  GTOL  ■  .1000 

♦OR-  KTOL  -  . 696 1E-0 


ROOTS  OF  INTEREST 
X  ■  — .30269E-03 

X-  -. 15816E-0  1 
X  »  -3.4927 

X  -  -3.4927 


Y  -  0. 

'Luis 

Y  ■  -4.1618 


REGION  OF  CALCULATION-REAL! 

IMAJ  s 

BRANCH  STARTING  AT  (-.015595)  ♦ 


CC-  -5.00 
DO  -3.00 
J(0.) 


TO  AA»  1 .00 
TO  BB-  3.00 


K 


LV.4I 


r.  v 

!»  V 


Appendix  F 

Least  Squares  Curve  Fittin 


To  find  the  best  equation  for  the  gain  schedule  of  Kq, 
the  least  squares  method  was  used.  A  curve  must  be  fit  to 
the  data  points  in  Fig  10.  The  basic  system  equation  is  the 
starting  point: 

Ax  =  b 

Premultiplying  by  AT  gives: 

AtAx  =  ATb 

Solving  for  x,  the  result  is: 

x  =  (ATA)"1ATb 

In  this  case,  the  vector  x  contains  the  coefficients  for  the 
equation  of  the  curve  which  fits  the  points.  The  work  done 
in  Method  III,  which  was  discussed  in  the  section,  "Selec¬ 
tion  of  Gains,"  will  be  shown. 

For  this  case,  the  points  in  Fig  10  were  divided  into 
two  groups  and  a  curve  was  fit  through  each  group.  The 
lower  six  points  were  in  one  group  and  the  upper  five  were 
in  another.  For  group  1: 


133 


A 


and 


‘l  Q1 

s 

1  69 

1  Q2 

1  99 

1  q3 

1  109 

1  q4 

1  134 

1  q5 

1  246 

1  Qfi 

1  275 

S' 

II 

00 

o 

• 

CM 

1 _ 

Xq2 

3.00 

K*3 

3.18 

3.65 

5.25 

00 

CTi 

• 

in. 

In  this  case 


x 


-1 


liiiiii 

'1  69' 

r  i  i  i  i  i  i] 

'2.08' 

69  99  109  134  246  275J 

1  99 

(69  99  109  134  246  275j 

3.00 

1  109 

L 

3.18 

1  134 

► 

3.65 

1  246 

5.25 

1  275 

W  W 

5.98 

1.174 

0.01727 


The  vector  x  has  given  us  the  coefficients  for  the  equation: 


Kq  =  C  +  DQ 


We  found  that: 


Kq  =  1.174  +  0.01727Q 
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I 


I 


This  same  procedure  is  repeated  for  the  upper  five 


points.  Now: 


x 


II  1 

1 

1 

1 

1 

11535 

539 

682 

951 

982 

v-1 


'l  535' 

f  1  1  1  1  il 

8.55' 

1  539 

1535  539  682  951  9821 

10.25 

1  682 

► 

10.08 

1  951 

12.21 

1  982 

13.45 

4.839 

0.008226 


This  results  in  the  equation: 

Kq  =  4.839  +  0.008226Q 

At  the  intersection  of  the  two  lines,  one  finds  the  Q 
at  which  the  gain  schedule  will  change  from  one  equation  to 
the  other.  Setting 

1.174  +  0.1727Q  =  4.839  +  0.008226Q 

2 

one  finds  that  the  crossover  Q  is  405  lb/ft  . 

Therefore,  the  gain  schedule  for  Kq  is: 

For  Q  less  than  or  equal  to  405  lb/ft2 

Kq  =  1.174  +  0.01727Q 

And  for  Q  greater  than  405  lb/ft2 


K„  =  4.839  +  0 . 008226Q 
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